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Abs 1 ra  ■ t 


This  paper  examines  the  problem  of  building  a machine 
to  reed  uncontrolled  type  fonts  set  with  essentially  no 
space  between  letter:-,  (within  words).  The  consequence  of 
this  type  of  data,  which  represents  the  usual  format  of 
printed  text,  is  that  the  data  vectors  produced  by  the  opti- 
cal scanner  contain  multiple  letters  and/or  fragments  of 
letters  that  cannot  be  easily  separated. 

An  algorithm  based  on  a variant  of  running  cross- 
correlation  between  prototype  letters  and  successively  "win- 
dowed" fragments  of  the  sentence  is  employed.  The  algorithm 
computes  the  Euclidean  distance  between  prototypes  and  the 
sentence  fragment  in  a filtered  Fourier  domain, 

Tt  is  shown  that  appropriate  normalisation  and  windowing 
techniques  allow  perfect  recognition  of  touching  letters 
within  words.  This  occ urs  even  when  no  apriori  knowledge  of 
letter  location  within  the  word  is  available,  provided  that 
suitable  prototypes  can  be  established . 

Eulti pie  alphabet  prototypes  wore  then  built  and  used 
to  examine  widely  differing  type  fonts.  Techniques  to  set 
acceptance  thresholds  were  evaluated  and  the  behavior  of  the 
resulting  recogni tier,  system  tabulated.  A number  o'"  false 
triggers  did  occur  in  this  case  and  these  were  discussed. 
Recommendations  for  further  improvements  in  the  system  are 
suggested. 
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Back -round 

A reading:  machine  (to  read  uncontrolled,  linotype-set 
fonts)  based  on  correlation  mathematics  (Refs  3 an d h)  has 
never  been  successfully  operated,  apparently  because  of 
the  nature  of  the  noise  in  the  data  stream  that  such  systems 
are  required  to  process.  The  noise  arises  from  multiple 
mechanisms,  including  shape  variations,  and  has  never  been 
well  defined  because  the  processes  producing  it  seem  to  have 
no  useful  description.  Empirical  data  for  this  noise  does 
not  seem  to  exist  in  the  literature.  Recognition  algorithms 
are  unable  to  reliably  operate  a threshold  detection  machine, 
even  for  isolated  letters  where  the  best  prototype  has  been 
employed.  In  the  more  realistic  case  of  touching  (or  nearly 
touching)  letters,  as  set  by  linotype  or  equivalent  machines, 
an  additional  noise  consequent  to  multiple  complete  and  in- 
complete signals  in  the  observation  window  further  reduces 
the  accuracy  of  all  known  algorithms. 

In  a previous  thesis  by  Captain  Paul  Kabler  (Ref  2), 

150  alphabets  with  widely  varying  fonts  were  digitized,  and 
these  serve  as  the  basis  data  set  for  this  thesis.  The 
Fourier  transforms  of  these  vectors  were  obtained  in  the 
unpublished  work  of  Captain  Troy  Sponaugle  (Ref  5)>  and 
were  used  to  test  various  isolated  letters  in  a correlation- 
type  recognition  system  against  the  kinds  of  noise  mentioned 
above. 
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Captain  Sponaugle  chen  constructed  various  data  sen- 
tences from  the  two-space  letter  data.  These  sentences 
were  each  made  from  a distinct  alphabet  with  fixed  spacing 
between  the  letters.  Further  improvements  of  the  computer 
algorithms  enabled  reading  the  sentence  by  passing  a v/indow 
across  the  sentence  by  single  element  increments.  This 
program  and  the  data  files  are  the  basic  research  tools  for 
this  thesis. 

Object 3 vfs 

The  major  objective  of  this  thesis  is  to  solve  the 
segmentation  problem  of  touching  English  letters.  A method 
of  shaping  the  windows  to  de-emphasize  the  energy  in  the  right 
and  left  edges  of  the  letters  is  used.  This  method  eliminates 
some  of  the  noise  introduced  by  the  neighboring  letters,  and 
in  fact,  produces  more  reliable  recognition. 

The  computer  performance  is  evaluated  by  comparing 
the  computer  output  with  the  actual  content  of  the  window; 
this  enables  acceptance  thresholds  to  be  set  for  each  letter. 
Therefore,  a simple  decision  rule  can  be  implemented  in  this 
recognition  system. 


Approach 


This  section  begins  by  summarizing  the  logic  executed 
by  the  actual  pattern  recognition  algorithm.  It  then  con- 
tinues by  discussing  the  various  window  shaping  functions 
that  can  be  applied  by  the  programs.  The  types  of  window 
shaping  functions  used  includes 

Unshaped  (or  Full) 

Trapezoid 
G auosian 

This  section  concludes  with  a summary  of  the  data 

files  and  computer  programs  (Ref  5)  used  in  this  thesis. 

The  tv/o  data  files  are: 

LETTERS 

SENDATA 

The  three  computer  programs  used  ares 

KAKFLTR 

PROTOBLD 

CKPRSEN 

Figure  1 gives  a flowchart  of  the  way  that  the  programs 
and  data  files  interact.  In  the  figure,  the  circles  represent 
the  data  files  and  the  rectangular  boxes  represent  the  com- 
puter programs.  The  triangular  boxes  are  files  built  by 
the  indicated  programs  for  use  in  later  programs. 


Pattern  Recognition  Algorithm 

A brief  description  of  the  pattern  recognition  algorithm 
used  in  this  thesis  follows.  This  algorithm  is  used  in  the 
KAKFLTR  and  CKPRSEN  program,  which  is  discussed  later  in 
this  thesis. 

A window  is  brought  in  as  a digitized  32  x 32  element 
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Fig  1 . Interaction  of  Data  Files  and  Computer  Programs 


array,  consisting  of  zeros  and  ones.  This  array  is  then 
input  to  a discrete  Fourier  transform  routine,  FOURT  (Ref  1), 
which  obtains  the  two-dimensional  Fourier  transform  of  the 
input  array.  The  discrete  two-dimensional  Fourier  transform 


is  calculated  as  follows! 


fts 
L 


where  j T(k,I)  = input  array,  and  0(r,s)  = output 

array. 

The  array  that  is  output  from  the  FOURT  routine  is 
then  further  processed.  First  it  is  filtered  in  order  to 
eliminate  the  spatial  harmonics  greater  than  the  fourth 
(Ref  6).  This  allows  reduction  of  the  signal  space  to 
Ht  - space. 
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The  array  is  then  energy  normalized  to  insure  that  the 


energy  of  each  vector  is  unity.  Energy  normalization  is 
computed  as  follows: 

N(r.s)-  0Cr,s)/[||:02(r,S)]'/2 

where  0(r,s)  is  output  array  from  FOURT  routine  and  N(r,s) 
is  the  normalized  array. 

The  normalized  array  is  then  reordered  into  a one- 
dimensional array.  This  new  array  contains  the  L/C  term 
followed  by  the  real  part  of  the  first  harmonic,  the  imag- 
inary part  of  the  first  harmonic,  the  real  part  of  the 
second  harmonic,  etc.  This  is  the  same  format  used  for  all 
data  vectors. 

The  next  step  in  the  algorithm  is  to  calculate  the 

« - 

Euclidean  distance  between  the  input  and  a prototype.  This 
is  simply  done  using  the  following  equation: 

DisW=[|[V/Vt)-Pz(0]]'/2 

where  V/(i)  = one-dimensional  window  array  and  P(i)  = one- 
dimensional prototype  array.  Thus  when  the  input  array  is 
exactly  the  same  as  the  prototype  array  the  distance  is  equal 
to  zero  - or  a perfect  match. 


’v'indow  shapes 

The  three  window  functions  used  in  this  thesis  arc 
the  unshaped  (or  full),  trapezoidal  and  Gaussian  shapes. 

These  various  functions  are  applied  to  windows  in  both  the 
MAXI LTG  and  Ci/.iRSEM  programs  (see  the  1 rogram  section  below). 


i'aeh  of  those  window  fund'd  on.i  are  calculator  to 


obtain  32  normalized  values,  with  the  highest  value  being 
unity,  The  values  are  put  Into  an  array  such  that  the 
unity  values  are  centered  within  it.  Then  each  of  the  32 
columns  within  the  window  are  multiplied  by  the  cor responding 
array  value.  This  process  causes  a de-emphasis  of  the  window 
edges,  thus  reducing  the  noise  introduced  by  the  neighboring 
letters.  This  process  also  reduces  the  signal  but  not 


enough  to  significantly  change  the  recognition  distances. 

In  the  case  of  the  unshaped  (or  full)  window  function 
the  above  described  array  contains  all  ones.  This  allows 
the  window  to  bo  read  in  directly,  with  no  change  made.  Tt 
is  used  only  to  see  the  effect,  the  neighboring  letters  have 
on  recognition  distances. 

In  this  thesis  the  trapezoid  plateau  length  was  varied 
in  order  to  study  window  shape  affects.  The  various  values 
used  for  the  plateau  length  were  1,  5,  10,  1%  TO,  25,  and 
^0  columns,  Examples  of  those  window  functions  arc  given  in 
Figure  2. 

The  Gaussian  window  function  was.  obtained  using;  the 
following  equation* 

| exp[-(x-iG)/2.s2]  x=«,v^ 
[exp[-(x-i7)/2s2]  xmvv?1 


where  W(x)  is  the  Gaussian  window  array,  x is  the  column 

position  along  the  window,  and  ■ is  t Vo  standard  deviation. 
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2.  Various  Trapezoidal  Window  Shape 


tandard  deviations  used  were  values  of  2 


and  10  columns.  Examples  of  these  window  functions  are 
given  in  Figure  3. 

The  above  described  process  is  executed  in  the  sub- 
routine TMF'.iLT . This  routine  is  listed  in  the  I'.AKELTR  and 


Cr.'PRSEN  programs  (see  Appendix  C,  pg.  69?  Appendix  E,  pg.  101 
respectively) . It  should  be  noted  that  these  listings 
contain  only  the  Gaussian  logic,  but  any  type  window  can  be 
implemented  by  simply  substituting  the  new  function  logic. 


LETTERS.  This  data  file  contains  the  digitized  letters 
of  150  various  font  alphabets.  Each  of  these  stored  letters 


are  isolated  and  centered  in  a 32  x 32  element  window 


'•'if:  3 Various  Gaussian  Window  ohapcs 

This  data  file  ir.  used  by  the  riAKFLTR  program  discussed 
below.  Some  examples  of  the  data  in  this  file  are  shown  in 
Appendix  A. 

SKN'dAT A . This  file  contains  various  numbered  sentences. 
Each  sentence  is  133  columns  lone;  and  3-3  rows  high.  These 
sentences  consist  of  letters  in  a line,  in  the  manner  of  an 
ordinary  sentence  but  with  no  spacing  between  letters,  Each 
letter  within  a sentence  comes  from  a single  alphabet  so  there 
is  no  mixing  of  fonts  within  the  sentence. 

This  data  file  is  used  by  the  CMFRSEN  program  discussed 
below,  A complete  list  of  the  sentences  used,  and  their 
characteristics  are  shown  in  Appendix  15. 


Programs 
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IAKVLTE . This  i>rogran  producer,  the  Fourier  transformed 


letters  necessary  to  "build  the  prototypes . It  applies  a 
shaped  window  to  each  isolated  letter  in  the  data  file  LSTTcP.S. 
This  shaped  letter  is  then  Fourier  transformed  and  stored 
in  the  format  used  by  the  PROTCBLD  program  discussed  below. 

A complete  listing  of  this  program  is  included  in  Appendix  C, 
PROTOPLD,  This  program  uses  the  file  built  by  the 
KAKFLTR  program.  It  uses  the  shaped  Fourier  letters  to  build 
average  prototypes.  The  averaging  process  used  is  the  summing 
of  the  Fourier  components  for  each  letter  and  dividing  by 
the  number  of  alphabets  used.  The  formula  used  is  as 
follows : 


P(K)=  UCK) 

Where  L^(K)  is  the  Kx  component  of  the  Fourier  representation 
of  the  i^  alphabet-version  of  a letter.  P(K)  is  the  com- 

ponent of  the  Fourier  representation  of  the  prototype. 

This  program  also  permits  the  user  to  specify  which 
alphabets  are  to  be  used  in  the  averaging  process.  A complete 
listing  of  this  program  is  included  in  Appendix  D. 

CP'PRSFN.  This  program  contains  the  recognition  algo- 
rithm. It  requires,  as  input,  both  the  prototypes  and  a 
data  sentence.  The  sentence  is  read  by  passing  a 32  x 32  element 
window  along  the  sentence.  The  input  at  each  window  position 
is  run  through  the  algorithms  discussed  in  the  previous 
sect! ons , 

These  Euclidean  distances  are  output  from  the  OF.FRSEN 
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The  first  format  in  a list  of  the 


d 


program  In  two  formats, 
distance  each  prototype  is  from  each  window  along  the  sen- 
tence . This  allows  easy  observation  of  each  distance  be- 
tween prototypes  and  the  contents  of  the  window.  The 
second  format  is  a listing  of  the  best  five  choices  (i.e. 
lowest  distances)  for  eaeh  window  position  along  the  sentence 
There  is  an  additional  capability  of  vising  cither 
complex  or  amplitude  Fourier  spectra.  However  the  same 
Euclidean  distance  logic  is  used  for  either  spectra,  A 
complete  listing  of  the  prof-yam  is  included  in  Appendix  E, 


Renul ts 


Si  n rI  c Alphabet  I rototy nf1:: 

Two  measurements  are  used  to  evaluate  the  performance 
of  the  recognition  system.  The  first  in  the  Euclidean 
distance  to  the  correct  choice.  The  second  measurement  is 
the  figure  of  merit  (PC'”),  which  in  defined  as  the  ratio 
of  the  distances  from  the  second  choice  to  the  first  choice. 
Therefore,  the  RON!  has  a range  of  values  greater  than  unity 
as  long  as  the  first  choice  is  also  the  correct  choice. 

The  first  prototypes  wore  built  from  a single,  unshapod 
alphabet,  identified  as  alphabet  #10.  This  is  the  same 
alphabet  used  to  construct  sentences  numbered  one  to  four 
(see  Appendix  B,  pg.  26).  Thus  we  have  perfect  prototypes 
for  each  letter  in  alphabet  //10j  in  the  sense  that  the  sen- 
tences are  constructed  from  the  same  letters  as  the  proto- 
types. 

Using  the  above  prototypes,  sentence  number  one 
("ABGD2PG")  was  read  in  with  the  unshaped  (or  full)  function 
being  used.  This  allowed  a standard  for  further  comparisons. 
The  only  noise  in  the  data,  in  this  case,  is  provided  by  the 
adjacent  letters  in  the  scan  windows.  This  sentence  was 
also  read  in  using  the  various  trapezoidal-shaped  function 
described  in  the  previous  section.  The  lowest  recognition 
distances  and  highest  FOR’,  to  the  correct  choices  occurred 
when  the  trapezoid  plateau  length  was  equal  to  five  columns. 

New  prototypes  were  then  constructed  with  a trapezoidal 


tl 


shape  of  plateau  length  equal  to  five  columns  and  various 
trapezoidal  shaped  functions  were  applied  to  read  sentence 
number  one  using  these  new  prototypes.  This  time  the  best 
recognition  distances  occurred  for  trapezoidal  function 
shape  v:ith  plateau  length  of  one  column.  These  empirical 
results  suggested  a function  shape  with  a sharper  peak. 

A Gaussian  function  w as  chosen  next.  The  various 
Gaussian  prototypes  were  each  run  through  the  different 
Gaussian  window  shapes.  From  these  runs  the  best  recognition 
occurred  when  the  prototypes  and  window  function  had  a stand- 
ard deviation  of  two  columns.  The  above  results  are  summarized 
in  Tabic  I.  Note  especially  the  very  low,  virtually  zero, 
distances  to  the  correct  prototype. 

The  data  in  Appendix  F are  the  actual  computer  output 
of  the  five  lowest  distances  for  each  window  position  along 
sentence  number  one.  The  three  outputs  are  for  both  proto- 
type and  window  standard  deviations  equal  to  2,  5»  and  10 
columns.  These  data  show  how  the  recognition  distances 
depend  upon  the  standard  deviations. 

These  prototypes  are  then  run  along  the  sentences 
numbered  one  to  four  with  a Gaussian  function  shaped 
window  of  standard  deviation  of  two  columns.  Tables  II 
and  TIT  summarize  the  output  from  the  complex  and  ampli- 
tude spectra,  respectively.  The  tables  give  a listing 
of  the  lowest  distance  achieved  by  each  prototype  filter 
on  each  of  the  four  sentences.  It  shows  how  the  correct 
prototype  distances  fall  to  nearly  zero.  It  should 
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Spectrum  - Lowest  Distance  Achieved  by  Prototype 
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bo  noted  that  the  way  the  sentence  is  read  in  causes  the 
first  letter  of  each  sentence  never  to  be  centered  in  the 
window.  This  is  the  reason  that  the  letters  A,!!  and  P do 
not  drop  to  zero.  This  was  not  corrected  in  this  project 
(by  adding  leading  columns  of  zeros  to  the  sentences)  be- 
cause of  time  considerations.  The  letter  L,  in  sentence 
number  two  is  incomplete  because  when  the  sentence  was  built 
the  lower  bar  was  omitted  (see  Appendix  B,  pg.  28).  Those 
data  errors  can  be  eas.i  ly  corrected  or.ce  the  sentence 
building  program  is  restored  to  the  computer  files. 

Ku  !ti  pie  Alphabet.  Prototypes 

The  next  step  was  to  introduce  multiple  alphabet 
prototypes.  The  alphabets  used  wore  Kabler's  numbers 
10,  11,  12,  13,  14  and  15  (Kef  2).  These  are  some  of  the 
alphabets  (number  15  was  excluded)  that  were  used  to  construct 
the  sentence  data  (see  Appendix  B).  Mote  that  with  multiple 
alphabets  an  additional  source  of  noise,  really  a "shape” 
noise,  is  added  to  the  problem.  As  a consequence  one  would 
expect  the  performance  of  the  system  to  be  degraded. 

In  the  single  alphabet  prototype  case,  it  was  found 
that  the  Gaussian  shaped  window  gave  perfect  results. 
Therefore,  the  various  Gaussian  shaped  prototypes  were 
built  from  the  six  alphabets.  Bach  of  these  prototypes  were 
run  along  sentence  number  one  with  the  various  Gaussian 
functions  applied.  In  this  case,  the  best  recognition 
(i.e.  lowest  distances,  highest  FOM)  occurred  when  the 


prototypes  had  a standard  deviation  or  ten  columns  and  the 
window  had  a standard  deviation  of  five  columns. 

These  best  prototypes  were  then  -un  along  all  22 
of  the  data  sentences.  Both,  the  complex  and  amplitude 
spectra  output  was  obtained.  However,  as  expected,  the 
distances  in  this  case  did  not  drop  to  zero. 

Some  examples  of  the  raw  data  output  from  the  above 
discussion  Is  contained  in  Appendix  S.  These  data  are  the 
lowest  five  distances  for  each  window  position  along  sentences 
numbered  one  to  four.  This  Appendix  contains  the  output 
for  both  the  complex  and  amplitude  spectra. 

Acceptance  thresholds  were  then  set  for  each  letter-. 

This  would  allow  a simple  decision  scheme  to  make  the  detection 
choice.  Therefore,  for  each  window,  where  the  letter  w as 
relatively  centered  within  the  window,  the  distances  were 
recorded  for  each  of  the  26  letters.  This  was  repeated  for 
each  letter  in  every  sentence  until  five  distances  (one  for 
each  alphabet),  were  obtained  for  each  letter,  and  for  each 
spectra.  The  highest  number  for  each  letter  was  then  set  as 
the  acceptance  threshold.  Therefore,  all  letters  were  detailed 
by  the  threshold  mechanism. 

Using  the  Euclidean  distances  output  for  each  prototype 
along  the  sentences,  every  number  below  that  letter  threshold 
was  marked.  hen  both  spectra  agreed  (i.e,  triggered  at 
the  same  sentence  position)  the  program  declared  that  letter 
was  in  the  window.  But  when  this  was  done  there  were  a large 
number  of  false  triggers. 
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It  was  then  observed  that  the  letters  from  alphabet 
“11  were  determining  thresholds  in  most  of  the  cases.  These 
sentences  were  made  of  atypically  thin  letters.  Thresholds 
were  reset  omitting  this  alphabet  and  there  was  a substantial 
decrease  in  th,e  number  of  false  triggers. 

A summary  of  these  false  triggers  is  given  in  Table  IV. 
Five  letters  (I,  J,  0,  Q,  W)  were  consistently  falsely 
triggered.  Other  false  triggers  did  occur  for  the  different 
alphabets,  and  these  are  entered  in  the  table. 


C on  • 1 u • i on: 


; and  .?oc'ori'".ptvi:i  t ; »it:, 

The  single  alphabet  prototype  results  showed  that  a 
perhaps  interesting'  fragment  of  the  segmentation  problem 
can  bo  solved.  Since  the  distances  fall  to  nearly  zero 
only  on  the  correct  letters,  a simple  threshold  detector  can 
be  implemented.  The  detector  can  be  used  on  either  the  complex 
or  amplitude  spectra. 

The  multiple  alphabet  prototype  results  turned  out 
fairly  well.  The  problem  of  false  triggers  that  occurred 
could  be  attacked  in  a number  of  ways.  T’irst,  ami  most  ob- 
vious, is  based  on  the  fact  that  the  prototypes  were  an 
average  of  the  letters  in  six  alphabets.  These  alphabets 
included  an  atypically  thin  font  (alphabet  ;f‘ 11)  and  a thick 
font  (alphabet  //l4) , This  wide  variation  in  font  might  very 
well  account  for  many  of  the  false  triggers  obtained.  A 
method  of  ameliorating  this  problem  is  suggested  below. 

It  would  be  possible  to  add  special  logic  for  the 
troublesome  letters  (i.e.  T,  J,  0,  Q,  W),  An  example  of 
this  occurs  when  the  B and  I filters  trigger  at  approximately 
the  same  location.  Since  the  T filter  would  normally  trigger 
in  multiple  neighboring  windows,  the  special  logic  would 
trigger  the  B filter  only. 

Another  example  of  special  logic  could  be  used  in  the 
following  case.  The  recognition  system  triggers  B,  I,  and 
E filters,  within  a small  number  of  window  locations.  The 
additional  logic  would  use  the  average  column  lengths  of 


B,  1 and  K,  taklnr  Into  account  Uie  fact  that  an  T could  never 


Bit  between  those  letters,  an.)  override  the  1 tri  j/yri  ny, . 

There  are  other  types  of  special  lopie  that  miftht  be 
implemented,  Perhaps  lo  se  leased  on  Kns.lish  eonstr:  tuts  could 
be  Included  such  as  the  rule  that  U always  follows  a 4.  Tt 
is  also  possible  to  adapt ivoly  vary  the  window  sire  and 
ro -evaluate  the  dee i s i an  computation  when  an  umbipui ty  arises. 
However,  to  keep  the  generality  of  this,  process  i t is  to  bo 
hoped  the  amount  of  special  lo,‘  io  can  be  kept  to  a minimum. 

The  obvious  next  step  in  this  pv*ojoot  would  be  to  p,o 
through  the  ISO  alphabets  and  divide  the  fonts  Into  separate 
categories  according';  to  the  thickness  of  line  strokes.  Then 
prototypes  for  each  letter  could  be  built  in  these  categories 
and  tested  nyainst  various  sentences  within  the  category. 

New  thresholds  could  then  be  set  as  discus  ('1  previously  and 
the  resul ts  stud i ed . 


It  is  also  necessary  to  develop  a font  detector  for 
the  overall  recognition  scheme.  This  device  would  enable 


the  machine  to  decide  the  most  suitable  prototype  to  use  for 
any  Input  sentence,  Tt  is  hoped  the  number  of  prototypes 
required  would  be  considerably  less  than  the  t'tal  number 
of  fonts  rend  by  the  machine,  Kor  instance,  by  observation 
of  the  first  Fourier  transformed  window  of  the  sentence,  the 


L)d  term  could  be  used  to  decide  the  most  appropriate  line 
thickness  prototype  to  use. 


Another  recommendation  would  be  to  Investigate  the 
behavior  of  this,  reroprtltion  ystem  with  arhl  *r  \ry  uvpienoer 


1 


con- 


1 


of  letters  within  the  sentence,  since  this  thesis 
corned  with  the  noise  introduced  by  ncl.fhboriu; . letters. 
However,  the  sentence  data  used  in  this  thesis  is  of  a set 
sequence  of  letters.  For  example,  the  letter  C is  always, 
between  letters  B and  D.  Therefore  the  many  combinations 
at  different  letter  triplets  should  be  investigated. 
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MAIN(INPUT,0UTPUT,TAPE3,TAPF-) 

ALPHSTR,  ALPHSTP,LTRSTRT,ITRST0P,LTRSI7E,NSPACE,  TAPFTV, 
$ LOGICAL  AC JUST  , RANOIN  * REAL  WKATIO,HRATIO 
FNORM,  OC 

LTR, VEXPANr,XFORM, IMA GF , FLIP, IN VERSE, ST  ORE 


NSPACE=81  $ 
ALPHSTP=  ISO 
PA  NOIN= • T. 
FLTR=.F.  f 


TAPEIM=3  f 
$ L TRS' 
S LTR=.F.  % 
INVERSE1 . F . ' 


APEDUT  =<* 

RT  = 1 $ L TRSTOP  = ?F 

VEXPANO=.F . $XFOPM=.F. 
STORE=.T. 


nc=c 


100 


110 


x-v 


PROGRAM 
INTEGER 
X TAPEOUT 
INT EGER 
LOGICAL 
FORMAT  (1HT) 

LTRSIZE=3? 

ALPHSTR=1 
AO JUST  = • F. 

IMAGE*  .'. 

ENORM=  0 ? 

WRITE  1 

IF  (RANOIN)  CALL  PRFPLTK (TAPEIN) 

CALL  PREPFL  (TAPECUT) 

CALL  FORTAPF(ALPHSTR,ALPHSTP,LTRSTR  , L T^STOP , l T RS I Z E, 

X NSPACE,  AOJUST  ,WP AT  10 , HR  ATI 0,T APEIN,  PA NOIN,  TAPEOUT  , L TR, 

X VEXPA no,xform,iwage,fltr,inverse,store,fnorm,dc> 

STOP 

FNO 

SUBROUTINE  FORTAPF(ALPHSTR,ALPHSTP,LTkSTR"  ,LTRSTOP,LTR31ZE, 

X NSPACE, AOJUST, WR A TIO, HP A TIG, TAPEIN,  E A NOIN, TAPEOUT, IT R, 

V.  VE XPA NO, XFORM,  I MAG E,F LTR, INVERSE, ST  ORE, FNORM, DC) 

INTEGER  ALPHSTR, ALPHSTP,LTRSTRT,LTRSTOP,L'  RSI7E,  NSPACE,  TAPEIN, 
XTAPEOUT  $ LOGICAL  AOJUST  ,RA  NOIN  $ REAL  WRAT 1 0 ,HP  ATI  0 
LOGICAL  LTR,VEXPANC,XFORM,IMAGE,FLTP, INVERSE, STORE 

INTEGER  L ETT  ER (3 2 , 32) , NUMLTRS, PECNUM 
COMMON/LETTERS/LETTER,NUMLTRS,RECNUM 
INTEGER  RECFLTR 
REAL  FLTRLTR(10O) ,PPOTO(180) 

COMMON/PROFLTR/FLTRLTR, PROTO, RECFLTR 
COMMON/X FORMS/ FOPTLTR 
COMPLEX  FORTLTR( 64 , 6A ) 

INTEGER  ALPHNUM, LTRNUM,  WIDTH, HEIGH' , NUMROWS , NUMCOLS 

FORMAT  (1 HI ) 

FORMAT  (1H  , “ALPHA8FT  ",I3,"  COMPLETE.-) 

WRITE  1 

NUMWROS*  NSPACE  *4 

no  110  ALPHNUM=ALPHSTR,ALPHSTP 

DO  100  ltpnum=ltpstrt,ltrstop 

CALL  GETLTR (ALPHNUM, LTRNUM, TAPEIN,  RP  NO IN) 

IF  (LTR)  CALL  PRNTLTR (ALPHNUM , LTr NUM, . F T F T. > 

CALL  XFOPMIT  (LTRSI 7F, WICT H, HEIGHT , AC JUST , WRATIO , HP AT  1 0,  NSPA CE , 
X NUMRO  WS  , NUMCOLS  » ALPHNUM  , LTRNUM,  VEX  PAN  0,  XF0RM,1MAGC', INVERSE) 

IF  (XFORM) 

y CALL  PRMTCPX ( F ORT LTR, NUMROWS .NUMCOLS , 1 , NUMROWS , 1 , NUMCOLS, 

XI, ALPHNUM  , LTRNUM, 0,.T.,.T.) 

CALL  FILTER  ( FORT LTR, NUMROWS , NUMCOLS , FL TRLTR, NSPACE) 

IF  (FLTR)  CALL  PRNTFLR ( ALFHNUM, LTRNUM, FLTf LTR, NSPACE, 3,  .F.) 

IF  (ENORM.EQ.l)  CALL  ENORMFL (FLTRLTR, FLTRLTR, NSPACE, OC) 

FLTRLT  R( NSPACE *1 ) =WIDT H 
FLTRLTR( NS PACE *2) = HEIGHT 
FL TRLTR( NSPACE *3 ) =NUMFOWS 
FLTRLT R(NSP A CE *4 ) = NUMCOLS 
XEV=(ALPHMUM-1)*  2E+LTRNUM 

IF  (STORE)  CALL  WRITMS  (T  APEOUT,FLT«l  W,  NUMWRDS,  KEY  ) 

CONTINUE 

IF  (STORE)  WRITE  3,  (ALPHNUM) 

CONTINUE 

IF  ( .NOT  .RANOIN)  REWIND  TAPFIN 
IF  (STORE)  RFWINO  TAPEOUT 

RETURN 


i i ii  w r-»  i'  it  i r*  \ I *_ » 1 1 y • * « ^ f i , w s>  t • y n , , y v/iiiwmw|  ir^irtviww  i 

COMPLEX  TEMP(NROWS,NCOLS) ,TA HP (NUMRD WS , NUMCOLS) 

INTEGER  NROWS,NCOLS, NUMROWS, NUMCOLS 
00  100  NC0l=l, NUMCOLS 
00  100  NROW=i, NUMROWS 
TAMP(NRDW,NCOL)=TEMP(NROW,NCOL) 

100  CONTINUE 
RETURN 
END 

SUBROUTINE  ARRANGE  (TEMP,  NROWS,NCOLS, TA MP, NUMROWS , NUMCOL S > 
INTEGER  TEMP(NROWSfNCOLS) ,TAMP(NUMROWS, NU^COLS) 

INTEGER  NPOWS,NCOLS, NUMROWS, NUMCOLS 
00  100  NCOLsl, NUMCOLS 
00  100  NR0W=1, NUHROWS 
T AMP <NRDW,NCOL)=TEMP (NROW, NCOL) 

100  CONTINUE 
RETURN 
ENO 

SUBROUTINE  ENOP.MFL  (FLTRLTR  »NORMLTR,NSPACE,DC) 

REAL  FLTRLTR (NSPACE) , NOE  MLTR (NSPACE) 

INTEGER  NSPACE, OC 

INTEGER  START, NVECTOR  $ REAL  SUMSORS 

START= 2 

IF  (DC.EQ.l)  START=1 
SUMSQRS=0 .0 

00  100  NVECTOR=ST ART , NSF ACE 
100  SUMS OR S= SUMS OR S* FLTRLTR  < NVECTOR)  **2 

SUMSQRS=SUMSQRS**  .b 
IF  (SUMSQRS.EQ.O .0)  SUMSORS=l.Q 
00  110  NVECTOR=START, NSPACE 
110  N0RMLTR(  NVECTOR) =FLTRLTR (NVECTOR) /SUMS1RS 
IF  (OC.EO.O)  NORMLTP (1) =1 
RETURN 
ENO 

SUBROUTINE  FILTER  ( FORT LTR , NUMROWS , NJMC )LS, FLTRLTR, NSPACE) 
INTEGER  NUMROWS, NUMCOLS, NSPACE 
COMPLEX  FORTLTRfNUMROWS, NUMCOLS) 

REAL  FLTRLTR (NSPACE) 

INTEGEP  NHRMNCS,NUMWROS, FIGHT, LEFT, DOW  4,Nf OW,NCOL 

NHRMNC S=  NSPACE*"  .c 

RIGHT=  NHRMNCS/2+1 

LEFTsNUMCOLS-RIGHT+2 

DOWNsRIGHT 

ELTRLTR(1)=REAL(F0RTLTR(1,1)  ) 

IF  (NHRMNCS • EO  *1 ) RETURN 
NUMWROS=0 

DO  100  NC0L=2, RIGHT 
NUMWRDSsNUMWR0S*2 

FLTRLTR(NUMWROS) =REAl ( FCRTLTR (1,  NCOL  ) ) 

100  FLTRLT  R( NUMWRDS  +1  ) =AIM AG ( FORTLTR ( 1 , NOOL ) ) 

00  220  NROW® 2, DOWN 
DO  210  NCOL=LEFT , NUMCOLS 
NUMWROS=  NUMWRDS*  2 

ELTRLTR(NUMWRDS  )=PEAL(FORTLTR(NROW,NCOL) ) 

FLTRLT R(  NUMWRDS* 1)  =AIMAG  (FORTLTP  ( NROW,  NCOL  ) ) 

210  CONTINUE 

IX)  200  NC0L=1, RIGHT 
NUMWROS= NUMWRDS* 2 

ELTRLTRf NUMWRDS  ) =PEAL ( FORTLTR (NROW ,NCOL ) ) 

FLTRLT R(NUMWROS*l)=AI MAG (FORTLTR( NROW,  NC Ot ) ) 

200  CONTINUE 
2 20  CONTINUE 
RETURN 
ENO 

SUBROUTINE  GET LTR  ( A LPHNUM,LTRNUM, TARE MUM, RANDOM) 

INTEGER  ALPHNUM, LTRNUM,T AFENUM 

A WJfcA  1 , illMi  MM!  ■ .r  rnmnmm  ' wii— — ««» 


INTEGER  LETTER (3  2,  32) ,NUNLTRS,RECNUM 
COMMON /l ETTERS/LFTTER,NUMLTFS,RECNUM 

INTEGER  RECRQO, RECSKP, MASKER  iLTRSTTl  , I , J , FROM, NC OL , AL°HPAC ( P 3 2) 
X » TEMP 

COMMON /l 0 A DMS/ALFHPAC 
DATA  RECNUM/i/ 

IF  (.Nor  .RANDOM)  GO  TO  7 
LTRSTRT*  (l  TRNUM-  1 ) * 32*1 
KEY=(AlPHNUM-i)*?f  U TRNl'M 

CALL  READMS(TAPENt)M,  ALPHPAC(LTRSTPT)  ,3?,Km 
GO  TO  or 
7 CONTINUE 

RECRQD  =8  L RHNUM 

IF  ( (RECRQO-RECNUM) .EO.-l)  GO  TO  60 
RECSKP=I A PS (RECRQO-RECNUM) 

IF  (RECROD-RECNUM)  10,3r,20 
10  00  15  1=1 » RECS  KP 

15  BACKSPACE  TAPENUM 

GO  TO  30 

20  00  25  1=  l.RFCSKP 

BUFFER  IN  (TAPENUM, 1)  (A  LPHPAC  <1 ) , AL  PHPAC  (8  32) ) 

CALL  CHCKTPE ( A LPHNUM,LTRNUM»  TAPENUM) 

25  CONTINUE 

30  BUFFER  IN  (TAPENUM, 1)  (ALPHPAC(i)  , AL  PHPAC  (832)  ) 

RECNUM=RECRQOM 

CALL  CHDKTPE(ALPHNUM,LTRNUM, TAPENUM) 

60  LTRSTRT* ( LTRNUM-1 ) * 32 

MASKER=MASK(1) 

00  80  NR0W=1,32 
TENP*A  LPHPAC (L  TRSTRT ♦NRC W) 

00  70  NCOL  = 1 ,32 
IF  (TEMP. AND. MASKER)  65,68 
65  LETTER (NROW, NCOL ) =1 

GO  TO  6? 

68  LETTER  (NROW,  NCOL)  = 0 

69  TEMP*SHIFT  (TEMP,  1) 

7 0 CONTINUE 

8 0 CONTINUE 

RETURN 
FNO 

SUBROUTINE  PREPLTR  (LTRTA PE) 

INTEGER  LTRTAPE 
INTEGER  LINOEX (3901) 

CALL  OPENMS (LTRTAPE, LINDE X,  3901,0) 

RETURN 
FNO 

SU3R0UTINF  PREPFL < FLR1 AE ) 

INTEGER  FIRTAP 
INTEGER  MTNDEX  (3‘  Cl) 

CALL  OPENMS (FLRT A P, Ml NOE  X,3S 01,0) 

RETURN 
END 

SU9R0UT1  NF'  REPLArE(TRFLTR,EJRTLT^,N'lM^  )WS , NUMCOLS) 

INTEGER  NUMROWS,NUMCOLS 
DIMENSION  TRPLTP (32,32) 

COMPLEX  FDRTLTR( WUMRCWS, NUMCOLS) 

00  100  NR0W=1,NUMPQWS 
00  100  NCOL= 1, NUMCOLS 
100  FORTLTR(NRON, NCOL) *TRPLTR (NROW,  NCOL) 

RETURN 
END 

SUBROUTINE  SIZEIT  (LTRSI7E,M10TH, HEIGHT) 

INTEGER  LTRSI7E, WIDTH, HEIGHT 
COMMON /LFTTERS/LETTFR,NUMLTRS,RfcCNU« 

INTEGER  LETTER (32,32) , NUMLTRS .RECNUM 
TNTFGFP  Nrnj  .NRnH.StRICOJ  .TFMP.STRT?nW 


00  100  NCOL  a 1»  3? 

00  100  NROW* 1*  32 

IF  (LETTER (NROW, NCOL) .EG. 1)  GOTO  110 

CONTINUE 

STRTCOL*NCOL 

DO  120  NCOL*l»  3? 

TENP»3  3-NCOL 
00  120  NROH«l,32 

IF  (LFTT  FR  (NROW,  TEMP). EG. 1)  GO  TO  130 
CONTINUE 

WIDTHsTFMP-STRTCCL  +1 

IF  (WIOTH.GE.l)  GO  TO  139 

WIOTHa HFIGHT=C 

RETURN 

CONTINUE 

00  140  NPOW* 1 y 32 
00  140  NC0Lal,LTRSI7E 
IF(LETTER(NROW,NCOL)  .EO.l)  GO  TO  1>1 
CONTINUE 

strtrow=nrow 

00  160  NC0M=1,LTKSI7E 
TFMPaLTRST7E+l-NRCW 

00  160  NCCL  = l,L'rFSI7e 

T (LETTFMTEMP,NCOL)  .EO.l)  GJ  T .)  1 
CONTINUE 

HEIGHT  =TFMP-?T9TP0W  + 1 
: ETURN 

1 NO 

SUBROUTINE  CPXREAL  (REAL  TRY,  CPYA^ Y , NUMR  TW^,NUHr«JLS) 

' FAL  RFALARY  (NUM»  OWS,NUM.CLS)  $ INTEGER  NUMROWS , NUNS 0L6 

COMPLEX  CPXARY (NUMROWS,NUMCOLS) 

00  100  NROHs If  NUMROHS 
00  100  NCOLa 1, NUHCOLS 
REALARY(NROW,NCOL) =CPXAKY  (NROW, NCOL) 

RETURN 

END 

SUBROUTINE  FRMTNUM ( FMT , A , 0, C ) 

INTEGER  FMT , A , B, C 

GO  TO  (10, 20, 30,40,50) , FMT 

FORMAT  <7"  FOUR  IFF  TRANSFORM  OF  ALPHABFT  ",I3,M  L F T T F R **  ,12,7) 
WRITE  1, (A, B) 

RETURN 

FORMAT (7  " FOURIER  TRANSFORM  OF  SENTENCE  ”,I2,"  STARTING  AT  COLUMN 
X **»  13,  7) 

WRITE  2, (A, 9) 

RETURN 

FORMAT  (/,**  FOURIER  XFOFM  OF  ALPHABET  ”,I?,"  LETTER  **,12, 

X - FILTERED  TO  M , 1 3 , **  SFACE**,7) 

WRITE  3, (A, B,C) 

RETURN 

FORMAT  <“  ALPHABFT  “ , 1 3 , “ LETTER  **,12,7) 

WRITE  i, (A,B> 

RETURN 

CONTINUE 

RETURN 

FNO 

'"BVOUTI  N:  GLPRT  ( MjN,  AfSTORt  , Nt  G,  Nei  F l,SK'  °AG:  ) 

HI  MENS  ION  PEN ( ?C  * ) , A (M , K ) , o TORE ( N,  ) 
f'TNtN^IPW  FMT1  (1)  , FwT  ? ( ?)  ,FMT3(2) 

LOGICAL  SKI  PAGE 

TM~  HA  I N PROGRAM  *UST  SUFFL>  FMT  1 , T2  , 1 TMT3 

FOR  EXAHPLF  IF  MS32,  USF  THt  FOLLOWING  PA  A CARPS  ” 1A]N  F C - • “ . 
DATA  F mt  1/?0H{  1H  , 3212)  / 


PRT  A F l,T  ? , F MT  T 71  H(1H  , 3 2A 1 ) , 1 OH  ( 1H*  , 3 “*  A 1 > / 
iF  M IS  LT  *3,  THt  FOLLOWING  CARO  MAY  IE  USlO. 


t.G.  t*  32  IS  UStP. 


oooonoooc 


UA  » A I 'M  * / Cl/MIAM  ' • 

NEG  = 0 IS  USEO  TO  PRESENT  THE  MAX  VAL'JE  AS  BLANK  (WHITE)  ? 
WHILE  NEG  NOT  = 0 IS  TO  PRESENT  THE  MAX  VALUE  AS  BLACK. 

THE  CHOI SE  USUALLY  OEPENOS  UPON  THE  OPERATOR, S PREFERANC E 
RECOMMENO  NEG  BE  CHOSEN  SUCH  THAR  THE  BACKGROUND  (HOST 
OF  THE  FI  EL  0)  IS  BLANK.  IN  OTHER  WO*DS,  IF  THE  PICTURE  IS 
HIGH  VALUES  ON  A LOW  VALUE  (BLACK)  BACKGROUND,  RECOMMEND 
NEG  NOT  =0.  THIS  WILL  CAUSE  THE  HIGH  VALUES  TO  BE  BLACK 
AND  THE  BACKGROUND  TO  BE  BLANK. 

INTEGER  NEG , NFORM 
INTEGER  STORE, OEN 
INTEGER  OUMMY 


10 

20 


10  0 


140 

130 


DATA  FMT 1/10H( 
DATA  FMT2/20H( 
HAT  A FMT3/20H( 
FORMAT (1H  ,=(“ 
FORMAT  (1H0) 
OATAtOEN  (1,  J) , J= 
t)ATA(OEN(?,J),J= 
OATA(DEN(3,J),J= 
OATA(DEN(4,J),J= 
OATA(DEN(E,J),J= 
OATA(DEN(6,J),J= 
DATA  (D  EN  (7  , J ) » J= 
OAT  A(DEN(8,J),J= 
OAT  A <0  EN  (°  , J) , J= 


DATA  (OEN  (it, 
DATA  (HEN  <lf  , 


1H  ,6412)/ 

1H  ,M#", = <A1) 
1H*,1X  ,=  ( Al) ) 
#”) ) 


'#*•) 


,7>/7 (1H 
, 7) /1H— , 
,7 ) / 1H= , 
,7) /1H+, 
,7)/lH)  , 
,7) /1H1 , 
,7)/lH7, 
,7) /1HX, 
, 7 ) /lHA , 


) / 
6 ( 1H 
6 (1H 
6(1H 
6(1H 
6 (1H 
6 ( 1H 
6 (1 H 
6 ( 1 H 


)/ 

)/ 

)/ 

)/ 

)/ 

)/ 

)/ 

)/ 


J) 

,J 

= 1 

, 7) / 1FM 

, 6 ( 1H  )/ 

J) 

, J 

=1 

,7 ) /1HC 

»1H- 

,5  (1H 

) 

/ 

J) 

,J 

= 1 

,7) / 1H C 

,1M  = 

(1M 

) 

/ 

J) 

,J 

= 1 

) / lML 

,1H  + 

,!  (1H 

) 

/ 

J) 

= 1 

,7 ) /1HO 

>1M> 

(1H 

) 

/ 

J) 

t J 

= 1 

,'  ) /ihl 

, 1H  + 

,1H., 

( 

1H 

> / 

J) 

♦J 

-1 

,7 ) /IHO 

flH  + 

,1H.  , 

1J 

= , 

T ( 1 1 

) / 

J) 

1 J 

= 1 

, J/iHC 

, 1HX 

,1H*  , 

1H 

» ) 

1H-, 

2 

(1H 

) / 

J) 

> vj 

= 1 

,7)  /1HC 

,1,JX 

* lH  . i 

1H 

t-4  . 

1HC, 

2 

( 1H 

) / 

J) 

» J 

= 1 

, )/lHC 

, 1HX 

,1H.  , 

i j 

H , 

1HP  , 

2 

(1H 

) / 

J) 

f J 

= 1 

,7 ) /1HO 

,1HX 

,1H.  , 

IP 

lH", 

1 

HV, 

1HA 

OMAX=A  (1,1) 

OMIN=A  (1,1) 

OO  100  1*1, M 
00  100  J*i,N 
0MAX=AMAX1(DMAX, A (I,J)) 

0MIN*AMIN1 (OMIN,A (I, J) ) 

CONTINUE 
WRITE  20 

IF  (SKIPAGE)  PRINT  1 
PRINT  4 

IF  (NFORM  .EQ.  0)  GO  TO  200 
WRITE  1C  , (N  + 2) 

00  120  1*1, M 
L = 0 

OO  130  J=1,N 

If (NEG .EQ • 0)  OUMMY  = 20.  - 19.* (A (I,  J) -OMIN) /(DMAX-DMIN) 
IF(NEG.NE.O)  OUMMY  * 1.  ♦ 19  .*  (A  (I  , J )-DMIN)  / (DMAX-DMIN) 

L = MA  XO  ( DUMMY,  L ) 

IF (OUMMY. GT. 20)  DUMMY=20 
IF ( OUMMY . LT . 1)  OUMMY=l 
00  140  K=1 ,7 
STORE (J, K) =DEN< OUMMY, K) 

CONTINUE 

PRINT  FMT 2,N,(ST0RE(J,1) ,J=1,N) 

If  ( L .LE,  10  ) GO  TO  120 
PRINT  FMT3,N,(STORE(J,2) , J=1,N) 

If  ( L .LE.  14  ) GO  TO  120 
PRINT  FMT3,N,CSTORE(J,3) , J*1,N) 

« i l -ic-  4«  > cn  rn  4»  ■ 


PRINT  FMT3,N,<ST0RE(J,A)  ,J=1,N) 

2F  ( L • L t • lb  ) GO  70  32 I 
r-KINT  FKT 3 , N , ( ST  Pc  E ( J , : ) , J=  1 , i) 
F ( L .LE.  19  ) GO  TC  1 2 L 
PRINT  r>T^fN,(q^  E(J,  6>  , J=1,N) 
PRINT  F-T3,N, (STORE <J,7  > ,J=1,N) 
CONTINUE 
WRITE  1 , (N  + 2) 
r-0  TO  22  R 
n0  223  1=1, m 
^0  21*  J = 1 , N 


' F (NEO.FO.P)  0UMM  Y = 9C  .*98.  ( A (I  , J ) - P M IN  ) / ( DMA  X -DM’  M) 

IF  (NEG.NE.O)  OUMM Y=  1 .+98  (A (I, J) -OMIN) / (DMAX-DMTN) 

TF  (DUMMY  .GT.  99)  0UMMY=99 
IF  (OUMMY  .LT.  1 ) DUMMY =1 
STORE(  J,  1)  =DUMMY 
PRINT  FMTl, (STOREt J,l) , J=1,N) 

CONTINUE 
PRINT  3 
FORMAT (1H1) 

FORMAT  (1H  ) 

FORMAT  (1HT) 

RETURN 

END 

SUBROUTINE  I VRTCPX ( LETTER ,NUM90WS , NUMCOLS ) 

COMPLEX  LETTER (NUMROWS, NUMCOLS) 

INTEGER  NUMROHS, NUMCOLS 
REAL  TEMP1.TEMP2 
00  100  NROW= 1» NUMROWS 
00  100  NC0L=1, NUMCOLS 
J=NR0W*NC0L 

IF  (MOD( J,2) .NE.l)  GO  TO  100 
TEMP1=P.EAL(LETTER(NROW,NCOL)  )M-i.O) 

TEMP2= AIMAG( LETTER (NROW,NCOL) )*(-!. 0) 

LETTER (NR0W,NC0L)=CMPLX(TEMP1,TEMP2) 

CONTINUE 

RETURN 

ENO 

SUBROUTINE  PRNTCPX ( ARRAY , ROWSIN, COLSIN , STPTROW , ST0PR9W,  STRT CDL, 

: STOPCOLf FMTfA,B,C,SEQHRMC,SKIPAGE) 

INTEGER  ROWSINfCOLSIN,STRTROW,STOPR:}W,STR~COLfSTOPC:)L  ,’MT  fA,9,G 

OMPLEX  ARRAY  (ROWS IN, COLS 

LOGICAL  S r-  0 w T M C , S X I c A G E 

’iTEGl  NS1  r,  OWSTOP*COLSTRT»COLCTO®,T  MP,NPAG  fMSPLI7S,H<  r. 


’ lT  EG E W S 1 r , WC7CF,CCL  T 

•<ttN0Ey  (;  C ) , RINDEX  <5  ) 

FORMAT  (1WT) 

FORMAT  (1  MM 


10 

FORMAT  </, 

>v,  =(2X#I3,3X>) 

20 

Mri'u 

,I2,?X,=  (F~  .3,  IX)  ) 

30 

FORMAT  (1  H 

,rX,  = ( F 7 . 3 , 1 X ) ) 

AO 

FORMAT ( /) 

50 

PoRMAT ( 1 H 

, “ARRAY  PPINTED  IS  PR01 

r w ",TT, 

ICLM 

% I3f **  TO  " 

,13) 

70 

FORMAT  (1H 

»“  INDEXING  REQUESTED  IS 

STRAIGHT 

SEQUENTIAL. 

”//) 

SO 

FORMAT  (1H 

,**  INDFXING  REQUESTED  IS 

HARMONIC 

NUMBER.”//) 

WRITE  1 
NUMCOLS* 16 

IF  C.NOT.SEQHRMC)  GO  TO  120 
DO  100  NCOL  = STRTCOL f STOPCOL 
CINOEX (NCOL) =NCOL 
00  110  NROW=STRTPOW  fSTOPROW 
RINDEX (NROW) =NROW  . 

GO  TO  17  0 
CONTINUE 

TEMP=STOPPOW-STPTROWp1 


130 


140 


150 


ISO 

170 


100 


190 

200 


$ 


210 


HO  130  NROW*STRTROW,HALF 
RINOEX (NROW) *T EMP 
TEHP«T EMP-2 
CONTINUE 

HALF*HALF*1  S TEMP*TEHP+2 
OO  140  NROW=HALF,STOPROW 
TENP=TENP«-2 
RINOEX (NROH)aTEMP 
CONTINUE 

TEMP»STOPCOL-STRTCOL*i 

HALF*TENP/2*STRTC0L 

DO  150  NCOL* ST RTC OL , HALF 

CINOEX (NCOL ) *T  EMP 

TEMP*TEHP-2 

CONTINUE 

HALF*HALF*i  } TENP*TEMP>2 

00  160  NCOL=HALF,STOPCOL 

TEHP*T E MP+2 

CINOEX (NCOL) *TEMP 

CONTINUE 

CONTINUE 

NSPLI7S=  ( CT  OPCOL  -STRTl  01 ) /NUMCOLS+1 

COLSTRT=STR7  00 L 

00  22u  M PA GE=1»N SPLITS 

COLSTOPsCOLSTKT+MJNCOLS-l 

IE  (NGPLITS.EQ.NFAGF)  CC l STOP=ST JPCDl 

TEMPsC  OL  ETOP-C  Ol^TT  T + l 

IF  (SKTPAGF)  WRITE  r 

«<■  ITE  E ( STFTROW,STOPROW,COLSTRT,CC  l .'TOP' 

if  (StOHR^D  ir  ,ico 

WRITE  "P 

GO  TO  2PP 

WRITE  PC 

CONTINUE 

CALL  FRMTNUM(FMT,A,p,C) 

WRITE  10, (NUMCOLS, (CINOEX (NCOL) t NCOL =COLSTRT,COLST OP) ) 

00  210  NROW*ST  RT  ROW, ST OPROW 

WRITE  20, (RINOEX (NROW) ,TEMP, (REAL (ARRAY (NPOW,NCOL) ) , NCOL=COL S T RT , 

X COLSTOP)) 

WRITE  30,  (TEMP,  ( AI  NAG  (ARRAY  (NROW,  NCOL)  ) , NCOL*  COLSTPT  , 2 OLST  0 »>  ) 
CONTINUE 

COLSTRT=COLSTOP*l 

WRITE  40 

CONTINUE 

RETURN 

ENO 

SUBROUTINE  PRNTE XP ( FORTLTR, NUNROWS, NUMCOL* , ALPHNUM , LTRN UN , 

X WRATIO,HRATIO,WIOTH, HEIGHT) 

COMPLEX  FORTLTR( NUNROWS, NUHCOLS) 

INTEGER  ALPHNUH,LTRNUM, WIDTH, HEIGHT,  NUMCOL S, NUMROWS 
REAL  WRATIO,  HRATIO 

FORMAT  (1H0, "EXPANSION  OF  ALPHABET  “13”  LETTER  "I2/1H  “WIOTH  RATIO* 
X-F6.2”  HEIGHT  R A TIO=“F6 • 2, / , 1H  “LETTER  SIZE  IS  ”12“  COLUMNS  WIDE 
XBY  “ 12“  ROWS  HIGH“/1H  , “WINDOW  SI*E  IS  “12“  COLUMNS  WIDE  BY  “12 
X*  ROWS  HI GH“/// 

X , * ( 1H  ,*(Ii>/>) 

FORMAT  (1HT) 

WRITE  2 

WRITE  1,  ( ALPHNUM, LTrNUM,WRAT ID, HRATIO, 
XWIDTH,  HEIGHT, NUMCOLS, NUNROWS, NUMROWS,  ( NUMCOLS,  (IFIX( PEAL  (FORUTR 
X (NROW, NDOL) ) ) , NCOL*l , NUNCOLS ) , NROW=l , NUNROWS)  ) 

RETURN 

ENO 

SUBROUTINE  PRNTFLR  (ALPHNUM, LTRNUM, ELT RLTR , NSPACE ,FMT , S KIPA G* ) 
INTEGER  ALPHNUM.LTRNUM, NSPACE, FMT  I LOGICAL  SKIPASE 

pf*i  fiTBiTDfygpirr.  — 


i 


' N'TEGE  r N^RMNC*, ZEROES, CCWN,INPEX,  <1W,  ' TART  ,STO° 

£ F4 L t^PTV(?) 

INTcGc R TNPFX(17) 

1 FORMAT <1H') 

5 RRRMAT  (1 w i ) 

10  FORMAT  </,  X * 1 3 ( X,I2,»X)> 

2 0 FORMAT  (1  “-,12,  , 13  (F9.‘  , IX)  ) 

33  FORMATfiM  X,  13  ( F«  .*>  , IX  > ) 

<*a  FORMAT  (/) 

A 0 FORMAT  </"  NSPACF  = ",  13, “ WHICH  TS  I r'C1  (NOT  • N Tj[  HJ**  ) r 
X OR  AL0HS  RrT  ",I3,"  LETTFR  M2,/l 
WRIT E 1 

F (MOr(*J?P.\CF  , 2)  .EO.l)  GC  ^0  9 • 
w IT£  • * , (NFPACt  , ALFHNUFjl  TFN'IM) 

RFTURN 

9 0 CONTINUc 

IF  (SKI»AGE)  WRITE  E 

OATA  (INDFX(I) ,1  = 1 , 13 > /1 3, 11 ,9, 7 , 5 , 3 , 1 , 3 , 1 ,7,9,11,13/ 

DATA  (EM»TY(I)  ,1=1,7) /7*C.n/ 

CALL  FRMTNUM(FMT , A LPHNUW , LTRNUM, NSPA CE> 

nhrmncs=nspace»*  .F 

START=  (13-NHRMNCS)  /2*1 

ST0P*START*NHRMNCS-1 

WRITE  1C,(IN0EX(I> , I =ST ART ,STOP) 

P0WN=NHRMNCS/2U 
N7ER0ES=  NHRMNCS/ 2 
INDXsl 

WRITE  20, (INOX  , (FMPTY (1 > ,I=l,N7EROFS> , FlTRLTR(l), 

X ( FLTRLT  R ( N VECTOR ) ,NVEC70R=2, NHRMNCS ,2) ) 

N7ER0ES=N7ER0ES*i 

WRITE  30, ( (EMPTY(I) , 1= 1 ,NZER0E S)  , (FLTRLTR (NVECTOR) , 

X NVECTDR=3, NHRMNCS, 2) > 

START=NWRMNCSM 
no  100  NROW  = 2, DOWN 

INOX  = INDX  ♦ 2 f STCP  = ST ART  + NHRMNCS*  2-1 

WRITE  20, (INOX, (FLTRLTR (NVECTOR) ,NVECTOR  = ST ART, STOP, 2)) 

ST  ART=ST  APT  *1 

WRITE  30, ( FLTRLTR(NVECTOR) , NVECTOR=START, STOP, 2) 

STARTsST OP*l 
100  CONTINUF 
WRITE  A 9 
RETURN 
END 

SUBROUTINE  PRNTI MG ( FORTLT R, NUMROWS , NUMCOLS , ALPHNUM , LTRNUH) 

COMPLEX  FORTLTR(NUMPOWS,NUMCOLS)  t INTEGER  NUMROWS, NUMCOLS 
INTEGER  ALPHNUM, LTRNUM 
REAL  STORE  (6 A,  7) 

REAL  BRIGHT  (6A,6M 

30  FORMAT (1H  “FOURIER  TRANSFORM  OF  ALPHABET  "13"  LETTER  "12"  CCNVERTF 
xrr  TO  BRIGHTNESS  LEVELS"/, 1H  "ROWS*"!  3"  COLUMNS="I3) 

DO  92  NR OW=i, NUMROWS 
DO  92  NCO  L= 1 , NUMCOLS 

92  BRIGHT (NROW , NCOL ) = (REAL ( FORT LTR ( NROW , NCOL ) ) '+?+ AIM AG ( c J RTLV P ( NROW. 

X NCOL)  )'  *2)**. 5 

CALL  ARRANGE  (BRIGHT  , £.it  ,FA  , 8k  IGHT  , N’HSCMS  , t UMC0LS) 

CALL  SLRPT ( NUMRC WS, NUMCCLS, BRIGHT, ST  07 H, 1,1, .F.) 
wr  ITE  3j  , ( A L P H NU V , L T R NUF  , NUMROWS,  N J ■■  ■ S o L ^ ) 

return 

ENO 

U9POUTT  Nr  PRNTLTR  (ALPHNUM, LTRNUM^l^oiT  ^RINTIT  ,ALRH><P,S«?;  » GD 
lNTEGFn  A l PHNUM,  LTF'NUM 
LOGICAL  niiMF  JT,OFI^TIT,ALFHSKP,S<I',’ 

INTEGER  LFTTER  (3?,  32)  ,NUMLTRS,RE3N')h 
COMMON /LF*tfrs/LETYFR,NI'MLT  h,RF0NU‘ 

NTEGf  A n JCCL  ,NU  ■‘ROW  , NUMCOL  ,COLSl  Tt  , N JHP  * Gt  ,DL0AL°H,  TTlL'(  .) 
TUTS££&  


* » ^ y 


s i i \ C.  w * i.»v  9 

1 FORMAT  (1H-) 

9®  FORMAT  ( 1 H ♦ , 3(21X,I3,6X.I2,7X) ) 

99  FORMAT  (/) 

100  FORMAT (/  ,lX,3(8X,3Al,lX,8Al,lX,3X,Al,lX,f Ai  , IX ,2X, IX, 3A 1, 3< > ) 

101  FORMAT (iH  ,3(7X, 33Al,2X) ) 

102  FORMAT (1H*,3(7X,32A1,1X,I2>  ) 

103  FORMAT  ( /, IX , 3 (7X , 16 ( IX , 1 1) , 3X > ) 

104  FORMAT  (1H1) 

106  FORMAT  ( IX  , 3 (7 X,  4 ( IX,  A 1)  , 12(  IX,  1 1)  , 3X)  ) 

107  FORMAT  ( 1M», 3 ( 7X , 33A1 , 2X ) ) 

OATA  NUMLTRS,NUMFAGE/(l,0/ 

WRITE  1 

IF  (NUMLTRS.EQ.O)  OLOALPH=ALPHNUN 
IF  (SKIpAGE ) NUMPAGE=0 

IF  ( (NUMLTRS.EQ. 0) .ANO. (NUMPAGE.EQ.D) .OR  .SKIPAGE)  WRITE  10 
IF  (PR INTIT ) GO  TO  105 

IF  (OUMPIT .OR.  (ALPHSKP.ANO. (OLOALPH.NE.ALPHNUM) ) ) 120,105 

105  NUMLTRS=NUMLTRS  ♦ 1 
AOJCOL*33MNUMLTPS-l>  +1 
COL  SIZE*  A OJCOL  *31 

00  117  NUMROW* 1, 32 
J=  0 

00  116  NUMCOL=AOJCOL,COLSIZE 
U=J>1 

IF  (LETTER  (NUMROW,  J ).EQ.O)  115,113 

110  KFEPER  ( NUMROW, NUMCOL)=-#“ 

KEEPER  1 ( NUMROW  ,NUMCOL ) =**0,‘ 

GO  TO  116 

115  KEEPER(NUMROW,NUMCOL)=KEEPEU(NUMROW,NUMCOL)=,‘  *• 

116  CONTINUE 

117  CONTINUE 

T ITLE ( NUMLTRS*  2*6  3) =LTRNUM 
TITLE(NJMITRS* 2*62)=ALPHNUM 
IF  (PRINTIT)  GO  TO  120 
IF  (NUMLTRS.LT. 3)  RETURN 
120  IF  (NUMLTRS.EQ.O)  RETURN 
COLSIZF*NUMLTRS*33 
00  140  1=1,32 
00  130  J = 1,C0LSI7E,33 
130  KEEPERKI,  J)  =“0" 

HC  135  J=72,C0LST7E,33 
13i:  KEEPERKI,  J)  ="0" 

140  CONTINUE 

00  1 4C  1 = 1,32,31 
r0  145  J = 1 » COLSI  7C' 

145  KFEPER  J ( I , J)  ="0'* 

00  16*  r=2 ,32,2 

00  lo1  J = 3?,C0LS17c,33 

KEEPER  (T  ,J)  ="  " 

165  KEEPERKI,  J)  =1 
00  170  1=1,32,2 
00  170  J = 33 , COLS  1 7E , 3 3 
170  KEEPER  (I,J)=KEEPFR1(I,J)=**  " 

OATA  ( TI  TL E ( J)  ,J  = l,63)/3*  (“*  ’V*", , "A" , “L’VP", "N“#"A“ , "P«  , 

X-,E"t  WT'*|  ••  y 

OATA  (KEEPERK  33,NUMCOL)  , NUMC0L=2 ,95  ,2)  /3*  ( 2 ,4 , 6 ,6 , 1, 1, 1 , 1,  1,  2,  2 
72,2,2,  3,3)/ 

OATA  (KEEPERl(34,NUMC0L),NUMC0L  = 2,9t,2)/3*  C*  ’V  *V  “»H  ",n,2,4 
76,8,0,2,4,6,8,0,2)/ 

COL  SIZE*  NUMLTRS*  21 

WRITE  100, (TITLE (NUMCOL) , NUMCOL=l ,C0LSI 7E ) 

COLSIZ E*NUMLTRS*  2*63 

WRITE  99 , (TITLE(I) ,1=64 ,COLSIZE) 

WRITE  99 

COLSI7E*NUMLTRS*  33 
On  180  T a 1 « 3 ? 


V 


WRITE  101, (KEEPER  ( I , J) , J = 1 , COLS  I 7E) 

IF  (M00(I,2) .EQ. 0)  GO  TO  175 
WRITE  107, (KEEPER1 (I , J) ,J=1»C0LSIZE) 

GO  TO  19  0 

WRITE  102, (K EE PERI ( I , J) , J=l, COLSI ZE) 

CONTINUE 

C0LSIZE=32*NUMLTRS 

WRITE  103, (KEEPER1  (33  ,KUMCOL) ,NUMCOL= 2,COLSIZE, 2) 

WRITE  10E,  (KEEPER1  (3'*  , KUhCOL)  , NUMC0L=2,  CCLSIZE,  2) 

WRITE  09 
NUMLTRS=  0 
NUMPAG  E=  NUMPAGE* 1 

IF  (ALPHSKP.ANC.  (OLOALPH.NE. ALPHNUM)  ) 100,195 

NUMPAG  E=  0 

OLDALpH=  ALPHNUM 

WRITE  10t 

GO  TO  10  5 

IF  ( (NUWPAGE.EQ.2) .OR.DUMPIT)  NUMPAGE=0 

RETURN 

END 

SUBROUTINE  XFORM I T ( LTRSI ZE , WIDTH , HEI GHT , A n JUST , WRATIO , HRATI 0, 
XNSpACE .NUHROWS, NUMCOLS, A LPHNUM, LTRNUM, VEXF ANO,X FORM, IMAGE, TRIERS EJ 
INTEGER  L”RSIZE,  WTOTH,  HE  I GHT  , NSP AC  E,  -JU  -(ROWS , NJ  MC  OLE,  t lpHNUm  , . TP  NU -< 
LOGICAL  ADJUST  f REAL  WRAT10»HRATT0 
LOGICAL  VFXP  AND,  XFORM,  I ►'ACE,  IN  VERSE 
INTEGER  L ETTER  (3  2,  3 2)  ,NL'MLTpS,RECN'l>i 
roHMON/LFTTEPS/LFTTER,NUMLTPS,RECN')J 
INTEGER  RECFLTR 
REAL  FLIPLTR (ISO) , PP0T0<16G) 

COMMON /PROF  LTR/FLT  P LTP, E ROT U ,RECFLT^ 

COMMON/X FORKS/ FO  TLTR 
COMPLEX  FORTLTR( ft*  ,C4) 

COMPLEX  WORKFR  (f  ) 

INTEGER  POWERS  (7  ) 

INTEGER  NN (2) , NUMCOLS, NUMF OWS ,NC DL , N ROW, FLAG, TOP, BOTTOM, ITEMP, 

V.  COL, ROW,  LEFT, RIGHT  , IFOR  R,  I SIGN 

FORMAT (1HT) 

DIMENSION  TRPLTR (32,32) 

DATA  POWERS/2,4,  6, 16,32,64,128/ 

WRITE  2 

FLAG=0 

CONTINUE 

If  (ADJUST)  GO  TO  100 
NUMROW  S=  NUMCOL  S= LTRSI ZE 
WI0TH=32 
HEIGHT  =22 

CALL  TRPMLT (LETTER, TRPLTR) 

CALL  REPLACE  (TRPLTR , FOR  TLTR  , NUMROW5  , NUMCOLS ) 

GO  TO  235 
CONTINUE 

CALL  SIZE  IT (LTRSI7E, WIDTH, HEIGHT) 

NUMROW S= HEI GHT*HRAT 10 
NUMCOL S= WIDTH* WRATIO 
NRINGS=NSPACE**.5 

TF (NUMPOWS.LT. NR INGS)  NUMROWS=NRINGS 

IF( NUMCOLS .LT.NRINGS)  NUMCOL S=NR. NGS 

IF  ( MO D( NUMROWS, 2)  .EQ.l)  NUMROWS  = NUMROWS +1 

IF  ( MO 0( NUMCOLS, 2)  .EO.l)  NUMCOLS=NUMCOLS  + l 

TOP  =IA9S (NUMR0WS-LTRSI7E)/2 

IF  (NUMROWS.LT  .LTRSI7E)  GO  TO  150 

OOTTOM=T OP ♦l+LTRS I 7E 

ITEMP=TDP 

00  140  NCOL=l,LTPSIZE 
TOP=ITEMP 


00  110  NR0W=1,T0P 
pn»Ti  Ter  wpnu  ..mrm  l sD 


i i n 
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T 120 


DO  120  NROW*l,LTPSI7E 
TOP»TOP*l 

FORTLTR(TOP.NCOL)  =LFTTFR (NROW, NCOL) 

00  130  NR0W»80TT  ON,  NUMROWS 
FORTLTR (NROW  ,NCOL ) »0 
CONTINUE 
GO  TO  170 
CONTINUE 
TOP*TOP*i 

^OT  T OM=T  OP+NUMROWS-1 
00  160  NCOL*l, LTRSI7F 
ROW  = 0 

no  160  NR OW*TO°, BOTTOM 
ROWsROWM 

FORTLT M ROW, NCOL >=LETTFF  (NPOW.NCOL) 

CONTINUF 

LEFT=IA9*'(NUMCOLS-LTRSI',E)/? 

IF(  NUHCOLR.LT  .LTRSI7E)  C-0  TO  22C 

fight=left>ltrsi^e*i 

no  210  NROW=l , NUMFOWS 

no  180  NCOL=RIGHT,NUMCOLS 

FORTLTR (NROW, NCOL) =0 

ITEMP=RIGHT 

no  190  NCOLsl»LTPSI7E 

ITEMP* IT  FMP- 1 

COL=LTRSI7E-NCOL  *1 

FORTLTR ( NPOW,ITEMP) =FORT LTR < NROW, C 0L ) 

00  200  NCOL  * 1 » LFFT 
FORTLT  R( NROW , NCOL ) =0 
CONTINUE 
GO  TO  23? 

CONTINUF 

LEFT »LFFT  >1 

RIGHT=  LF  FT ♦NUMCOL  S -1 

DO  230  NROW  = 1 , NUMROWS 

COL*0 

00  230  NCOL=LFFT, RIGHT 
COL*COL*i 

FORTLT R( NROW, COL) = FORTLTR (NROW, NCOL) 

CONTINUF 

CALL  A RNGFPX  (FORTL TR, 64 , 6*.  , FORTLTR,  NUMROWS , NUMCOLS) 

CONTINUE 

IF  (FLAG.GT.O)  GO  TO  24f 

IF  (VEXPAND)  CALL  PRNTEX P (FORTLTR, NUMROWS , NUMCOL S, ALPHNUM ,L T * NUM, 
XMRATIO.HRATIO, WI OT H , HEIGHT ) 

IF  (XFORM. OR. IMAGE) 

X CALL  IVRTCPX ( FORTLTR » NUMROWS, NUHCOL  S) 

NN ( 1)  = N'JMROWS  2 NN  (2)  =NUMCOLS  fNniH  = ? $ ISIGN=-1  $ IFORM  = 0 
COL*RO W*  0 
00  240  ITEMP=1,7 

IF  (NUMROWS. FO.POWERS(ITEMP) ) R0W=1 

IF ( NUMCOL S.EO.POWFRS(ITEMP) ) COL*l 
CONTINUF 

IF  ( ROW.EQ.l.ANO.  COL.EQ.l)  GO  TO  250 

CALL  FOURT  (FORT LTR, NN,N01M, I SIGN, IFORM , WORKER) 

GO  TO  25  F 

CALL  FOURT  ( FORT LTR , NN, NOIM, ISIGN, IFORM, 0 . 0 > 


THIS  IS  THE  PL  AC  F TO  PUT  IN  AN  INVERSE  PRINT  ROUTINE  HAT 
FILTERS  THE  FOUKIFR  TRANSFORM  FROM  T HF  MIPOLE  AND  THFN  CALL 


jL  PPT 


CONTINUF 

IF  (FLAG.GT.O)  RETURN 

TF  (IMAGF ) CALL  PRNTIMG ( FORTLTR , NUMROWS, NUMCOL S , ALPHNUM , L TR NJ N) 
iF  (XFORM) 

imcxu  .unnriii  c. . 


■joooooodoooooooooooooooo  ooooonooooooonoooooooononnnnooo  non 
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X 1 , ALPMNUM  ( LTRNUM , p , ,T . , ,T.  ) 

FLAG»FLAGM 

IF  CXFORM.OF' .INAGR)  GO  TO  9i 

Kt  TURN 

END 

SUBROUTINE  FOURT  ( D A T A , NN  , NOI  M , I S I GN,  1 F IN  M , WORK) 


t *.  sj  m \j  y 


THE  COOLFY-TUKEY  FAST  FlUKITR  TRANSFORM  IN  USAS1  BASIC  FORTRVN 


TRANSF0RM(K1 ,K2, . . . ) = S l)M  <0  A T A ( J 1 , J 2,  . . . ' * E XP  ( I S I GN*  2*  P I S0?T(-1> 
*( (Jl-1)  * (Ki-1)/NN(1)  ♦( J2-i)*(K2-l) /NN(2> ♦ ...)>> , SUMMFD  FUR  ALL 
Jl,  K1  BETWEEN  1 ANO  NN ( 1 ) , J 2,  K?  BETWEEN  i ANO  NN<2),  L TC . 

THERE  IS  NO  LIMIT  TO  THF  NUMBER  OF  SUBSCRIPTS.  OATA  TS  A 
MULTIOIMFNSIONAL  COMPLEX  ARRAY  WHOSr  REAL  ANO  IMAGINARY 
PARTS  A5F  AOJACFNT  IN  STORAGE,  SUCH  AS  FC'TKAN  IV  PLACES  THEM. 

IF  ALL  IMAGINARY  FAPTS  ARF  7FRO  (PAt  A ARE  OISGUISro  REAL),  SET 
IFORM  TO  7FR0  TO  CUT  THF  RUNNING  TIME  BY  UP  TO  FORTY  PRPCEN’ . 
OTHERWISE,  IFORM  = »1.  THE  LENGTHS  OT  AH  DIMENSIONS  ARE 
STORED  IN  ARRAY  NN,  OF  LFNG1H  NPIM.  THEY  MAY  RF  ANY  POSITIVE 
INTEGERS,  THO  THf  PROGRAM  RUNS  FAST  F P ON  COMPOSITE  1NTFGFRS,  ANO 
ESPECIALLY  FAST  ON  NUMBERS  RICH  IN  F ACTOR*"  OF  TWO.  ISIGN  TE  M 
OR  -1.  IF  A -1  TRANSFORM  IS  FOLLOWEO  BY  A *1  ONE  (OR  A *1 
HY  A -1)  THE  ORIGINAL  DATA  REAPPEAR,  MULTIPLIED  By  NT OT  <=NN(J)' 

NN ( ?) * ...)  • TRANSFORM  VALUFS  ARE  ALWAYS  COMPLEX,  AND  ARE  RFTURNFP 
TN  ARRAY  OATA,  REPLACING  THE  INPUT.  IN  AOOITI ON,  IF  ALL 
DIMENSIONS  ARF  NOT  FOWt.FS  OF  TWO,  A'PAY  WORK  MUST  BE  SUPPIIFD, 
COMPLEX  OF  LENGTH  FOUAL  TO  THF  LARGEST  NOf  2**  K DIMENSION. 
OTHERWISE,  RFPLAf  WORK  BY  ZERO  IN  IMF  CALLING  SFQUFNCt  . 

NORMAL  FORTRAN  OATA  ORDERING  IS  EXPECTED,  FIRST  SUBSCRIPT  VARYING 
FASTEST.  ALL  SUBSCRIPTS  BEGIN  AT  ON F , 

RUNNING  TIME  IS  MUCH  SHORTER  THAN  THE  NAIVE  NTOT*- 2,  BEING 

GIVEN  BY  THE  FOLIOWING  FORMULA.  DECOMPOSE  NTOT  INTO 

?**K2  * 3»*K3  * ■♦♦KB  * ....  LET  SUMS  = 2*K2,  SUMP  * 3*K3  ♦ 

■4  ...  ANO  NF  = K 3 ♦ K5  ♦ ....  THE  1 IMF  TAKEN  BY  A MULTI- 
DIMENSIONAL TRANSFORM  ON  THESE  NTOT  DATA  ’S  T = 10  ♦ NTOTMU* 
T2*SUM2fT3*SUMF*T<i*NF)  . ON  THE  CPC  3300  (FLOATING  POINT  ADO  TIME 
OF  SIX  MICROSECONDS),  T = 3J  00  ♦ NTO T • ( 5 0 " 3* SUM2*I 8 * SUMFf 
320*NF>  MICROSECONDS  ON  COMPLEX  DATA.  TN  ADDITION,  THE 
ACCURACY  IS  GREATLY  IMPROVED,  AS  THE  RMS  ILLATIVE  ERROR  IS 

BOUNDED  BY  3*2**  ( -B) *SUM  < R ACTOR ( J ) * l.*»)  , WHER  L B NUMBER 

OF  BITS  TN  THF  FLOATING  POINT  FRACTION  AND  FACTOR(.I)  ARF  TMT 
PRIME  EATTORS  CF  NTOT. 

PROGRAM  PY  NCRMAN  BRENNTF  FROM  THF  "ASIC  PROGRAM  BY  CHARLES 
KADER.  RALPH  ALTER  SUGGESTED  THE  nr; A FO  THf  DIGIT  .EVERS®. . 

MIT  LI  NC  01  N LABORATORY*  AUGUST  in'.  tHIS  IS  The  pacTpsT  AN3  MOST 
VERSATILE  VERSION  OF  THf  FFT  KNOWN  :0  THF  AUTHOn.  SHORTER  050- 
GRAMS  rpUPi  ANT  FOUR 2 RESTRICT  01 MENSI TN  IFNGTHS  1 0 POWERS  Or  lWp. 
SEE—  TEFF  AUDIO  TRANSACTIONS  (JUNE  1BN7),  SPFCIAL  ISoJR  ON  nr, 


THE  OJPCRETE  EOUMFR  I Kfi  NSFOKM  PLACES  THK’  f RFST  RIOT  IONS  UPON  T MF 

DATA. 

1.  THR  NUMBER  OF  INPUT  DATA  AND  THt  NUMB!  R OF  1 RANSFDRM  VALJ'S 
MUST  BF  THE  SAMF. 

2.  BOTH  THE  INPUT  DATA  ANO  THF  TRANSFORM  VALUES  MUST  RFPRf ?r  NT 
(DUISPACrn  POINTS  IN  THfTF  RRSPFCTIVr  DOMAINS  OF  TINT  ANO 
FREQUENCY.  CALLING  THESE  SPAOINGS  OILTAT  AND  ITR.LTAF,  IT  MUSI  Bf 
TRUE  THAT  OElTAF="*rI/(NN(I)*OFlTAT)  , OF  COU5SR  , OILTAT  NL  f D NUT 
BE  THE  SAME  FOR  F VERY  DIMENSION. 

3.  CONCEPTUALLY  AT  LFAST,  THF  INPUT  OATA  ANO  THF  TRANSFORM  OUTBID 
REPRESENT  SINGLF  CYCLES  OF  PERIOTIC  FUNCT I ONS. 


LVAMOI  f « 


TMcrf.ntMUicinw.M  tiuwAtn  cnnui-B  Tunwcrn^M  nf  « 
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COMPLEX  ARRAY  DIMENSIONED  32  BY  25  ?Y  13  IN  FORTRAN  IV. 
DIMENSION  DATA  (3  2,26,  13)  , WORK (5  0)  ,NN<3) 

COMPLEX  DATA 
DATA  NN/32,25,13/ 

(70  1 1=1,32 
00  1 J =1 , 25 
00  1 K=l, 13 


1 OATACI ,J,K) =C0MPLEX  VALUE 

CALL  F0URT(0ATA,NN,3,-i,l,W0RK) 

EXAMPLE  2.  ONE-OIMFNSIONAL  FORWARD  TRANSFORM  OF  A REAL  ARRAY  OF 
LENGTH  64  IN  FORTRAN  II. 

OIMENS ION  DATA  (2,6**) 

DO  2 1=1,64 

DAT A(1»I)  = REAL  PART 

2 DAT  A (2 ,1 ) =0 • 

CALL  F0URT(DATA,6h, 1,-1, 0,0) 


DIMENSION  DATA  (1)  ,NN(1) ,IFACT(32) ,W0RK(1) 
TW0PI=6. 28310530? 

WR=0. 

WI  = 0. 

MSTPI=0. 

WSTPR=n. 

IF  (NDIM-1)  920, 1, 1 
NTOT=2 

DO  2 IDIM=1,NDIH 
IF(NN(I0IM) )023,°2Q,2 
NTOT=NTOT'  NN  < I PI  M) 


MAIN  LOOP  FOR  EACH  DIMENSION 


NP1  =2 

no  BIO  miM=  1 , NDIM 
N=NN  (I  DI  M) 
NP2=NP1*N 
TF(N-l)52Q,90a,6 


FACTOR  N 


M=N 

NTW0=NP1 
TF=1 
IDI V=2 

10  IOUOT=M/IDIV 
IREM=M-iniV*IOUOT 
IF(IQUOT-IOIV) 50,11,11 

11  IF  ( IRE  M)  20, 12,  20 

• 12  NT WO=NTWO  +NT  WO 

M*IQUOT 
GO  TO  10 
20  IDI V=3 

30  IOUOT=M/IOIV 
IRFM=M-IOIV*IOUOT 

IF  ( IQIJOT-IOI  V)  60, 31, 31 

31  IF(IREM)40,32,40 

32  IFACT ( IF) r I0IV 
1F»IF*1 
MxIQUOT 

GO  TO  30 

* 40  101  V=I  DI  V*2 

GO  TO  30 

50  IF  ( IRE  M)  60,51,60 

««  wTUn-MTUOiklTUn, 


✓ X 


60 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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GO  TO  7C 

IFACT(in  = M 

SEPARATE  FOUR  CASES  — 

1.  COMPLEX  TRANSFORM  OR  REAL  TRANSFORM  FOR  THE  4TH,  5TH,FTC. 
DIMENSIONS. 

2.  REAL  TRANSFORM  FOE  THE  2NO  OR  3RD  DIMENSION.  METHOD— 
TRANSFORM  HALF  THF  DATA,  SUP°LYING  THE  OTHER  HALF  BY  CDN- 
JJGATE  SYMMETRY. 

3.  REAL  TRANSFORM  EOF  THE  1ST  DIMENSION,  N ODD.  MET  HOO-- 
TRANSFORM  HALF  THE  DATA  AT  EACH  STAGE,  SUPPLYING  THE  OTHEP 
HALF  BY  CONJUGATE  SYMMETRY. 

4 . REAL  TRANSFORM  FOR  THE  1ST  DIMENSION,  N EVEN.  METHOD- 
TRANSFORM  A COMPLEX  ARRAY  OF  LENGTH  N/2  WHOSF  REAL  PARTS 
ARE  THE  EVEN  NUMBERED  REAL  VALUES  AMD  WHOSE  IMAGINAkY  PARTS 


C ARE  THE  OOP  NUMPEREO  REAL  VALUES.  rEPA9Al  E AND  SUPPLY 

C THE  SECOND  HALF  BY  CONJUGATE  SYMMETGY. 

C 

70  N0N2=NP1*  (NP2/NTWO) 

1C ASE=  1 

IF ( IOI M-M 71,90,P0 

71  IF(IFORM)  72,72,90 

72  ICASF=  2 

rF(IOI  M-l)  73,73, °fl 
7 3 IC  ASE=  3 

IF(NTWO-NPl) 90,90, 7M 
74  ICASE=4 

NTW0=NTW0/2 
N = N/2 
NP2  = NP  2/ 2 
NTOT=NTOT/2 
1=3 

DO  80  J=  2, NT OT 
OAT  A ( J ) = D AT  A (I) 

80  1=1+2 

90  I1RNG=  N°1 

IF (ICASE-2) 100, 9r  ,100 
95  I1RNG=  NPO  * ( 1 +NPR  FV/2) 

C 

C SHUFFLE  ON  THE  FACTORS  OF  TWO  IN  N.  AS  TWE  SHUFFLING 

C CAN  BE  DONE  BY  SIMPLE  INTERCHANGE,  NO  WORKING  ARRAY  IS  NEEOFD 

C 

100  IF(NTWO-NPl) 600,600,110 
110  NP2HF=NP2/2 
J*1 

PO  150  I2=l,NP2,NON2 


IF  ( J-I  2>  120,130,  1 30 
120  IlMAX=T2+NON2-2 

PO  125  11=12, UMAX, 2 
00  125  13=11 , NTOT , NP2 
J3=J+I 3-1? 

TEMPR=OATA (13) 
TEMPI=DATA  (13+1) 

OAT  A (I  3)  =OAT  A ( J3 ) 

OAT  A (I  3+ 1 ) =OATA( J3  + 1) 
DATA  ( J 3)  =TEMPR 
125  OATA(J3+l)*TEMPI 
130  M*NP2HF 
140  IF(J-M)150, 150,145 

145  J= J-M 

M*M/2 

IF  ( H-NON  2)  150,1mo,140 
150  J= J+M 


r . 

*.  c MAIN  LOOP  FOR  FACTORS  OF  TWO.  PERFORM  FOURIER  TRANSFORMS  OF 

* C LENGTH  POUR,  WITH  ONE  OF  LENGTH  TWO  TF  NErOEO.  THP  TWIDDLE  -ACTOt 

C N=EXP ( I3IGN*  2*  PI”  SORT  ( -1 ) *M/  (4*MM AX)  ) • CHECK  FOR  W*ISI  GN- SUl  ( -1) 

C ANO  REPEAT  FOR  W = ISIGN*SQRT ( -1)  'CONJUGATE (W) . 

* C 

N0N2T=N0N2*N0N2 

IPAR=NTWO/NPl 

310  IF(IPAR-2)  350, 330, 320 
320  IPAR=I  PAR/4 
GO  TO  31 P 

330  OO  340  11=1,11 RNG, 2 

DO  340  J3=I1 ,NCN?«NP1 
00  340  K1=J3,NT0T,N0N2T 


K2=K1*N0N2 
TEHPR=  OAT  A ( K 2) 

TFHPI= DATA (K 2+1) 

DATA  (K  2)  =DATA(K1)-TEMPR 
DATA (K2+1) =0ATA< Kl*l) -TEMPI 
DATA(Kl)  = DATA(K1)  +T  FMPR 
340  OATA(KIM)  =0ATA(K1+1) +TEMFI 
350  MMA  X=N  ON  2 

360  IF  ( MMA  X-NP2HF)  PG,60U,b00 

370  LMAX*MAX0 ( NON2T, MMAX/2 ) 

IF  (MMA  X-NON2)4  05,4  05, 38  G 
3 80  THETA=-TWOPI*FLOAT (N0N2) /FLOAT (4 ♦MM A X) 

* IF  (ISI GN) 4 00, 390 , 390 

390  THETA=  -THETA 

400  NR=COS  (THETA) 

«I=SIN  (THETA) 

NSTPR= -2.' WI*WI 
WSTPI=2.*WR*WI 

405  DO  570  L=N0N2,LMAX,N0N2T 
M=L 

IF (MMAX-N0N?)4  20  ,4 2 0,4 IP 
410  N2R=WR*WR-WI*WI 
W2I=2.-WR' WI 
W3R=W2R»  WR-W2I  *WI 

• M3I=W2R*  WI*W2I*WP 

420  00  530  11  = 1, II RNG, 2 

DO  530  J3=I1,NCN2,NP1 
KMIN=J3+IPAR*M 
IF  ( MMA  X-NON2)  4 30,4  30,4<*0 
430  KMIN=J3 
440  KOIF=IP«P/  MMAX 

450  *STEP=4*K0IF 

DO  520  K1=KMIN,NT OT,KSTEP 

K2=Ki*KDIF 

K3=K2*KDIF 

» X4=K3+KDIF 

IF ( MMA X-NON2) 460 ,460,480 
460  U1R=DATA (Kl) +DATA (K2) 

U1I«0ATA (Kl+1) «OATA(K2M ) 

U2R=DATA (K3) +OAT  A ( K4 ) 

U2I=0ATA (K3*l) ♦OATA(M  M> 

U3R=DATA (Kl) -DATA (K2) 

U3I =DATA (Kl +1) -0  A T A (K2+1 ) 
rF(ISIGN)470,475,475 
470  ULRsOATA  (K3  + 1)  -DATA(K*.*1) 

U4I  =OA TA  ( K4 ) -0  ATA  ( K3) 

GO  TO  510 

r 475  U4R*0ATA  (K4*l)  -DATA  (K3+1) 

U4I«DATA (K3) -OAT  A (K4) 

GO  TO  510 


• IW  l\  w 
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j 

490 

500 
510 

> 

5 20 

5 30 

540 


550 

560 

565 

570 


T2I*W2R*  DATA (K2+1) + W2I*  DA  TA ( K2) 

T3R*WR*0ATA  <K3)-WI*0ATA  <K3*1> 

T3I»WR*DATA  <K3U>  «-WI*0 AT  A <K3 ) 

T4R  = W3R‘ OATA  <K4> -W3I*DAT A <K4 +1) 

T4I  = M3R»  OATA  <K4*1> ♦H3I*P ATA ( K4> 

U1R=0ATA<K1)  ♦T2R 

U1I=DATMK1*1>  *T  21 

U2R=T3R+T‘‘R 

U2I=T3I*T4I 

U3R=0ATA (Kl) -T 2R 

U3I  = OA  TA  ( K1  ♦ 1)  -T2I 


lFCISIGN^gO^OO  ,K00 
U4R=T3I-T4I 
U4I=T4R-T3R 
GO  TO  510 
U4  R=T4 I-T  31 
U4I*T3R-T4R 
OATA(Kl)  * U 1 R *U  2K 
DAT  A (K 1*1)  = U II  HJ  ?T 
OAT  A (K  2)  ®U3R4U4R 
OATA  (K  2*1)  *U3I  HI4I 
OAT  A (K  3)  ®U1R-U2R 
DATA  ( K3*  1)  = U1I  -U2I 
DATACKO  =U3R-U4R 
OATA(K4+l) »U3I  -U4I 
KMIN®4  MKMIN-J3)  +J3 
KOI F*KST  FP 

IF (KDIF-NP2) 45 0,r  3 0,5  30 

CONTINUE 

M=MMAX -i 

IF ( ISI GN) 54  0 , 550  ,*60 
TEMPR=W* 

NRs-HI 
«I  = -TEMPR 
GO  TO  550 
TEMPR® WR 
«R»WI 
MI®TEMPR 

IF ( N-L  MAX) 565,565,410 
TEMPR=WR 

WR*WR* WSTPR-WI*WSTPI+WR 
MI  = HI*  WSTPR*TEMPF*WSTPI*W1 


IPAR®3  -I  PAR 
MMAXsMMAX+MMAX 

GO  TO  35  0 
C 

C MAIN  LOOP  FOR  FACTORS  NOT  EQUAL  TO  TWO.  APPLY  THE  TWIDDLE  FACTOR 

C W*EXP(  ISIGN*  2*  PI  SQRT  (-1 )'  (J2-1)  • ( J1-J2)  / (NP2*  IFP1) ) , THEN 

C PERFORM  A FOUR IFT  TRANSFORM  OF  LENGTH  IFACT(IF),  MAKING  USE  3F 

C CONJUGATE  SYMMETRIES. 

C 

600  IFCNTWO-NP2)t05, 7 00,700 

605  IFP l*NON  2 
I F=  1 

NP1HF*  NP1/2 

610  IFP2®IF°l/IFACTnF) 

J1RNG* NP2 

IF(ICASE-3)612,fcll,612 

611  J1RNG®  (NP2*IFP1> /2 
J2STP=NP2/IFACT(TF) 

J1RG2®  (J2STP  + IFP?)  /2 

612  J2MIN®  1*IFP2 


j 


* 


1 


y 

i 


H 


' 
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THETA* -TWOPI*FLOAT< J2-1)/FL0AT(NP2) 
IF<ISIGN>625  ,620,620 
1 620  THETA*  -THETA 

625  SI  NTH*  SI  N ( THET  A/  2 . > 

HSTPR*  -2.*  SINTH*  SINTH 
WSTPI=SI N (THET  A) 

WRsHST  PR  + 1 • 

HI*HST  PI 
JiMIN* J2+IFP1 

DO  635  Ji=JiMIN, J1RNG,IFP1 

UMAX*  J1+I1RNG-2 

DO  630  I1=J1,I1MAX,2 

DO  630  13=11, NTOT.NP2 

U3MAX* I3+IFP2-NP1 

00  630  J3=I3, J3MAX,NP1 

TEMPR=DATA<J3) 

OATA(J3)=DATA(  J3)*  WR-OATA  ( J3  + 1)  * HI 
6 30  OAT  A(J3+1)=TEMPR*WI+0ATA  (J3+D+WR 


6 35 
640 

645 

650 


* u 


655 


660 


665 


670 


675 


TEH PR* W* 

WR=WR* WSTPR-WI*WSTPI+WR 
WI*TEMPRBWSTPH-WP  WSTPR+WI 
THETA=-TWOPI/FLOAT (IFACT (IF) ) 
IF(ISIGN>650,645,645 
THETA* -THET  A 
SINTH*  SIN(THETA/2.) 

HSTPR*  -2.’  SINTH*  SINTH 
HSTPI*  SI  N ( THET  A) 

KSTEP=2*N/IFACT(IF) 

KRANG=KSTEP* (I  FACT ( IF ) / 2 ) +1 
00  698  I1=1,I1RNG,2 
00  698  13=11, NT0T,NP2 
00  690  <MTN=1, KRANG,KST£P 
J1MAX=I3  + J1RNG-IFP1 
00  680  J1=I3,J1MAX,IFF1 
•J3MAX*  J1+IFP2-NP1 
00  680  J3=J1, J3MAX,NP1 
■J2HAX=J3  + IFP1-IFP2 

K=KMIN+( J3-J1+(J1-I3)/IFACT (IF) )/NPlHF 
IF(KMIN-l)  655, 6*  r , 665 
SUMR=0  . 

SUHI=0. 

00  660  J2= J3, J2MAX,IFP2 
SUMR=SUMR+OATA  (J2) 

SUHI=SUMI +OATA (J2+1) 

WORK (K ) = SUMR 
WORK(KH)  = SUMI 
GO  TO  690 

KCONJ=K+2' (N-KMIN+i) 

J2= J2M AX 
SUMR*D AT  A ( J2) 

SUMI*D AT  A ( J2  +1 ) 

OLOSR*  0* 

OLDSI*0. 

J2=J2-IFP2 
TEMPR=SUMR 
TEMPI* SUMI 

SUHR*T  WO WR* SUMR-OL OSR  +DA1  A ( J2) 

$UMI*T WCWR  + SUMI-OLOSI +DA  T A ( J 2 + 1) 

OLOSR* TEMPR 
OLOSI* TEMPI 
J2*J2-IFP2 

IF ( J2- J3) 675 ,676 ,67  0 
TEMPR* WR*  SUMR-OLDSR+OATA (J2) 

TEMPI* WI+ SUMI 
uno*  nr^trMOfc.TCMPT 


f * ) kl  tl 


WORK(KCONJ)  = TEMPR*-TEMPI 
TEHPR=WR*SUMI-OLDSI+DATA (J2+1) 

TEMPI=WI*SUMR 
M0RK(K+1)=TEMPR+TEMPI 
W0RK(KCDNJ+1)=TEMPR-TEMFI 
660  CONTINUE 

IF(KMIN-1>685,685,686 
' 685  «R=WSTPR*1. 

NI=WST  PI 
GO  TO  690 
6 86  TEMPR=WR 

WR=WR* W5TPR-WI*WSTPI+WR 
WI-TEMPR*WSTPI +WIVWSTPR+WI 

690  TWOWR=WR*WR 
rF(ICASE-3>692,691,€92 

691  IF(IFP1-NP2)695,P92,692 

692  K=1 
I2MAX=I3+NP2-NPl 

00  693  12=13, I2MAX,NPl 
OAT  A (I  2)  =WORK  ( K) 

OATA(I2  + 1)=WORK(K*1> 

693  K=K+2 

GO  TO  698 
C 

C COMPLETE  A REAL  TRANSFORM  IN  THE  1ST  DIMENSION,  N ODD,  BY  30V- 

C JUGATE  SYMMETRIES  AT  EACH  STAGE. 

C 

695  U3MAX= I3+IFP2-NP1 

00  697  J3=I3,J3MAX,NP1 
U2MAX=  J3+NP2-J2STP 
DO  697  J2=J3, J2MAX, J2STP 
J1MAX= J2+J1RG2-IFP2 
U1CNJ=  J3+J2MAX+J2STP-J2 
00  697  J1=J2,J1MAX,IFP2 
K=1+J1-IT 
OATA(Jl)  =WORK(K) 

OAT  A ( J 1+1) =W0RK(K  + 1) 

IF(J1-J2)697,6  97,69C 

696  0ATA(J13NJ)=W0RK(K) 

DATACJ13NJ*1)=-W0RK(K+1) 

697  U1CNJ=  J1CNJ-IFP2 

698  CONTINUE 
I1F=IF+1 
IFP1=IFP2 

IF<IFP1-NP1>700,700,610 

C 

C COMPLETE  A REAL  TPANSFOPM  IN  THE  1ST  DIMENSION,  N EVEN,  BY  CDN- 

C JUGATE  SYMMETRIES. 

C 

700  GO  TO  (90C, 800,900, 701) ,ICASE 

701  NHALF=N 
N=N+N 

THETA=-T WOPI /FLOAT (N) 

IF(ISIGN)703,702,"’02 
* 7 02  THETA=-THETA 

703  SINTH=SIN(THETA/2.) 

WSTPR=  -2 •* SI  NTH*  SI  NTH 
NSTPI=SIN(THETA) 

«R=WSTPR*1. 

WI*HST  PI 
IMIN*3 

JMIN=2*NHALF-1 

GO  TO  725 

Tin  ll=  IMTW 


< A V 


V ‘ MAII 

DO  720  I*IMIN,NTnT,NP2 
SUMR*(0ATA (I) ♦OAT A (J) )/?. 
SUMI*(0ATA  (I ->1)  ♦DATA  (J*l  ) )/2. 

01 FR* ( OAT  A ( I ) -OAT  A ( J) ) /?, 

0IFI*(DATA(IM)-0ATA(JM)  )/2. 
TEMPR=WR*SUMI*HI*niFR 
TEMPI*  HI* SUM I -MR' 0IFR 
OATA(I)=SUMR*TEMPR 
OATAd  ♦1)  = 0IFI*TFMPI 
OATA (J )=SUMR-TEMPR 
OATA (J*l)= -OIF I+TFMPI 
7 20  ‘J*  J+NP  2 

IMIN*I MI  N*2 
JMlN=JMIN-2 
TEMPR* HR 

HR* HR* H5TPR-HI*HSTPI*WR 

mi*tempr*nstpi*hi*hstfr*ni 

725  IF(IMIN-JMIN)710  ,730,?4C 

730  rF(ISIGN)731,740,740 

731  OO  735  I = IMIN,NTOT,NP2 

735  DATA(I*1)=-0ATA(I*1) 

740  NP2=NP2*NP2 

NTOT*NTOT*NTOT 
J=NT0T *1 
IMAX*NTDT  72*1 
‘ 745  IMIN=IMAX-2*NHALP 

I* IMIN 
GO  TO  7?5 

750  OAT  A(J)=0ATA(I) 

OATA(J+l)=-DATA(I*l) 

755  r=I*2 

'J= J-2 

1F(I-IMAX)750,760,760 
760  OATA(J)=OATA (IMIN) -OATA (IMIN+1) 
OATACJ+l)  =0. 


< 765 
770 


775 


700 


C 

C 

C 

C 

800 
0 05 


810 
0 20 
8 30 


IF(I-J)77C,780,780 
DATA(J)=0ATA (I) 

OATA(J*l)=DATA  (1*1) 

1*1-2 
J*  J-2 

IF(I-IMIN)  7 75,77^,765 

OAT  A ( J ) = DAT  A (IMIN)  + 0ATA  ( 1 MI  N*l) 

T)ATA(JU)  = 0. 

IMAX=IMIN 
GO  TO  745 

OATA(1)=OATA (1)+DATA(2) 
PATA(2)*0. 

GO  TO  P00 


COMPLETE  A REAL  TRANSFORM  FOR  THE 
CONJUGATE  SYMMETRIES. 


2ND  OR  3RD  DIMENSION  PY 


IF  (I1RNS-NP1)  8 05,900,300 

00  060  I 3*  1 , NT  OT  , NP2 

I2MAX* I3*NP2-NPl 

TX)  860  12=13, 12MAX,NP1 

IMI N*I 2*11  RNG 

IMAX*I2*NPl-2 

JMAX*2*I3+NP1-IMTN 

IF(I2- 13)620, 820, 810 

JMAX=J MAX*NP2 

IF ( IOI M-2) 850,050,830 

J=  JMAX  ♦NPn 

nn  nun  T*TjjiJN.mAX.Z 


«UOUS<  O'  O'  O OOO  out) 


• W 


r'  9 


OATA(I ) = PAT  A (J) 

PAT  A <1  H)=-OATA(  J*l) 


640  JsJ-2 
8 50  JsJMAX 

00  660  I=IMIN, IMAX,NP0 
CATA(I)  = PATA(J> 

OATA(I*l»=-0ATA(  J*l> 

60  J= J-NP  0 

FNO  OF  .OOP  ON  EACH  DIMENSION 

00  NP0=NP1 
NPi=NP2 
10  NPREV=N 
20  PFTURN 
FNO 

SUBROUTINE  TRPMLT (LETTEF ,TRPLTR> 

INTEGER  L ETT FR  (3 2 * 32) 

OI MENS  ION  TRPLTR  ( 32 , 32) ,7RPWDW<32) 

1 FORMAT (1H1) 

2 FORMAT  (1H  , 24  ( **)  , /1H  , ***  TRAPEZOID  PARAMETERS  ♦“,/lH  ,2  C*‘ 

3 FORMAT (1H  LENGTH  OF  TRAPEZOID  PLATEAU  = ”,12) 

4 FORMAT  (1H  ,/,"  VARIANCE  IS  = ",F5.2> 

5 FORMAT  (1H  ,/,"  TRAPEZOID  WINDOW  ARRAY  TSt  **) 

6 FORMAT  (“  ",  1 0F6 • 3) 

7 FORMAT  (IH  ,/,"  THIS  IS  A NORMALIZED  WINDOW  ") 

IF  ( PS  KI P • E Q • **  OK  ” ) GO  T 0 45 

DATA  NORMALISED  WINDOW  SET  I PLEN  = 0 ANO  TP PHT=VARIANCE 
IPLEN  =0  $ TRPHT  =5.0 

THIS  PART  BUILDS  THE  TRAPEZOID 

INCLN= (32-IPLEN) /2 
IF  (IPLEN. NE.O)  GO  TO  9 
00  8 1=1,16 

8 TRPWOW (17-1) =TRPWOW(lfc*I ) =EXP(-( (1-1) **  2) / ( 2.* ( (TRPHT)  **2> ) ) 
WRITE  1 

WRITE  7 
GO  TO  4+ 

9 IF  (INCLN.EQ.O)  GO  TO  3E 
00  10  1=1,1 NCL  N 

TRPWOW (I )=TRPWDW( 33-1) = < 1/ ( I NCLN+1 . ) ) TRPHT 
10  CONTINUE 

00  20  1=1, IPLEN 
TRPWDW (I ♦INCLN)=1 • 

20  CONTINUE 

TF(  (32-(2*  INCLN)  ) .EQ. IPLEN)  GO  T3  43 

iodd=in:ln*i 

DO  30  I=l,IOOD 

TRPWOW (33-1) =(1/ (INCLN+2. ))  TRPHT 
30  CONTINUE 
GO  TO  43 

35  00  40  1=1,16 

40  TRPWOW  (D=TRPWDW(33-I)=1. 

43  CONTIN  UE 
WRITE  1 
WRITE  2 
WRITE  3,  IPLEN 

44  WRITE  h, TRPHT 
WRITE  5 

WRITE  6, TRPWDW 

THIS  PART  DOES  THE  MULTIPLICATION 


■*>  \J*J  U U X * X f o c. 

00  50  J=  1 » 3 2 

50  TRPLTR (J»I)=TRPWCW(I)* (FLOAT (LETTER (J,  t ) ) ) 
60  CONTINUE 
PSKIP= "OK" 

RETURN 

END 


U 


t 


PROGRAM  MAIN  (I  NPUT  , OUT  PUT  , T A PEA , T APF  ID 
C PROGRAM  MAIN(INPUT=R02B,0UTPUT=lJ02o,TAPE  = 3 Oi.  B , TA  •'El  -*  3u  19  ) 

1 FORMAT  (1HT) 

INTEGER  ALPHSTR,  ALPHSTP,  E XCPTNS ( 20) , NEXCP* S »TOTAL°H,LrKST?T, 

X LTRSTOP  » NUN SUB,  RANGE  (2, 1 0) , NSPACE ( 1 0 >, TAPEIN, TA PEOUT, ENORN , I P 
X »LSPAC E 

LOGICAL  RA NOIN, NFWFILE, ADD, CHANGE, Of  LETE 
WRITE  1 

ALPHSTR=i  S ALPHST  P = 15C  1 LTRST RT  =1  $ L' RST0P=2b  f TAPEIN-* 

RANDIN=.T.  $ NEWFILE  = .T.S  AOO=.F.  $ CM1NG‘=.F.  t 0<-LET£=.F. 

ENORM=  1 $ OC=0  f TAFECUT  =17  $ NUMSUP=1  $ NEXCPTS'OS  LSPASE=9  1 
XXXX=99  ? ENORM=l  5 NUMSUB=1S  LSPACE=615  NrXCPTS=6 
DATA  EX:PTNS/2*i,6?,8  2,i0i,106/,XXXX/l.r 
OAT  A NSPAC  E/ 1 0 *6 1 / 

DATA  RANGE/1 ,32/ 

CALL  PREPFLRU) 

CALL  PRE PRTC(TAPFOUT) 

CALL  PRTO°LD(ALPHSTR, ALPHSTP, FXCPTf  S,NEXCPTS,TOT ALPH, 

X LTRSTRT , LTRST OP, NUMSUB, RANGE, NSPACr ,EMORf,OC, TAPE  IN, RAN OIN,. SPA C 
X ,T  A PEOUT ,NEWFILE, ADO, CHANGE, OELETE) 

CALL  PRTOTBL (LTRSTRT , LTF>  STOP , TOT  °RT:  ,.  T E .,  .T  .,.  F.  , MSP  ACE  , 

XEN0RM,0: ,T A PEOUT , .F.,0) 

STOP 

END 

SUBROUTINE  ENORMFL  ( FLT  KLTR , NORiLT k,  NSPACE. , DC) 

PEAL  FLTRLTR(NSPACE) , NOR MLTR ( NSPAC E) 

INTEGER  NSPACE  ,DC 

INTEGER  START, NVFCTCR  t REAL  SUHSORS 

START=2 

IF  (DC.EQ.l)  STA  RT  = 1 
SUMSDRS=  P • 0 

00  100  NVFCT OR =S TART, NSPACE 
100  SUMSQRS=  SUMSQRS+FLTRLTR ( N VECTOR) * * 2 

SUMSQRS=SUMSORS*"  .5 
IF  (SUMSQRS. EQ.O .0)  SUMfORS=1.0 
DO  110  NVFCTOR=START, NSPACE 
110  NORMLTR(  NVEC  TOR)  =FLTRLTR  ( NVECTOR)  /St'MSORS 
IF  (OC.EO.O)  NORHLTR(l)=l 
RETURN 
FNO 

SUBROUTINE  GETF1.TR  < ALPHNUM,  LT  RNUM,  l SP  ACF  , NSPACE  , T APENUM , RA  NO.  M) 

INTEGER  ALPHNUM,LTRNUH,L SPACE, NS PACE ,T  ApE*  UM 

LOGICAL  RANOOM 

INTEGER  RECFLTR 

REAL  FLTPLTR  (180) ,PROTO(180) 

COMMON/PROFLTR/FLTRLTK,FROTO,RECFLT'? 

INTEGER  RECRQD,RFCSKIP,KEY 
l DATA  RECFLTR/1 / 

NUMHRDS=LSPACE+2 
IF  (.NOT. RANDOM)  GO  TO  95 
KEY  = ( ALPHNUM-1 )• 2F+LTRNUM 

CALL  REA  DMS (TAPfNUM, FLTPLTR. NUMWROS,KET) 

GO  TO  ISO 
95  CONTINUE 

KECRQD=( ALPHNUM-1)*  26 ♦LTRNUM 
RECSKI P=  IA RS  (RECFOD-RfcCFLTR) 

IF  (RECROD-RECFLTR)  100,1*0,120 
100  DO  110  1 = 1,  R EC  SKIP 
110  BACKSPACE  TAPENUM 
, GO  TO  14  P 

120  DO  130  1 = 1 , RECSKTP 

130  REAO  (TAPFNUM) 

i i»n  pfao  (taptmhmi  r FJ  tr.i  tr  (t  f .tsi  .NiiMurnsi 

i...  — -----  ■■  --  - 


RECFLTP=RFCR00+1 

150  CALL  R ED UC E(FL TR LTR, LS PA CE,FLTRLTR,N SPACE) 

RETURN 

END 

SU9R0UTI NE  PPEPFLR ( FORTA PE) 

INTEGER  FORTAPE 
INTEGER  FTNDEX  (3901) 

CALL  OPENMS (FORT APE, F1NDEX, 3501,3) 

RETURN 

ENO 

SUBROUTINE  PREPRT0 ( PROTA PE) 

INTEGER  PRO TAPE 
INTEGER  PI  NO  EX  (2^8) 

CALL  OPENMS( PROT A PE , PIND EX,2 88 , 0) 

RETURN 

END 

SUBROUTINE  PRT09L0 (ALPHSTR, ALPHSTP, FXCPTNS  > NEXC  PT  S, T QT  A LPH , 

V.  LTRSTRT ,LTRST OP, NUMSUB, RANGE, NSPACE,EMOF^, DC, TAPE IN, RANDIN,. SPACE 
X ,T  A PEOUT ,NEWFILE, ADD, CHANGE, DELETE) 

INTEGER  ALPHSTR, ALPHSTP, EXCPTNS(2C)  , NEXCP' S ,TOTALPH, LTRSTR I , 

X LTRST OP,NUMSUB,RANGE(2,10)  .NSPACEdO)  , TAPEIN, TAPED JT , ENORH  , 3C 
y. , LSPAC  E 

LOGICAL  RA NOIN, NEWFILE, A OD, CHANGE, DELETE 

INTEGER  OLDLTR,LTRNUM,CNTALPH(27) , PROTNUM, SUBLTRS, KEY,PSPaCF, 
y NCOL, NROW , A LPHNUM , W IDTH , NWORD, N VECTOR,  NEXI S T S 

REAL  PROTOS (169,10) 

INTEGER  MTAPLE(27,5)  ,TAPLF(15,1G)  , ET  A^L  F ( 26  U , 14 ) 

COMMON/TABLES/  MT ABLE, TABLE, ETA9LE 

INTEGER  RECFLTR 

REAL  FLTRLTR (180) , PROTO (180) 

COMMON /PRO FLTR/FLTRLTR, PROTO, RECFLT? 

OATA  MTA BL E/135*  0/ 

DATA  LTRINDX^LOLTRjCNTALPH^T^T^-’'*  0/ 

OLDLTR=LTPSTRT 

IF  (NEWFILE)  GO  TO  90 

CALL  READMS (TA PEOUT, HTAPLE, 135,1) 

IF  (ADD. OP, CHANGE)  GO  TO  90 
IF  (.NOT  .DELETE)  GO  TO  9 0 
OO  80  LTP.NUM=LTRSTRT  ,L  TP  SI  OP 
00  8(1  NCOL  = l,  5 
80  MT  ABLE  (LTRNUM,  NCOL)  =0 

90  DO  230  LTRNUM=LTPSTRT,LTPSTOP 

IF  (CNTALPH(OLDLTR) .GE.MTA8L E( OLOLTR ,3 > ) GO  TO  10F 
C IF  THE  OLDLTR  EXISTS  SAVE  THE  OLOLTR  PROTOTYPE  INFORMATION 

CALL  WRITMS (TA PEOUT , TABL E ,NUMWRDS, KFY) 

DO  100  PROTNUH=l ,SUBLTRS 

100  IF  (TAB_E(i,PROTNUM) . EQ.l)  CALL  WRIThS ( T APEOUT , PROTOS ( 1 , PRO  T N'JM ) 

*.T  ABLE  (6,PROTNUM)  , TABLE  ( 2 , PROTNUM)  ) 

105  IF  (CNTALPH (LTRNUM) .EQ.C)  GO  TO  123 

C IF  THE  NEW  LETTER  ALREADY  EXISTS  BUT  IS  NOT  COMPLEX  RE  AO  IT 

C INTO  CORE  TO  BE  WORKED  ON 

KEY=MTA9LF (LTRNUM, 2) 

SU9LTRS  = MTA  BLE (LTRNUM, 3) 

NUMWR0S=SU8LTRS* 15 

CALL  R EA  OMS ( TA  PEOUT ,TABl E ,NUMWR03 ,KFY ) 

00  115  PR0TNUM=1,SUBL1RS 
TF  (TABLF(1, PROTNUM) .EQ.O)  GO  TO  iir 

CALL  READMS(TAPEOUT,PKOTOS(1, PROTNUM) ,TABl  E (8, PROTNUM) , 

7.  TABLE  (2, PROTNUM ) ) 

GO  TO  115 

110  DO  112  NROW=l,PSPACE 
112  PROTOS (NROW, PROTNUM) =0 
115  CONTINUE 
GO  TO  160 

4 


L M 


130 

140 

150 

160 


170 

175 


C# 

c# 


180 

19  0 

20  0 


205 


210 


0 


MTABLF  (ITRNUM,  1)  = 1 

MTABLE (ITRNUM, 2) =KFY=LTF  NUM  + 1 

MTABLE (LTRNUM, 3) = SUPLTRS=NUM SUR 

NUMHRUS*SUBl  TRS*  15 

MTABLE(LTRNUM,5) =T  OT  AL  PH 

00  150  NCOL=l»SURLTRS 

LTRINnX=LTRINOXM 

TABLE(1,NCOL>=0 

TABLE( 2.NCOL) =LTRINDX 

1ABLE(3, NCOL)= LTRNUM 

TABLED, NCOL)  = NCOL 

TABLE(E,NCOL)=RANGE(l,NCOL) 

TABLE(6,NCOL)* RANGE (2,NCOL) 

T ABLE( 7,  NCOL ) =RANG£ (2 , NC OL  > -RANGE ( 1, NC OL ) +1 
TABLE ( 8, NCOL)=PSPACE=NSFACE(NCOL) 

TA8LE(9,NrOL) = ENORM 
T ABLE  ( 10 , NCOL)  =DC 
DO  130  NROW  = 11,15 
TABLE ( NR OW,NCOL)=0 
DO  140  NROW*l,  PSPACF 
PROTOS  (NROH.NCOU  = 0.0 
CONTINUE 
CONTINUE 

DO  200  ALPHNUM=AIPHSTR,ALPHSTP 
IF  (NEXCPTS.FQ.O)  GO  TO  175 
DO  170  NROW  *1 ,N  EXCPTS 

TF  (ALPHNUM.EQ.FXCPTNS(NROW  ))  G3  TC  20! 

CNTALPH(LTRNUM)=ONIALFH(  LTRNUM)  +1 
CALL  GETFLTR(  ALRHNUM,LTRNUM,LSPACE,L  SPAPF  ,TAPEIN,D.6NDI  N) 
WinTH=FLTRLTR(LFPACE+l) 
no  190  PROTNUM*l , SURLTRS 

IF  ( (WIOTH.LT.TAPLE  (5  ,PE  OTNI'M)  ) .OR  . ( WIOT  H . GT  .1  A RLE  ( 6 , PROTNUM)  ) > 

X GO  TO  19C 

PSPACE  =t AHLF (8 .PROTNUM) 

TABLE(12,PROTNUM)=TABLE(12, PROTNUM) M 
MASKER=M ASK ( 1) 

NWORO=EO-MOO(ALPHNUM-l,f  P) 

MAS  KER=S  HI  FT (MASKER, NWORO) 

NMORD=  (ALPHNUM-1 ) /GO+13 

TABLEt  NWOKD, PROTNUM) = MAS KLR. OR. TABLF ( NWORO , PROT NUM) 

CALL  R EDIIC  F ( FL  TRl  TR , LSPACF , DROTO,  PSPACE) 

XX  XX=99  • $ IF(ENORM.FC.l)  CALL  ENOKlFL(r*ROT0,  PR  OI  O, PSPACE  , DC ) 
XXXX*9P«  CALL  LFARKLTk'KUM, PROTNUM, FNORM, PC, PSOTO.PSFACL) 

DO  180  NVECTOR=l, PSPACF 

PROTOS  (NVFCT  OR  ,PF  OT NUM)  «PPOT  0 (NVECTPP)  RPRPTOS(NVECTOr  , PROTNUM  > 

TABLE( 1,  PRO INUM) =1 

CONTINUF 

CONTINUE 

OLOLTR=LTPNUM 

IF(CNTALPH(LTRNUM)  .GE.TOTALPH)  GO  TO  21r< 

GO  TO  230 
NFXIST  S*  0 

DO  220  PRO T NUM  *1,SUBLTPS 

IF  (TAPLF  ( 1 , PROTNUM)  . E Q . L ) GO  TO  22r 

NFXISTS«NEXIS7SM 

PS  PACE =TABLE (8, PROTNUM) 

on  210  M VFCT OR =1 , PSPACE 

PROTOS (N VF C T OR, PROTNUM )=PKOT OS (NVEC I OP, PROTNUM) /T  ABL  E ( 1 2 , PR  Cl NUM) 
v*XX»99.  « IF(ENOFN.EO.l)  CALL  ENORMF  L ( PRO' OS  ( 1 , PROT NUM)  , 

< PROTPS( 1 .PROTNUM) , PSPACE , DC)  JXXXX=99. 

' A ml  PC  11  .PPOTNUM) =TABLF( 12, PROTNUM) 

*i**»4Mr  CALL  LE AP 2 (PRO! OS ( 1, PROTNUM) , PRO, NUM, PS  PACE  » DC , ENORM) 

».L  NR1  TMS  (TA  Pr  i'UT, PR  OT  OS(l,  PRO!  NUM), nSPACE, TABLE  (2,  PROTNUM)  ) 
i*n*iT  xnhf 

( TRNU**,  4) -NEXISTS 


IT  A O-.OllT  . T AIM  f.  . NPMW  Rns.KCV) 


2 30 


C? 

C? 


10 

11 

12 

1 

i 

25 

30 

35 

40 

41 
45 
50 


80 

85 


15 

20 


CONTINUE 

MT  ABLE (27,1) =LTPTNDX 
CALL  WRITMS (TAFFOUT.MTAPLF, 135, 1) 
XXXX=99$  CALL  PRT OS7  X (LTF STRI  ,LTRSTOP 
• F.  , »F ,, NS  PACE ,FNORM,CC,TAPFOUr , .F 

return 

END 

SUBROUTI NF  PPTOTPL (LI RST  KT .LTRSTO0, 
SKIPAGF, 

OVERIOF,OSPACE,OFNORM,OrC ,T  A PENUM , Nf 
INTEGER  LTRSTRT, LTRST UP , T OTPRT 0, T APF 
,oenorm, one 

LOGICAL  PF INTI T, ^K I PAGE , OVER  IDE  , NT 
INTEGER  L ETT  EKS ( ) , MASKER, LTRNUM, Kf 

nspace,pspacf,nccl,nrow,nalphs,nte^ 

ALPH(  15  0)  ,ENORM,OC,NINrFX(  300) 
INTEGER  MI  ABLE (2  7 ,5) ,TAPLF (15,10) , ET 
COMMON n A PL FS/  MT A 8L E ,T A PL E , ETARL E 
INTEGER  RECFLTR 
REAL  FLrRLTRC18O),PROTO(180) 

COMMON /PPOFLTR/FLTKLTKjFROTU.PECFLTF 
FORMAT ( 1 H 0 ) 

FORMAT  ( 1H1 ) 

FORMAT  ( 1 H ,41X,54 ("*") ,/,lH  ,41X,"C-< 
F ”12" , LETTERS  "12”  70  "12,/tlH  ,41 
FORMAT ( 1HC  » 4 2 ( ”) / 1H  , ” L T R SUBLTR  R 
“/  1H  ,42  ("*’*)  ) 

FORMAT  (1H  ."CLASS  #”12",  LETTER  "A1 
FORMAT  (=("*”)) 

FORMAT  (1H  , "CLASS  *"12",  LETTER  ”A1 
FORMAT (1HP, "TOTAL  # OF  F LPHABFTS="I? 

* OF  (POSSlPlt)  SUBCLASSF 
"SUBCLASSES  IN  ASCENDING  ORUER  IS«"I 
FORMAT (1H  ,”#  OF  (ACTUAL)  PROTOTYPE? 
"THESE  PROTOTYPES  ARE  SHOWN  RELOW. "/ 
FORMAT (1H  , ”LT R SUBLTR  FANGt  SPACE  F 
IN  THE  PROTOTYPE") 

FORMAT (1H  , IX, A1 , 3X  , 1 2 , 3 X , 12  , , I 2 , 

22  (13  ,",”)  /,f>  (1H  ,44X, 22(13,",")  /) ) 
FORMAT (1H  , IX, Al« 3X, 12, ?X,  12,  "-M, 12, 
FORMAT  (1H  ,42("-")) 

FORMAT (1H  ,1X,A1,3X,I2,3X,I2,"-",I2, 
"THIS  RURLTR  OOFS  NOT  EXIST,  NO  ALP 
PAT  A L ETTERS/"A”,  "B”,  **C”  , "0"*  "E",  **F" 
"N", "0  ”, "P","0", "R", "S", "T", "U", “V", 
T OT  PRT  0=  0 

IF  (OV  E?  IDE  . AND. NEWT APF ) CALL  OPENMS 
CALL  READMS(TAPENUM,MIAPLE, 133,1) 
MASKER=M ASK ( 1) 

IF  (SKIP AGF)  WRITE  5 

(PR INTIT .AND. SUMMARY)  PRiNTIT  = ,F. 
(PRINTIT . OR  .SUMMARY ) WRIT  F 10,  (Tr 
(SUMMARY)  WRITE  11 
160  LTRNUM=LTPSTRT,L7RSTOP 
( .NOT . PR I NT  IT ) GO  10  85 
WRITE  1 
WRITE  15,(20) 

IF  (MTABLE(lTRNUM.l).FQ.l)  GO  TO 
WRITE  12,  (LTRNUM, LETTERS  (LTRNUM)  ) 
15,(20) 

16  0 

20  , (LTRNUM.L  ETT  LRS  (LTRNUM)  ) 

15, (20) 

IF  (MT ABLE (LTRNUM, 1) ,r0.0)  GO  TO 
KEY*MTA3LF (LTRNUM, 2) 


,LSPACF,TAPEIN,PAND1N,.F.  , 
«f.F.,.F .,1,1) 


T0TPf?'0,PRINTn,  SUMMARY , 

WT  APE  ,T  APFTMP) 
NUM,TrPETMP,OS°ACF(lC) 

APE, SUMMARY 
Y , 3URI  TRS, PROTNUM, 

P,  ALPHNUM, ALPHSTP, U PHST’ , 

ABLE (26 0,14) 


ARACTTRIsriCS  Or  PRGTOTY’F.  FJi 

X,  34  ("»••)) 

ANGE  SPACE  ENORM  DC  TERM  AALPH 
•*  *♦  DOES  NOT  FxiST  * »'"> 

) 

S = ‘*12  , / RANGE  PF(POSSIB.F)  " 

2,  "-"T  2 » 9(  ”,  ”,  12,"-"  12)  ) 

="I 2” • CHARACTERISTICS  " 

) 

NO IM  nc  TERM  AALFH  ALPHABETS 

3X,I3,3X,A3,4X,A3,4X,I?,3X, 

3X,I3,3X,#3,4X,A3,4X,I3,3X) 

3X,13,3X,A3,4X,A3,4X,IT,JX, 

HA  BETS  IN  THF  RANGE") 

* "G ",  *’H  " , "I  ** , **  I”,  "<  ",  "L  ”,  "M", 
V* 


IF 

If 

IT 

00 

IF 


WPITE 
GO  TO 
WRITE 
WRITE 


*W”,  "X",  "Y  ”,"Z“/ 


(TAPS'  MP.NINDEX,?  A<3,  0) 


PFNUM , LTRSTRT, LTRST OP) 


81 


lhl 


, NUMWRDS5  SUBLTRS* 1 5 

CALL  RFA0MS(TAPENUM,TA8LE,  NUMWKDS , KF  Y ) 
M * IF  ( .NOT • PRINTIT ) GO  TO  90 


WRITE  25  , ( MT  A8LF (LTRNUM,5 ) , 

X SUBLTFS,  (TABLE(  5 , NCOL ) , T6  RL  E (F , NCOL  > , NCOL  = 1,  SURL  dS ) ) 

WRITE  30,  (MTA  RL  F ( LTRNUN  ,<> ) ) 

WRITE  15,(71) 

WRITE  35 


WRITE  15, (71) 

90  CONTINUE 

DO  ISO  PROTNUM=i,SUBLTRS 
IF  (TA9LE(1, PROTNUM)  .FO.C)  GO  TO  1!-? 

TOTPRT  0=T0TPRT0  +1 
IF  (.NOT.  OVERIOF)  GO  TO  100 
PSPACE=T APL  E (8  , PROTNUM) 

TARLE( 8, PPOTNUM) =NSPACE=OSPACE(PROTNUM) 

T ABLE ( 9,  PROTNUM)  =OENOKM 
TABLE (10, PROTNUM)=OOC 
IF-  ( .NOT  .NEWTA  PE)  GO  TO  100 
KEY=TA«LE(  2,  PROTNUM) 

CALL  REA  DMS  ( TA  PENUM, PKOT O, PS PACE, KEY) 

CALL  REDUCE  (PROTO, PSPACE, PROTO, NSPACE) 

IF  (OENORM. FO. 1)  CALL  ENCRMFL ( PROTO, PROTO , NSPA C E , ODC) 

CALL  WRITMS ( T A PETMP, PROT O , NS PACF, KEY ) 

100  CONTINUE 

DO  110  NCOL=2, IS 

110  CTABLE (TOTPRTO  , NCOL- 1) =TA8L E(NCOL , PROTNU*) 

IF  (SUMMARY)  GO  TO  142 
IF  (.NO!. PRINTIT)  GO  TO  150 
NALPHS=0 

OO  140  NWR0S=1,3 
NTEMP=TABLE( 12*NWROS,PRCTNUM) 

ALPHSTR=  (NWROS-1) ♦60+1 
ALPHSTP=ALPHSTR*r9 
IF  (ALPHSTP.GT.1F0)  ALPHSTP=150 
00  130  ALPHNUM=  ALPHSTR, ALPHSTP 
IF  ( NT  EM  P.  A NO.  MASKER)  125,127 
125  NALPHS=NALPHS«-1 

ALPH(NALPHS) =ALPHNUM 
127  NTEMP=SHI FT  ( NT  EMP,  1 ) 

130  CONTINUE 
140  CONTINUE 

142  NCOL=TABLE(12, PROTNUM) 

ENORM=  DC  = "OUT‘* 

IF(TABLE(9,  PROTNUM)  .EO.l)  ENORM=**IN  ” 

IF(TABLE(10, PROTNUM) .EO.l)  DC="IN  " 

IF  (TABLF(1, PROTNUM) .EQ. 0)  145,147 

145  WRITE  50, (L  ETT  ERS (LTRNUN ) ,PKOTNUM,T A PL E (F , PROTNUM) , TAIL F (C, PROTNUM 
X) ,TABLF(8, PROTNUM) , ENORM , DC , NCOL) 

IF  (SUMMARY)  150,149 

147  IF  (.NOT. SUMMARY)  GO  TO  148 

WRITE  41,  (LFTTERS(LTRNUF)  , PROTNUM,  TABLED,  PROTNUM)  , TABLE  (b,  PROTNU* 
X) ,TABLE(8, PROTNUM) ,ENORM,OC, NCOL) 

GO  TO  150 

148  WRITE  40, (LETTERS ( LTRNUM ) ,PkOTNUM , T A BL E (5 , PROT NUM) , T A BL  E ( b, PROTNUM 
XI ,TABLF( 8, PROTNUM) ,FNORM,OC, NCOL  , ( ALpH ( NROW) , NRO W=1 , NA L PHS) ) 

149  WRITE  45 
15  0 CONTINUE 
lfcO  CONTINUE 

RETURN 

END 

SUBROUTINE  REDUCE  ( INFLTR , I NSPACE, OUTFL TR , OTSPAC E) 

INTEGER  INSPACE, OTSPACF 


REAL 


INFLTR (INSPACE) ,OUTFLTR(OTS°ACF) 

i>inirMaMgug|glgADQHN.QUT\/CTR|INVCTRtSTART,ST  OP 


00  85  NVECT0R=1,PTSPACE 
85  OUT FLTR(NVEC TOR) = I NFLTRiN VECTOR) 

GC  TO  125 
90  CONTINUE 

INHRMNC=  INSPACE,J,.f  J OUTHRMC=OTSPACE  *.5 
SKIP=( INHRMNC-OUTHRHC) *c 
D0HN=0UTHRMC/2*1 
OUTFLTR(l) =INFLTR(1) 

TF  (OTSPACE.EQ.l)  RETURN 
00  100  OUTVCTR=2,OUTHRMC 
100  OUTFLTP(OUTVCTR) =INFLTR (OUTVCTR) 

INVCTR=START=OUTPRMC*i 
DO  120  NR0W=2,  DOWN 
INVCTR=INVCTR+SKIP 
STOP=START ♦OUTHRMC* 2-1 
HO  110  OUTVCTR=START,STCP 
OUTFLTP( OUTVCTR) =1 NFLTR ( 1 NVCTR) 

1NVCTR=INVCTR+1 
110  CONTINUE 

START=ST0P*1 
120  CONTINUE 
125  CONTINUE 

00  130  1=1,4 

13  0 OUTFLT R( OTSP AC E + I ) =INFLTR ( I NSPACE+I ) 

RETURN 

END 

SUBROUTINE  PRTOS TX ( LTRSTRT , LTHSTOP, L S PACE , LTRT APE, RA NDOM , PR  I N I . , 

*/.  SKIPAGE,OVERIOE»OSPACE,OENORM,OOC»PROTAor,COPOLTN,HSTGGRM» PA  •?«?, 

X NDEV,NOIVS,WRAT10,HRATIO) 

INTEGER  LTRSTRT,  LTRSTOF , OSP ACE < 10 ) , PROT APE  ,LS PACE , NDEV , NOT  VS , 
•/.LTRTAPE,  OEN ORM  ,0 DC 
REAL  WPA T 10 , HR ATT 0 

LOGICAL  PRINTIT, SKI PAGE, OVER IDE , COROLT  N , HSTOGRM , BA  PS, RANDOM 
INTEGER  MT ABLE (2/ ,5) , T ABLE ( 15 , 10 ) ,ET AOLE( 260,14) 

INTEGER  RECFLTR 

REAL  FLTRLTR  (160) ,PROTO(ieO) 

COMMON/PROFL  TR /FLTRLTR, FRO TO, RECFLTR 
COMMON /TABLES/  MT ABLE, TABLE, ETABLE 

INTEGER  ALPHSTR,  AL  PHSTP,  NWRDS,PROTNUM , NT  E>'P,  TCTPRTO,  RLPHNUM 
X , PE NORM, PPC»T0TALPH,PSPACF»SU8LTR 
MASKER=M ASK ( 1) 

CALL  PRTOTBL(LTRSTRT,LTFSTOP,TOTPRTO,PRTKTIT,.F.,SK1PAGE,0 VERIOE, 
•/.OSP  ACE,  DENORM,  OOC,  PROT  APE,.  F.,0) 

00  120  PPOTNUM=l  , TOTPRTO 
LTRNUM  = ETABLE  (PROTNUM,  2) 

TF  ( (LTRNIIM.LT. LTRSTRT  ).  OR.  (LTRNUM.  GT  .LTPSTOP)  ) GO  TO  120 
SUBLTR=  ETABLE  (PPOTNUM, 3) 

TOT ALPH=  ET  A RLE (PROTNUM,ll) 

PSP ACE=  ETABLE  (PPOTNUM,?) 

PENORM=ETABLE( PR0TNUM,8) 

POC=ETA3LE(PPOTNUM,B) 

CALL  REA OMS ( PROT AF E, PROT  C,  PS  PACE, ETA  RL  E( PPOTNUM, 1) ) 

IF(IO. NE. 23) PRINT' ,"  PROT  “.PROTO  S TO=23 
DO  110  NWkOS=l,3 
ALPHSTR5  (NHROS-1) 

ALPHSTP  = ALPHST  R +59 

NT EMP= ETABLE (PROTNUM, ii+NWROS) 

DO  100  ALPHNUM=  ALPHSTR, ALPHSTP 
IF(NTEMP. AND. MASKER)  95,97 

C 5 CALL  STATS(ALPHNUM,ALFHNUM,T0TALPH,LTRNUM,SUBLTR,LSPt3E,PSPA:E, 

CCCCCXCLTRT APE,RANOPM,PENORM,FOC,COROLTN,HSTOGRM,BARS,NncV,NOIVS,  TCCCCCr 
CCCCC7.CPR0TA  PE,  WPATIO,  HP  AT  10)  CrCCCCCCCCCCrCCCCCl'CCCDCCCCOCrCCCCCCDCCCrCO 
95  CALL  GETFLTR  (ALPHNUH, LTRNUM, LSPACE^SPACE, LTRTAPE, RANDOM) 

CALL  STATX( FIT RLTR,PSFAf  E, TOTAL  PH) 

97  NTEMP=SHIFT(NTEMP,1) 

inn  rnMTTMiiF  


1 


110  CONTINUE 

CALL  M AH  PA  K ( PRGT  0 , PSP ACE ) 

CALL  WRIT NS  ( PROT A PE , PROT 0, PSPACF , ETA BL  E(PrOTNUM, 1)  ) 

120  CONTINUF 
RETURN 
END 

SUBROUTINE  ST ATX ( F VFC , NENT , NTOT) 

C**********  INSFRT  TO  FROTO  BUILD  2/3/7F  JRL 

COMMON /STX/  STND0EV(81,3) ,FMEAN(31,3) 

DIMENSION  FVFC  (NFNT) ,FT(2,81> 

IF  (INUM.FQ.O)  GOTO  60 
DO  50  J=  1 , 169 
50  FT (1, J)=FT (2, J)  = 0.0 
60  CONTINUE 

INUM=INUM*1 
DO  10  0 J=  1 » NENT 

CALL  DOIT (FVFC  (J) ,STNDOEV(J) , FMEAN ( J > , FT ( 1 , J ) , I NUM) 

10  0 CONTINUE 

IFtlNUM.  EO.NTOT)  I NUM  = 0 

RETURN 

FNO 

SUBROUTINE  DOIT (FIN, FSD, FM,FT»INUM) 
c*»***«»**«  INSFRT  TO  PROTO  BUILD  ?/3/Tf  JRL 

DIMENSION  FT ( 1 ) 

FT  ( 1)  * FT  ( 1 ) ♦ FIN**2 
FT ( 2) = FT ( 2 ) ♦ FIN 
FM  = FT  ( 2)  /INUM 

FSD  = (-T(1)/INUM)  - ( (FT (2) /INUM) ♦* 2) 

FSD  = SORT (A  BS (FSD) ) 

RETURN 

END 

SUBROUTINE  MAHPAK ( PRO! 0, PSPACE) 

INSFRT  TO  PROTO  BUILD  2/8/7E  JRL 
COMMON/STX/  STNDDFV (81,3) ,FMEAN(31,3) 

INTEGER  PSPACE 
DIMENSION  PROTO(RSPACE) 

IF ( IQ . NE • 23)  PRINTS,”  FN  ", FMEAN 

IF ( IQ. NF • 23 ) PRINT*,”  SDEV  ”,STNOOEV  * IP=23 

DO  20  1=1, PSPACE 

20  PROTO ( I) =PACK( FMFAN(I ) ,STNDOEV(I>  ) 

RETURN 

FNO 

REAL  FUNC TION  PACK(A,B) 

C**********  INSERT  TO  FROTO  BUILD  2/8 /?f  JKL 

IA=A*1E5  t IB  = B* 1E5 

PA CK*(  SHIFT  ( IA  ,3  0 . AND. M ASK ( 30) ) .OR.  . NOT  . N ASM  30  > . AND . I P 
IF(IA.LT.O)  PACK=PACK.OF .MASK(l) 

IF(IB.LT.O)  PACK  = PACK. OP. SHIFT ( MAS K ( 1 ),' 0 ) 

RETURN 

END 

SUBROUTINE  LEAR1 (N EXT , I SET , NORM, I DC,  FVFC , NOIM) 

C**»”**»”  INSERT  TO  PRTOBLD  2/17/78  JRL 
COMMON /STX 1 /FT (2 , 8 1 , 3 ) , J NUM ( 3) 

DIMENSION  FVFC(NDIM) 

COMMON /STX/SDEV(61 ,3) ,FN (81, 3) 

DATA  IGO/O/ 

IF(NEXT.NE.IGO)GOTO  20 
10  CONTINUE 

C *** * • • **  UPDATF  ONE  OF  3SFFCIFIE0  HEANANB  SPEV  VECTORS 
INUMdSET  ) = INUM(TSET)  *1  f I X=INUM ( I SET ) 

TF ( NORM. EQ, 1 ) CALL  FNORMFL (FVEC.FVED , NOIM, I DC) 

DO  15  J=  1 , NOI M 

CALL  DOIT(FVfC(J)  , SOE  V ( J,  ISFT) , F H ( J,  I SET)  , FT  (1 , J , I SET ) , I X) 
15  CONTINUE 
RETURN 

?n  rnwTTMiiF 


REINITIALIZE  FOR  EACH  NEW  SET  OF  3 VECTOR  STREAMS 


IGO=NE  XT 
00  30  J=  1 » 8 1 
00  30  K=l,3 
INUM(K ) = 0 

30  FT(l,J,<)sFT (2,J,K>=0.0 
GOTO  10 
END 

SUBROUTINE  LEAR2 (PROTO, ISET, PSPACE, OC, NORM) 

**********  r NSERT  TO  PRTOBLO  2/17/78  JRL 
C0MM0N/STX1/FT  (2 , 8 1 , 3) , I NUM ( 3) 

COMMON /ST X/SD( 81 , 3) ,FM(8i,3) 

INTEGER  PSPACE, START 
DIMENSION  PROTO(l) 

START=2f  IF(OC.EO.l)  START  = 1 

IF ( NOR  EO  . 1 ) C AL  L ENORMFL <FH ( 1 , 1 SET)  ,FM( 1,ISET) , PSPACE,  DC) 
00  10  I=START, PSPACE 

S0<I,ISET)=FT(i,I,ISET)/INUM(ISET)-FM(I,ISET)**2 
10  SO (I, I SET) = SORT(ABS(SOd,ISET) )) 

IF(IQ.NE. 23) PRINT*,"  PR  ", (PROTO (K) , K=i, PSPACE) 

IFCIQ. NE. 23) PRINT" ,"  FM  " , (FM( K, I SET) , K= 1 ,PSP ACE) 

IF(IQ.NE.23)  PRINT*, ” SD  ", < SD( K, ISFT) , K= 1, PSPAC E) 

11*11+1  $ IF(II.GT.l)  IC  = 23 
00  20  1=1, PSPACE 


20  PROTOd)  *PACK  ( FM  ( I,  ISET ) , SO ( I, ISET)  ) 
RETURN 
END 


PROGRAM  MAIN (INPUT , OUTPUT , T A PEI 0 , T AP  D.  0 , T A PE  31 , T AP Er f ) 

INTEGER  STRTCOL,STOPCOL,NSKIPS,SENTUUM,SlhTAPE,TAPftmP, 

X OSPAC  E(  IP)  ,OENORM,OOC,LTRSTRT,LTKSTOP,PPOTAPE,TOPPTOS,  SENT  MP 
LOGICAL  OVERIDE, BETWEEN, WINOOWS.PRNT  ALL  ,PRNTSUM, ADJUST 
REAL  WRA  T 10 » HRAT  TO 
INTEGER  LWIOTH,HWIOTH,SKWIDTH 
LWI0TH=2  $ HWIDTH=32  S SKWID7H=? 

STRTCOL5 1 $ STC.PC0L  = 133  $NSKIPS=1  S SrNTNUM=A  l SENTAPE=5  0 

COATA  NEXT  LINE  VARIES  WHICH  SENTENCE  USER  FOR  DATA 

SENTNUM=1 

OVERIDE=.T.  I 0EN0RM=1  $ OOC.=  0 5 LTF  STRT  = 1 $ L T R STD  P - 26 

PROTAPE=iO  $ PRNTALL=.T.  $ BETWEEN*. T.  * WINDOWSNT. 
NRTWT0P=5 

TOPRT  OS  = E $ T a P ET  N P=  3 1 $ AOJUST=.T.  $ WPATID=1.5 

HRATI0=1  .5$  PRNTSUM=.E.  % SENTEMP=16 
OATA  0S° ACE/ 10*8 1 / 

CALL  PREPSEN(SENTAPE) 

CALL  PREPRTO (PROTAPE) 

CALL  CMPRSEN(STRTCOL,STOPCOL,NPKIPS,SENTNUM,SENrAPE,OVEKlOE, 

XOSPACE,OENORM,OOC,LTRSTFT,LTRSTOP, P?OTAPF,PRNT ALL, P^NT3UM,Rfcr  WEEN, 
XMINDOWS,  TOPRTOS,TAPFTMP, AOJUST, WRATIO, HRAtIO,SENTEMP,lWIOTH 
X ,HWIDTH,SKWIDTH,NBTWTOP) 

STOP 

ENO 

SUBROUTINE  CMPRS EN ( STRTCOL , STOPCOL , NSKI PS , S ENTNU M, SENT  APE ,D VE RI"E , 
XOSPACE,OENORM,ODC,LTPSTRT  ,LTRSTOP,P^OT  APF ,PRNT  ALL,°RNTSUM,BFT WEEN, 
X W1 NOOWS,TOPPTOS,T A PETMP, AC JUST, WRATIO, HR  A ~ 1 0 , S ENT  t MP,  L W I DTH 
X ,HWIOTH ,SKWIOTH , N9TWT0F) 

INTEGER  STRTCOL»STOPCOL,NSKTPS,SENTNlJM,SEI.TAPE,TAPrTMP, 

X OSPAC E( in)  , OENORM,CCC,ITRSTkT,LTRSTOd,PRPTAPF,TOPPTDS,SENi  EMP 
LOGICAL  OVERIOE, PFTWFEN, WINDOWS, PRN1 ALL  , CRNTSUM, AO JUST 
REAL  WRATIO, HRATTO 

INTEGER  PRNTC0L»LTRMASK»SU3MASK,l)XHASK,NSFACE,EN0Rm,PC, 

X LSTPRT0,TOTPRTO,PROTNUN,NCOL,FLAG,I APENUM,KFY,PSPACF, 

X PE NORM, POC,START,STOF,PRTOCNT,LASTROW,I,J, XFRMCOL , NS PL  I TS, N°AG  F, 

X STRTROW,STOPROW,SUBLTR(5  2) ,LTR(52) , DX(r  2)  ,TEMP,LOW,HISH,SUnTNT 
X ,T  OPNUM,  LETTERS(26)  , ENDCOL  ,SINOEX  ( 134) 

X , WMASK  , WIDTH, LWIDTH, HWI DTH ,SKWIDTH,LRANGE,HRANGE 
COMMON /uETTERS /LET  TER,  Nt'MTRS,RE3N'J1' 

INTEGER  LFTTER (32, 32) ,NUMLTRS,RECN'M 
COMMON /RFSULTS/RFSULTS 
INTEGER  RFSULTS(52, 133) 

COMMON/PROFLTR/FLTRLTK,PROTO,RECFLTF 

INTEGER  RECFLTR 

REAL  FLTRLTR(18D  ,PROTO(180) 

COMMON /TABLES/  MTABLE, TABLE, ETABLE 

INTEGER  MTABLE (27,5) ,T ABLE(iO,lE) , ET ABL E ( 2fc U , 14) 

REAL  WINDOW(ieC) 

EOUIVALENCE ( WINOCW, FLTRLTR) 

5 FORMAT  (1HT) 

10  FORMAT (1H1) 

15  FORMAT (1H  , 136 ("* ”) ,/lH  , 24X , “RESULT S OF  ALGORITHM  TO  SOLVF  T Hf  ST 
XGMENTATI ON  PROBLEM  FOR  NON-I SOLATED  TEXT  CHARACTERS”/ 

X1H  , 136 ( "* ”) ,///) 

16  FORMAT (1H  , "SENT  FNCF  ”13”  COLUMN  ”13"  TO  "13”  EVALUATED  AGAINST”/ 
X“  PROTOTYPE  SET  ”12"  LETTERS  ”12"  TO  ”T2"  WITH  A LETT  ER-TO-WT  NDOW 
XWIOTH  RATIO  OF  "R».  ,2,  /”  AND  A L ETTF R-T 0 - WI NOD W HEIGHT  TATIO  OF  " 
XFL.2,/,”  SENTENCF  WINDOWS  RANGE  FROM  "12"  TO  ”12"  COLUMNS  WIDE" 

X"  WITH  SKIPS  IN  WINDOW  SITE  OF  ”12”  COLUMNS",/,”  AND  SKIPS  BETWEEN 
X WINDOWS  WITHIN  THE  SENTENCE  OF  ”12”  S3ACES”///> 

17  FORMAT  (1H  , "FINAL  RFSULTS  OF”) 

18  FORMAT  (1HC , “USER  REOUFSTEO  PRINTOUT  OF  THE  CLOSEST  "13"  PROn’YPE' 
v fop  the  ftmai  nrriFTnw  hatpty .*•/••  unwcwrr.  fhp  tmf  p^nimvc:.  Flit 


- - ....  ..i  - — - 


oonn 


-i-  ■«»»>« — rrir*  -r-irmi'w  1 . ifrrfxTf^tf  V ' ""T rLrVrL1 '»  vi"'ny”  *’tr,C»Tx  I ' • ■ , \»  t v > < • -.  . « «. 

X AND  PANGE  OF  LETTERS  SPECIFIED,  ONLY  "13"  PROTOTYPES  EXIST."/ 

X~  USER  REQUEST  FOR  THE  FINAL  DECISION  NAT f IX  RtOUCED  TO  THE  CLOSES 
XT  -13**  PROTOTYPES”//) 

20  FORMAT  (1HC  , "PROTOTYPES  IN  THE  MATRIX  ARE  I ”26 ( A1 , 1 1 , 1 X)  , / , 3 
X (1H  ,30X,26(A1,I1,1X)/)) 

25  FORMAT (1H  , "RESULTS  MATRIX  SHOWING  THE  PROTOTYPE-T O-WIN DOW  " 

X/"  DISTANCES  FOR  PROTOTYPES  WHOSE  WIDTH  VARIATION  ENABLES" 

X/"  THE  PROTOTYPE  TO  BE  COMPARED  TO  WINDOWS  OF  WIDTH  "12"  COLUMNS”) 

26  FORMAT (IH  , "MATRIX  IS  SORTED  BY  PROTOTYPE.  PROTOTYPE-TO-WT NOHW" 

X /"  DISTANCE  IS  SHOWN  AS  IT  VARIES  CONTINUOUSLY  OVER  THE  SENTtNCE" 
X /"  COLUMNS  USING  WINOOWS  OF  WIDTH  "12"  BEGINNING  AND  ENDING  A T " 
X/"  THE  COLUMNS  SHOWN  BY  THE  MATRIX  LEFTMOST  COLUMN") 

27  FORMAT C1H0.60X, "MATRIX  "12”  OF  ”12) 

30  FORMAT ( IH  , 1 0 ("*“),= >),/"  COLUMNS  * ", 

X =(3X,Al,Ii,4X),/,lH  , 1C  (”*”> , = (9( ))) 

32  FORMAT  (IH  ,I3"-"I3, = <F7.3,2X)) 

35  FORMAT (IH  ,I3"-"I3 = <A1 ,11,"-", F5. 3, IX) ) 

40  FORMAT  Cl H , "MATRIX  IS  SORTED  BY  DISTANCE  WITHIN  THE  SENTENCE" 

X/"  COLUMNS.  THE  "I?"  TO  "12"  CLOSEST  IDENTIFYING  PROTOTYPES  ARE" 
X/"  SHOWN  FOR  EACH  WINDOW  OF  WIDTH  "12"  PFC-INNING  AND  ENDING  XT" 

X/"  THE  SENTENCE  COLUMN  1NUICATEO  BY  THE  LEFTMOST  MATRIX  COLUMN.") 

41  FORMAT (IH  . "MERGED  RESULTS  MATPIX  SHOWING  THE  PROTOTVPE-TO-WI NDOW" 
X/”  DISTANCES  FOR  ALL  PRCTCTYEPS  WHOSE  RANCE  VARIATION  ENABLES  THE" 
X/"  PROTOTYPE  TO  RF  COMPARED  WITH  WINDOWS  OF  WIDTH  "12"  TO  "12 

X/  " COLUMNS.  THF  MATRIX  IS  SORTEO  BY  DISTANCE  WITHIN  THE  SENTENCE 
X“/"  COLUMNS  TO  SHOW  THE  TOP  "12"  PROTOTYPE  CHOICES  REGAROLESS  OF" 
X/"  WINDOW  SIZE.  ALL  WINDOWS  BEGIN  AT  THE  SAME  SENTENCE  COLUMNS" 
X/"  WHICH  ARE  SHOWN  IN  THE  LEFTMOST  MATRIX  COLUMN") 

42  FORMAT  (IH  , "MATRIX  OF  THE  "12"  TO  " 12"  CLOSEST  IDENTIFYING  * 

X "PROTO TY PES  FOR  WINOOWS  OF  WIDTH  "12"  SORT  ED  BY  DISTANCE") 

43  FORMAT  ( IH  11  ("*"),=  (9  ("*")),/,  IH  , "SENTENCE  1 ", 

X=("GUESS  "I  2, 1 X) , / , IH  , "COLUMNS  * "= ( "PPOTO/DX  "),/,lH  11  ("  ”), 

X = (9  ("*"))) 

45  FORMAT  <1H  2X,I3,2X"  * ",  = ( Ai , II,  F?  . 3,  IX) ) 

46  FORMAT  ( IH  11  ("*"),  = (16  ("*")>/,  IH  , "SENTENCE  * ", 

X =("  CHOICE  #"T2"  ") , / , 1 H , "COLUMNS  "= ( "PROT0-WI DTH- 0<  ")  , 

X / , IH  ,11  ("*-),  = (16  ("*")>) 

50  FORMAT (IH  , 2X, 13 , 6 X ,=  ( Al , II , 2X, 12 , 2X  , FS . 3 , 3 X) ) 

OAT  A L ETT  ERS/"  A"  , "B",**C" , "D“  ,"E",  "F" , **G",  *’H",  "I",  ” J",  **<",  "L  ”,  "M", 
X"N~,"0", ~P","0", "K","S","T","U","V",  "H",”X",”Y ",”?"/ 

LTRSIZE=  32  f L ENGTH  = 133 
PRNTCOL=13  $ NRNTCOL=7 
LTRMAS K=  MASK  (8  ) 

SUBMASK=SHIFT  (MASK  (8)  ,5  2) 

WHASK=SHIFT  ( MASK  (B)  ,4*.) 

OXMASK=SHTFT  (M  AS  K ( 2 2 ) , 22 ) 

WRITE  10 
WRITE  15 

CALL  OPENMS( SENT EMP, SI NT  EX, 134,0) 

CALL  GETSFN  ( 32, SENTNUM, SENT APE) 

CALL  PRTOTBL (L TRSTRT , LTRSTOP , TOTPRT J , P 2NTALL,PRNTSUM,.T. ,OVF?lDF, 

X OSPACE, OENORM,ODC,PROTAFE,OVERIOE,T  APSTMT) 

TAPENUM=  PROTAPE 

IF  (OVERIDF)  TAPFNUM=TAFETMP 

TOPNUM=TOPRTOS 

IF  (TOPNUM.LF.TOTPkTO)  GO  TO  80 
TOPNUM=TOTPRTO 
80  CONTINUE 

INITIALIZE  THE  TOPNUM  ROWS  OF  THE  RESULTS  ARRAY  TO  THE 
LARGEST  DISTANCE,  LTR=  7 AND  SMBLTR=f 

1 = 26 

T=SHIFT( 1,52) 

TEMPa2  00  0 0 0 0 


UU  I w • »w  v-t  - O MM  V UW  f W I SJ  % W W 9 II 

00  70  NR0W=1 »T OPNUM 
Rf SUITS (NR0W.NC0L) =I.0R.TEHP 
70  CONTINUE 

f 

) SORT  OF  THE  ETARLE  BY  W10TH  ANO  NS®ACE 


IF  (T0T9RT0 • EO  .1 ) GO  TO  123 
LSTPRT 0=T  OTPRT  0-1 
FL AG*0 

00  120  PR0TNUM=1 * LSTPRTC 

IF  (CETABLE(PR0TNUM,6) .GT.ET ABLE (PROTNUM ♦1,6))  .OR. 

X ( (ETABlE(PR0TNUM,6)  .EQ.EI  ABLE(PROTNUM>l,r  ) ) .ANO. 

X (ETABLE(PROTNtJM,7) .GE . ETARLE ( PROTNUM+1 , ~ ) ) ) ) GO  TD  120 
00  110  NC0L=1,14 
TEMP*ETARLE ( PROTNUM, NC0L ) 

ETABLE (PR0TNUM,NC0L) = ETA FIE ( PROTNUM* 1,  NCOL ) 

ETARLE (PR0TNUM*1,NC0L)=TEMP 

CONTINUE 

FLAG=1 

CONTINUE 

IF  (FLAG.EQ.l)  GO  TO  100 
CONTINUF 

BEGIN  LOOPING  ON  THE  PRCTYPES  CONTAINER  IN  THE  EXIST  TABLE 
PRODUCED  BY  PRTOTBL 

WIDTH=0 
OC=0 
NSPACE=P 
ENORM= 0 
PRT  OCNT=  0 
NL00KS=1 

DO  380  LHWI0TH=1,NL00KS 

NWIOTH  =H WIDTH- SK WIDTH* (LHHIDTH-1) 

DO  370  PROTNUM=l , TOTPRTO 
LRANGE*ETABLE( PROTNUM,* > 

HRANGE=ETABL  E ( FROTNUM ,5 ) 

IF  ( (NWIOTH. LT  .LRANGE)  .OR.  (NWIOTH. G1  .HRANCE)  > GO  TO  370 
KEY  = ET  ABLF ( PROTNUM , 1) 

PSPACE=E TABLE (PROTNUM, 7) 

PENORM=ETABLE (PROTNUM, 8) 

PDC=ET ABLE ( PRO TNUMjB) 

CALL  REA  OHS (TAPENUN , PROT 0 , PSPACE , KEY ) 

IF  THE  WIDTH  OF  THE  PROTOTYPE  CHANGES  aPINT0UT  THE  INTERMEDIA  F 
RESULTS  BEFORE  CONTINUING  WITH  THE  NEXT  PIOTOTYPE  AND  TKANSFD^M  T« 
SENTENCE  USING  THE  NEW  WINDOW  WIDTH 

IE  THE  ENERGY  NORMALIZATION  OR  THE  nc  TERM  IN  THE  NOPHA  LI  ZED 
WINOOWS  REQUIRES  A CHANGE  DUE  TO  A CHANGE  IN  THE  PROTOTYPE 
CHARACTERISTICS,  THE  ENTIRE  SENTENCF  MUST  BE  TRANSFORMED  AGAIN 

IF  ( (WIDTH. NE. NWIOTH) .OR . (Pt  NORM. NE. ENQRM) . OR. ( OC. NF. POC ) ) 

X CALL  FORTS EN ( LTF SI ZE, NWIOTH .ADJUST , WRAT 10 , HRA T 1 0, °SP®C  f , PE  NORM , 
XPDC ,STRTCOL,STOPCOL,NSK1PS,SENTNUM,SENTEMC) 

TF  ((WIDTH.NE. NWIOTH) .OR. (PENORM • N E. ENORM) .OR. ( DC. NE.  POC ) ) 
y.  NSPACEcPSPACE 

IE  THE  SPACE  OF  THE  PROTOTYPE  CHANGES  REDUCE  THE  SPACE  OF  THE 
WINOOWS  CONTAINED  IN  THF  TEMPORARY  SENTFNCE  FILE  OF  XFORMS 

IF  (NS PACE  .EO.PSPACE)  GO  TO  3*0 
00  340  XFRMCOL*STRTCOL,ENOCOL  ,NSKIPS 
CALL  R EA QMS ( SENT  FMP , WINDOW, NS PACE , XFRMC OL > 

CALL  REOUCE(WINOOW,NSFACE,WINOOW,PSPACE) 

CALL  WRIT  MS  (SENT  FMP,W  INC  OW, PSPACt ,XERMCOL> 
rnuTTunc 


ihi » 


» 1MUI 


UPDATE  ALL  VARIABLE  FOR  THE  CURRENT  PROTOTYPE 

IFMPRTOCNT.NE.O)  .AND. (WIDTH. NE.NWIDTH  ))  GO  TO  12r* 

CONTINUE 

«IOTH=NWIOTH  i NSPACE=PSPAC£  $ ENORH=PENOPH  $ OC=PDC 
PRT  OCNT=PRT  OCNT+1 
NROW=PRT  OCNT ♦TOPNUM 
ENOCOL=STOPCOL 

IF  ((ENDCOL+WIOTH-i)  .GT  .LENGTH)  EN3COL=L ENGTH- WI DIH+l 

STORE  THE  RESULTS  OF  THE  PROTOTYPE  AGAINST  ALL  THE  SENTENCE 
WINDOW  TRANSFORMS 

DO  360  XFPMCOL=STRTCOL,ENOCOL,NSKIPS 
CALL  REAOMS(SENT  EMP, WINDOW,N SPACE fXFRMQOL ) 

RESULT S( NROW . X FRMCOL ) =EUCLI D ( PROTO, WINDOW, NSPACE, DC)  *1010001 
I=SHIFT(ETABLE(PP0TNUM,2) ,52) 

J=SHIFT(ETABLE(PPOTNUM,3) ,4<  ) 

1=1. OR. J 

J=  SHI  FT  (WIDTH,  36) 

1=1. OR. J 

RESULT S(  NRO W, X FR MCOL > =1. OR. RESULTS (NPOH,XFRMCOL) 

CONTINUE 

IF  THERE  ARE  PLANK  COLUMNS  IN  THE  RESULTS  ARRAY,  SFT  THE  DISTANCE 
TO  2.0  ANO  THE  LETTER  TO  7 AND  THE  SUBLTR  TO  C 

1 = 26 

I =SHIFT (1,52) 

TEMP=2  CO  0 0 0 0 

DO  362  NCPL=XFRMCOl ,S , OPCOL , NSKIPS 
RESULT S(NROW,NCOL) =1.  K. TEMP 
CONTINUE 

IP  (PROTNUM.EQ.TOTPRTO)  GO  TO  125 
GO  TO  37  0 

CONTINUP  TO  LOOP  ON  THE  NEXT  PROTOTYPE 

STRTROW=  TOPNUM  *l 
STOPROW=STRTROW+PRT  OCNT-1 
LASTRO  W=  STOPRO  W- 1 

PRINTOUT  OF  THE  PROTOTYPES  OF  WIDTH  XX  AND  THF  WINDOWS  OF  THE 
SENTENCES  THEY  ARE  BEING  EVALUATED  AGAINST 

IF  (STRTROW.FO.STOPROW)  GO  TO  295 

SORT  BY  INCREASING  LETTER  PROTOTYPE  FOR  DISPLAY  OF  EACH  LETTERS 
PERFORMANCE  BY  SFNTENCE  COLUMN.  THIS  IS  THE  FIRST  OF  THE  OPT lONAi 
PRINTOUTS  FOR  THE  INTERMEDIATE  RESULTS 

FLAG=0 

OO  170  NROW=STRTPOW,LASTROW 

IF  ( (LTRMASK.AND. RESULTS (NROW) ) .LE . ( LT RMAS K. AND . RESULT  3 l NRO W* 1 ) ) > 

X GO  TO  170 
FLAG=1 

00  16U  XFRMCOL=STRTCOL,FNOCOL  ,NSKTPS 
TEMP=RESULTS (NROW,XFRMCOL) 

RESULT S( NROW f XFRMCOL )=RESULTS( NROW ♦l,XFRMCOL) 

RESULTS ( NROW *1 ,XFRMCOL)=TEMP 

CONTINUE 

CONTINUE 

IF  (FLA3.EQ.1)  GO  TO  15S 
DO  180  J=  STRTRO W, STOPRO W 

• TO  I \\  rSUTCT.U  XBJA/LS*  . A UD  - DC « t II  Kl  I.CTBTPftl  \ . 
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SUBLTR (J)=SHIFT(SU0MASK.AND. RESULTS ( J,STRTCOL) ,16) 

180  CONTINUE 
WRITE  10 

WRITE  25,  (WIDTH) 

WRITE  25,  (WIDTH) 

NSPLITS=  < PRTOCNT -1) /PkNTCGL*l 

START=STRTROW 

DO  210  NPAGEsl»NSPLITS 

ST0P=START4PRNTC0L-1 

IF  (NP AS E.EQ.N SPLITS)  ST0P=ST0PR0W 

SU8CNT=ST0P-START*1 

IF  (NPASE.GT.i)  WRITE  1C 

WRITE  2f  , (NPAGE»NSPLITS) 

WRITE  20,  ( LETTERS(LTR< J) ) ,SD8LTR(J) ,J=S1 ART, STOP) 

WRITE  30 , (SUBCNT , SUBCNT , (LETTERS (LTK  ( J ) ) , SUBLTR ( J) , 

*/.  J=START,  STOP)  , SUBCNT) 

DC  200  XFRMCOL=STRTCOL»  ENOCOL  ,NSKIC>S 
DO  190  J=START,STOP 

190  DX(J)=OXMASK.ANO. RESULTS (J  ,XFRHCSL) 

WRITE  32, (XFRMCOL.XFRMCOL+WIOTH-l, 

V.  SUBCNT,  (OX  ( J)  /iOOOOOQ.C,  J=START,STCP)  ) 

200  CONTINUE 

START=STOP  + l 
210  CONTINUE 
C 

C SORT  BY  DISTANCE  WITHIN  THE  SENTENCE  COLUPN  FOR  THE  SECOND  0J~PUT 

C OF  THE  I NTERMEDI ATE  RESULTS 

C 

CALL  S OPTIT(STRTCOL, ENOCOL, NSXIPS,STRr?OW,STOPROW) 

START=ST  RTROW 
TEMP=NBTWTOP 

IF  (NBTWTOP.GT.P«?TOCNT>  NBTWTOP=PRTCCNT 
NSPLIYS=  (NBTWTOP-l)  /PRNTCOL+1 
LOW  = l 
WRITE  10 

WRITE  25,  (WIDTH) 

WRITE  40,  (1,  SUBCNT,  WIDTH) 

WRITE  20,  (LETTERS  (LTR(J)  ) , SUBLTR  (J)  , J=STRT‘ROW,  ST  OPROW) 

OO  260  NPAGE=1,NSPLITS 
STOP=START4PRNTCOL-1 

IF  (NPAGE. EO  .N SPLITS)  STOP=STRTR0W  + NBTWTOP-l 
SUBCNT=STOP-START4l 
HIGH=LOW+SUBCNT-l 
IF  (NPAGE. GT.l)  WRITE  1C 
WRITE  2f  , (NPAGE,  NSPLITS) 

WRITE  42,  (LOW, HIGH, WIOTH) 

WRITE  43,  (SUBCNT, SUBCNT,  (NROW, NROW=L OW, HIGH) , SUBCNT, SUBCNT) 

DO  250  XFPMCOL=STRTCOL,ENDCOL  ,NSKI=>S 
OO  240  J-ST ART  ,STOP 

LTR(J) =SHIFT (LTRHASK.AND.RESULTS(J  ,XFRMCOL) ,6) 

SUBLTR (J)=SHIFT( SUBMASK. AND. RESULT S(J  , YFRMCOL ) ,16) 

DX  (J)  = DXM  ASK.  A NO.  RESULTS  (J  ,XFRMCOL) 

240  CONTINUE 

WRITE  35, (XFRMCOL,XFRMCOL+WI DTH-1, SUBCNT , (L ETT ERS (LTk ( J ) ) , 

X SUBLTR< J) , OX (J) /I 000000. 0,J=START,STOP) ) 

250  CONTINUE 

START=STOP*i 
LOW=HIGH4l 
260  CONTINUE 

NBTWTO P=  TEH  P 
295  CONTINUE 
370  CONTINUE 
380  CONTINUE 
RETURN 
END 

simpniiTT  mf  APNGrrx  fTFHP^NRnw«;.Nrni  s.  ta*P.  mihROWS  .n'IMCOI  S) 
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C DAT  A 
COAT  A 


C THIS 


2E 


C THIS 
50 


100 

110 


100 
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COMPLFX  TEMP(NRCWS,NCOLS) , TA MP (NUMRD WS , NIU'COLS ) 

INTEGER  NROWS, NCOLS, NUMROWS, NUMCOLS 
HO  100  NCDL=1, NUMCOLS 
00  100  NROW= 1 , NUMROWS 
TAMP(NRDW,NCOL)=T  EMP ( NRC W,NCOL) 

CONTINUE 

RETURN 

END 

SUBROUTINE  ENOEMFL  C FLT R LTR , NORMLTR , NS  PACE , DC ) 

RFAL  FLTPLTR  (NSPACE) , NORMLTR (NSPACE) 

INTEGER  NSPACE, DC 

INTEGER  START, NVECTOR  $ REAL  SJMSQRS 

ST  ART = 2 

IF  (OC.EQ.l)  START  =1 
SUMSQRS=0  .0 

OO  100  NVECTOR=START , NSPACE 

SUMSQR S=SUMS OR S* F LTRLTR ( NVECTOR)  **2 

SUMSQRS=  SUMSQRS*  * .5 

IP  (SUMSQRS.EO.O.C)  SUMSQRS=1.0 

00  110  NVECTOR=ST ART , NSPACE 

NORMLTR( NVECTOR) =FLTRLTR ( NVECTOR)  /SUMSQRS 

IF  (DC .EG. 0 ) NORMLTR(l)=l 

RETURN 

ENO 

REAL  FUNCTION  EUCL I D ( FLT R1 , F LTR2,  NSPACE , DD 
INTEGER  NSPACE ,CC 

REAL  FLTR1 (NSPACE) ,FLTR2(NSPACE) 

REAL  SUMSQRS  S INTEGER  START  $ LOGICAL  AMPSPEC 
TO  GET  COMPLEX  SPECTRUM  --  A MPSPEC=F AL SE 
TO  GET  AMPLITUDE  SPECTRUM  — AHPSPEC=TRUE 
AMPSPEC=. FALSE. 

SUMSQRS=  0.0 
START=2 

IF  (OC.EQ.l)  START=1 
IF  (AMPSPEC)  GO  TO  50 

BRANCH  FINDS  EUCL  EDI AN  DISTANCE  - COMPLEX  SPECTRUM 
90  25  NVECTOR=START, NSPACE 

SUM SQRS= SUMSQRS* (FLTR1(NVECT OR) -FLTR2 (NVECTOR) )* * 2 

CONTINUE 

GO  TO  110 

BRANCH  FINDS  EUCLEDIAN  DISTANCE  - A M°LITUD  E SPECT°JM 
PO  100  NVECTOR=START, NSPACE, 2 

AFLTR1=SQRT  ( (FLTR1  (NVECTOR)  ’**‘2)  *(FLTR1  (NV^CTOR+l)  * *2) ) 
AFLTR2=SOPT  ( (FLTR2  (NVECT  OR)*' *2)  ♦ ( FLT  °2  (NVECTOR  + 1)*  2) ) 
SUMSQRS  =SUMSQRS*  (AFL1R1-AFLTR2)  **2 
CONTINUE 

EUC  LID=SUMSQRS**  . 5 

RETURN 

FNO 

SUBROUTINE  FIL TER ( FORT LT R , NUMROWS , NUMCOLS , FLTRLT R, NSPAC E ) 
INTEGER  NUMROWS, NUMCOLS, NSPACE 
COMPLEX  FORTLTR(NUMROWS, NUMCOLS) 

REAL  FLTRLTR(NSPACE) 

INTEGER  NHRMNCS, NUMWRDS v RIGHT, L EFT , DOWN, NT OW,NCOL 

NHRMNCS= NSPACE** .5 

RIGHT=  NHRMNC  S/ 2* 1 

LEFT  = NUMCOL S-R IG HT  *2 

OOWN=RISHT 

FLTRLTR(l)=REAl .( FORTLTR ( 1 , 1)  ) 

IF  (NHRMNCS • EQ .1 ) RETURN 
NUMWRDS=  0 

00  100  NC0L=2, RIGHT 
NUMWRDS=NUMWPDS*2 

FLTRLTR(NUMWROS)=PEAL (FORTLTR (l,NCOL) ) 

FLTRLT  P(  NUMWPOS  *1  )=AIMG(FORTLTR(l  ,NC0L)) 
nn  ??n  gpnu=?.QnuM 


210 


i 

| 200 
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95 


100 

110 

120 

130 

140 

150 


00  210  NCOL=LEFT,NUMCOLS 
NUMWROS=NUMWROS*2 

FLTRLTP(NUMWRDS  ) =PEAL (FORI LTP ( NROW , NCOL ) ) 

FLTRLTR(NUMWRDS*1) = AIMAG ( FORTLTR ( NRP W, NCOl  ) ) 

CONTINUE 

no  200  NCOL=l»RIGHT 
NUMWR0S=NUMWR0S«-2 

FLTRLTR(  Nl'MWRDS  ) =F EAL ( FORT LTR (NROW , NCOL  ) ) 

FLTRLT  R( NUMWRO S+1)=AIHAG (FOR TLTR( NROW,  NCOL ) ) 

CONTINUE 

CONTINUF 

RETURN 

END 

SUBROUTINE  FORTS FN (LTkSI ZE , WIDTH, AD JUST , WF AT  ID , HRATI0 , N SP A3 F , 

X ENORM,OC  ,STRTCOL  .STOPCGL  ,NSKIPS , S ENTNUM , ‘ APENUM) 

INTEGER  LTRSIZE, WIDTH, NS PACE»STRTC9L»ST0pCOL,NSKIPS,SENTNUM, 

X TAPENU*,ENORMtDC 

LOGICAL  ADJUST  $ PEAL  WRAT 10, HRATIO 

INTEGER  RECFLTR 

REAL  FLTRLTR  180 ), PROTO ( 180) 

CO MMON/°RPF LTR /FLTRLTR, PROTO, RECFLTR 
COM MON /X  FORMS/FORT  L TR 
COMPLEX  FORTLTRC  fi4 ,64) 

INTEGER  ROW  * COL, NROW , NCOL , ENDCQL , XpRMCOL , L ENGTH 

LENGTH=1 33 

ENDCOL=STOPCOL 

IF  ( (ENDCOL  + WIOTH-1)  .GT. LENGTH)  EN  DC  OL  = L ENGTH-WI OT  H*1 
00  130  XFRMCOL=STRTCOL, ENDCOL  ,NSKTPS 
CALL  MIDDLE(LTRSIZE»XFRWCOL,WlOTH) 

CALL  XFORMIT (L TRS I ZE , NI TTH , HEIGHT , AO JUST , WRAT 10 , HRAT 1 0,  NS °A CE 
X ,NUMRO  WS»NUMCOLS,0, 0,.F.,.F.,.F.,.F. , SENT HUM) 

CALL  FILTER (FORT LTR , NUMF OWS , NUMCOLS, FLTRLTR»NSPAC£) 

IF ( (ENORM.EO.l).CR.  (OC.FQ.l) ) CALL  ENOR MFL(FLTRLTR, FLTRLTR, N 3 FACE, 
X DC) 

CALL  WRITMS (TAPENUM  , FLTRLTR, NS°ACF, XFRMCOL) 

CONTINUE 

RETURN 

END 

SUBROUTINE  GETFLTR  (ALPHNUM, LTRNUM, L SDACE , NSPf CE , T A PF NU M , RA NO  DM ) 

INTEGER  ALPHNUM, LTRNUM, L SPACE, NSP ACE , T A PENUM 

LOGICAL  RANOOM 

INTEGER  RECFLTR 

PEAL  FLTRLTR(180) ,PROTO(180) 

COMMON /PRO ELTR /FLTRLTR, FROTO , RECFLTR 
INTEGER  RECRQO,RFCSKIP,KEY 
DATA  RE3ELTR/1/ 

NUMWROS=LSPACE+4 

IF  (. NOT. RANDOM)  GO  TO  95 

KEYr (ALPHNUM-l)*2f ♦LTFNUM 

CALL  REA DMS (TAPENUM, FLTFLTK. NUMWRDS, KEY) 

GO  TO  15  0 
CONTINUE 

RECRQO=( ALPHNUM-1)*  2E+LTRNUM 
R EC SKI P=IABS (RECROO-RECFLTR) 

IF  (RECRQO-RECFLTP)  100,140,120 
DO  110  1=1, RECSKTP 
RACKSPA3E  TAPENUM 
GO  TO  14  0 

DO  130  I=1,RECSKIP 
READ  (TAPENUM) 

READ  (TAPFNUM)  ( FLTRLTR  ( I ) , I = 1, NUMW?  ITS ) 

RECFLTR=RECROD*l 

CALL  REDUCE ( FL  TRL TP, L SPA C F, FLTRL  TR,NSPACF) 

RETURN 

CNO 

— i | ji_ 


AD-A069  298  AIR  FORCE  INST  OF  TECH  WRIGHT-PATTERSON  AFB  OHIO  SCH~ETC  F/6  9/2 
• SEGMENTATION  OF  TOUCHING  ENGLISH  LETTERS. (U) 

MAR  79  R E BENTKOMSKI 

UNCLASSIFIED  AFIT/6E/EE/79-1  NL 


■JWJIVUV  I x I 


Oil  JL  |1  \ L • CA  O X 4 L)  Jtl 


1 1V>  If  I ►*  I U'lVM  1/ 


INTEGER  LENGTH, LTRSIZE,SENTNUM,TAPENUM 
COMMON /S ENT NCE/S AMPLE, HOLOIT 

INTEGER  SAMPLE  (3  2, 133) , HOLOIT (254),  NUMCHAF, TITLE  (3D,!.)  , 

V.  HOLDS EN  (133) 

EQUIVALENCE  (HOLOIT (1) ,NUMCHAR) , ( HOL 01 T ( 2 ) , TITLE (1 , 1) ) , 

7.  (HOLD  IT  (122),H0LDSEN(1)  ) 

INTEGER  TEMP 
LENGTH  =133 
NUMHRD  S=LENGTH+12(!  *1 
HASKER=MASK (1) 

CALL  REA  OMS (TAPENUM, HOLOIT, NUMWRDS, SENTNUM) 

00  *130  NC0L  = 1,  133 
TEMP=HOLDSEN (NCOL ) 

00  130  NR0W=i,LTRSI7E 
IF  (TEMP. ANO. MASKER)  110,100 
SAMPLE  (NROW, NCOL)  =0 
GO  TO  120 

SAMPLE (NR0W,NC0L)=1 
TEMP=SHI  FT  (TEMP,  1) 

CONTINUE 

RETURN 

END 

SUBROUTINE  MIOOLE(LTRSI7E,XFRMCOL,HIOTH) 

INTEGER  LENGTH, LTRSIZE, XFRMCOL, WIDTH 
COMMON/LETT ERS/L ETT ER ,NUMLTRS ,RECNUM 
INTEGER  LETTER (32, 32) , NUMLTRS.RECNUM 
COMMON/SENTNCE/SA  MPL  E, HCLOIT 

INTEGER  SAMPLE (3  2,133) ,HOLOI T (254 ) , NUMCHAF, TITL E ( 3 n , i ) , 

7,  HOLDS  EN  (133) 

EQUIVALENCE  (HOLOIT  (1)  ,NUMCHAR)  , (HO.OIT (2)  ,TITLE<1,1)>  , 
r (HOLOIT (122), HOLOSEN(l)) 

INTEGER  ST  OP, NROW,  NCOL,  L EDGE,  REOGE,  T EM** 

FORMAT (1H0, "WIDTH  ( *'1 3**)  PLUS  START  COLUMN  C*I3**)  EXCEEDS  TH;  SEN’ 
XENCE  LENGTH  OF  "1 3"  COLUMNS") 

LENGTH=133 

STOP=XFRMCOL  +WI0TH-1 
IF  (STOP. LE. LENGTH)  GO  TO  10  0 
WRITE  1, (WIDTH, XFRMCOL, LENGTH) 

RETURN 

LEDGE=  (LTRSI7E-WIDTH)  /2 
REOGE=LEOGE+WIOTH*i 
DO  140  NROW=l, LTRSIZE 
DO  110  NCOL=l, LEDGE 
LETTER (NROW, NCOL ) =0 
TEMP=L EDGE 

DO  120  NCOL=XFRMCOL » STOP 
TEMP=T  EMP  + l 

LETTER (NROW, TEMP) =SAMPLE (NROW, NCOL ) 

IF  (REDGE.GT .LTRSIZE)  GO  TO  140 

00  130  NCOL=REOGE, LTRSIZE 

LETTER(NROW,NCOL)=0 

CONTINUE 

RETURN 

END 

SUBROUTINE  PREPRTO ( PROTAFE) 

INTEGER  PROTAPE 
INTEGER  PINOEX  (28  8 ) 

CALL  OPENMS (PROTAPE, PINOEX, 288,0) 

RETURN 

END 

SUBROUTINE  PREPSFN ( SENTA PE) 

INTEGER  SENTAPE 
INTEGER  SINDEX(2*4) 

CALL  OPENMS (SENT APE, SINOEX, 244,0) 

RETURN 
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SUBROUTINE  PRTOTRL  (L  TRS T h T , ITRST OP,  T 0 TPR"  0 , PRINTI  T , SUMMARY  , 
XSKIPAGE, 

XOVFRIOE,OSPACE,OFNORN,OPC,TAPENUi.NrWJ  APE,TAPFTMP> 

INTEGER  LTRSTPT,  l TRSTGP, T OTPPTO. T APT  NUM,  T A PETMP,  OS*»ACE < 1 0 ) 

X , OENORM,  OOC 

LOGICAL  PRINTIT,SKIPAGE,OVERIOE,HEWI  APE, NUMMARY 
INTEGER  IFTTFRS<26> , MASKER, LTRNUM, KEY, SUBI  TRS.PROTNUM, 

X NSPAC  F,  P^PACE  ,NCOL»NPOW,NALPHS.NTEiP, ALPMNUM,ALPH^TR,  ALPHST5, 

X ALPH( 150) ,ENORMt0C,NlNPEX<300) 

INTEGER  NT ABLF<27,5> , TABLE (15,10)  ,Et  ABLE  ( rt>0,lw) 

COMMON /! APLFS/  MT A Bl E , T A 01 E , ET ARL E 

INTEGFR  RECFLTR 

PEAL  FLTRLTR <180 ) ,PROTO<180> 

COMMON/PPOFl T R/FLTRLTR, PROTO, RE CFLTi 
FORMAT  < 1M0 ) 

FORMAT  <1H1> 

FORMAT <1H  ,41X,5M"*"),/,1H  ,41X , “OR ARACT rRI ST IC S OF  PROTOT  Y»  F FIL 
XE  ”12”,  LFTTFRS  "12"  TO  ”I2,/,1H  ,^1  X,  55  <”*">> 

FORMAT  <1H0, *»?<***  ”)  / 1H  ,“LTR  SURLTP  FANGE  •'PACE  ENORM  03  T ERM  *AIPM 
X "/  1H  ,42  <”*”>) 

FORMAT  <1H  , "CLASS  #”12“,  LtTfER  ”A1”  **  DOFS  NOT  EXIST*  ** ”» 

FORMAT  (*<”*”)> 

FORMAT  <1H  ,MCLASF  »"l 2",  LETTER  “AD 
FORMAT  <1  HO,  “TOTAL  # OF  ALPHA  BETS=”I3 
X ”,  # OF  (POSSIBLE)  SUBCl ASSF  S= ”1 2 , /**  RANGE 

X”SUBCL  ASSFS  IN  ASCFNOING  ORDER  IS  » ”T  2, ”-”I2 , 9 (” , ”, le » 

FORMAT <1H  OF  (ACTUAL)  PROTOTYPES  ="I2”.  CHARACTERISTICS 

X“THESE  PROTOTYPES  APE  SHOWN  BELOW.”/  ) 

FORMAT (1H  ,"LTR  SUBLTR  RANGE  SPACE  F NORM  "C  TERM  #AlPM  ALPHABETS 
XIN  THE  PROTOTYPE”) 

FORMAT  < 1 H ,1X,  A1 , 3X  , 1 2 , 3X  , 1 2 , 1 2, 3X  , T 3 , 3X  , A 3 , 4X  , A 3, 4X  , 1 3 , 3 X , 

X 22<I3, ”,”)/, 6 <1H  ,44X, 22(13, ”,”)/>> 

FORMAT ( 1 H , IX, A1 , 3X, 12, 3X , 12 , I 2, 3X, I 3, 3X, A 3, AX , » 3,  AX , 13, 3 X) 
FORMAT  <1H  ,42  <”-”)) 

FORMAT  <1H  ,1X, A1 , 3X , 1 2, 3X  , 1 2 1 2 , 7X  , 1 3 , 3X  , A 3, 4X  , A 3, 4 X , 1 3, 3 X , 

X “THIS  SUBLTR  OOFS  NOT  EXIST,  NO  ALPHABETS  IN  THE  RANGE”) 

DATA  L ET T TRS /” A” ,” B” , ”C“ , “0” , "E” , "F”, "G ” , ”H " , ”1 ” , ” J” , ”K ” , "l ”, ”M” , 
X**N”,”0”,  ”P”,  ”0”,  ”P“,"S“,  “T”,  "U”,  ”V“,  ”W”,”X“,"Y  ",”7”/ 

TOTPRT  0=0 

IF  (OVERIOF.ANO.NFWTAPE)  CALL  OPENMS <T APFT MP ,N1N01X ,2 88 , 0 ) 

PALL  REAOMS (TAPE NUM.MTAPLb, 135,1) 

MASKER =M ASK  < 1 ) 

IF  (PR INTIT . AND. SUMMARY ) PRINTIT=.E. 

160  LTRNUM*LTPSTRT,LTRST()P 


or ( possi plF) 


- <>«? 


DO 

IF 

IF 

GO 


80 


( .NOT .PRINTIT ) GO  TO  85 
(NT ABLE (LTRNUM, 1 ) . EQ. 1)  GO  TO  80 
TO  150 
CONTINUE 

IF  (MTABLF (LTRNUM, 1) .EQ. 0)  GO  TO  160 
KEY»MTABLE (LTRNUM, 2) 

SUBLTRS=  MTABLF  (LTRNUM,  3) 

NUMWRDS«  SURL  TRS*1 5 

CALL  R EA  DMS ( TAPENUM , TABLE , NUMWRDS , KF  Y> 
DO  150  PROTNUM«l,SUBLTRS 
IF  (T ARLF ( 1 , PROT NUM) . ED • 0)  GO  T 3 1 ? 

T OTPRT 0* TOTPRT 0 »1 
IF  (.NOT.  OVFRIDF)  GO  TO  100 
PS«»ACE=TAPLF(8  ,PROTNUM) 

TABLE! 8, PROT NUM) =NSPACF=OSPACE(PROTNUM> 

T ABLE  ( 9,  PPOTNUM)  *OENOKM 

TABLE(10,PROTNUM)=OOC 

IF  ( .NOT .NEWTAPF)  GO  TO  100 

KFY«TARLF(2,PR0TNUM) 

CALL  REA  DMS (TAPE NUM, PROTO, PS PACE,  KEY ) 
CALL  REDUCF  (PROTO,  PSP  AC  L , PROTO, NS PA  CF ) 

i'.A(  I WRTTMSXrAB£lMP.PILDTn.NSPAr.F.KrY* 


10  Q CONTINUE 

00  110  NCOL * 2,  15 

110  ETA8LE  (TOTPRTO  , NC0L-1)=TA8LE(NC0L, PROTNUM) 

IF  C SUMMARY ) GO  TO  142 
IF  (.NOT.PRINTIT)  GO  TO  IE  0 
NALPHS*0 

00  140  NWR0S=1,3 

NT EMP* TABLE ( 12 +N MR DS, PROTNUM ) 

ALPHST  R=  (NWROS-l)*60d 
ALPHST  P*  Al  PHSTR*E  9 
IF  (ALPHSTP.GT.150)  ALPHSTP=150 
00  130  ALPHNUM*  ALPHSTR, ALPHSTP 
IF  (NTEMP.ANO. MASKER)  125,127 
125  NALPHS=NALPHSM 

ALPH(NALPHS) =ALPHNUM 
127  NTEMP*  SHIFT  (NT  EMP , 1) 

130  CONTINUE 
14  0 CONTINUE 

142  NCOL*T  A3LE  <12»  PROTNUM) 

EXORM*  DO  a "OUT" 

IF(TA8LE(9, PROTNUM) ,EQ.l)  ENORM="IM  " 

IF(TABLE(10, PROTNUM) .EO.l)  OC*"IN  " 

IF  (TABLEd,  PROTNUM)  .60.0)  145,147 
145  CONTINUE 

IF  (SUMMARY)  150,149 
147  IF  (.NOT. SUMMARY)  GO  TO  148 
GO  TO  13  0 
14  8 CONTINUE 

149  CONTINUE 

150  CONTINUE 
16 0 CONTINUE 

RETURN 

ENO 

SUBROUTINE  REOUCE  ( INFLTR , INSPACE ,OUTFLTR , OTSPAC E) 

INTEGER  INSPACE, OTSPACE 

Rf  «L  INFLTR (INS FACE) ,OUTFLTR (OTSPACE) 

If  EGER  INHRMNC, OUTHRMC, SKIP, OOWN, OUT VCTR , IN VC TR, START, STOP 
IF  (INSPACE.NE.OTSPACE)  GO  TO  90 
0<~  85  NVECT0R*1, OTSPACE 
85  Ot  ’TLTR(NVECTOR)  =INFLTR  (NVECTOR) 

G r TO  125 
90  CONTINUE 

I*'MRMNC=  I NS  PACE"  * .5  1 OUTHRMC=OTSP  AC  E * ,5 

S*  TPa ( INHRMNC” OUTHRMC ) * 2 
OOWN»OUTHPMC/2M 
OUTFLTR( 1) =INFLTP (1) 

IF  (OTSPACE. EO.l)  RETURN 
00  100  DUTVCTR=2, OUTHRMC 
100  OUTFLTR( OUT VCTR) = INFLTR ( OUT VCTR) 

INVCTR=START  = OUTHRMCd 
DO  120  NROH=2,OOWN 
INVCTR=I NVCTR*SKT  P 
STOP*STAPT+OUTHRMC*  2-1 
00  110  0 UT VC TRaS TART, STOP 
OUTFLT  P( OUT VCT  R) =TNFL1 R (INVCTR) 

INVCTR  = INVCTRM 
110  CONTINUE 

START* ST OP+l 
120  CONTINUE 
125  CONTINUE 

DO  130  1*1,4 

130  OUTFLTR(OTSPACEd)  =INFLTR  (INSPACE*I) 

RETURN 
• FNO 

SUBROUTINE  SI7EIT (LTRSI7E , MI DTH, HEIGHT > 
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COMMON /t  ETT  ERS /L  CTTER , NUMLTkS .RECNUH 
INTEGER  L ETT  ER 13  2,  32)  , NUMLTRS,  RECN'IM 
INTEGER  NCOL, NROW, STRTCGL, TEMP, STRT50W 
HO  100  NCOL*l,32 
00  100  NR0W*1,32 

IF  (LETT  ER  ( NRO  W,  NCOL ) .EQ.l)  GOTO  110 
100  CONTINUE 
110  STRTCOL*  NCOL 

00  120  NCOL  ■ 1,32 
TEMP*3  3-NCOL 
OO  120  NROW»l,32 

IF  (LETTER  (NROW,  TEHP).EQ.1>  GO  TO  130 
120  CONTINUE 
130  WIOTH* TEMP-STRTCOL  *1 

IF  (WIDTH. GE.l)  GO  TO  13? 

WIDTH* HEIGHT  =0 
RETURN 

135  CONTINUE 

00  1*0  NRO W* 1,32 
00  140  NCOL  * 1 , LT  RSI 7E 
If (LETTER(NROW,NCOL) .EQ.l)  GO  TO  15 C 
140  CONTINUE 
150  STRTROWsNROW 

00  160  NR0W*1,LTRSI?E 
TEMP*LTRSI7E ♦1-NROW 
00  160  NC0L=1,LTRSI?E 

IF  (LETT ER(TEMP, NCOL)  .EO.l)  GOTO  170 
160  CONTINUE 
170  HEIGHT *T  EMP -STRTROW  + 1 

RETURN 
END 

SUBROUTINE  SORT IT ( STRT COL ,ST OPCOL , MSKr°S, STKTROW, STDPRJ W) 

INTEGER  STRT COL, ST OPCOL, NSKI PS, STRT? OW,STOPROW 
COHHON /RESULTS /RESULTS 
INTEGER  RESULTS! r 2 , 13 3) 

INTEGER  FLAG,XFRMCOL  ,TEM P,L ASTROW, DXM ASK 
OXMASK  =S  HIFT (MASK ( 22) , 22 ) 

LASTRO W*ST0PR0W-1 

00  120  XFRMCOL*STRTCOL,STOPCOL,NSKIBS 
100  FLAG*0 

00  110  NROW* STRT ROW, LAST  ROW 

If  ( (OXMASK. ANO. RESULTS < NROW, XFRHCOL)  ) . L E . ( OXMASK.  MO. RESULTS  (NROW 
x*i,xfrh:od  ) ) GO  TO  110 
TEMP*RESULTS(NROW,XFRMCOL) 

RESULT S( NROW, XFRMCOL ) =RE SULTS (NROW+1 ,XFRMrOL) 

RESULTS(NROWM,XFRMCOL)=TEMP 

FL AG*1 

110  CONTINUE 

IF  (FLAG. EQ.l)  GO  TO  100 
120  CONTINUE 
RETURN 
END 

SUBROUTINE  X FORM IT (LTRS17E, WIDTH, HEIGHT, AT JUST, WRATI 0, HRATT D, 
XNSPACE,NUHROWS,NUHCOLS,ALPHNUM,LTRNJH,  VEXrANO,XFORM, I MAGE, INVERSE, 
XSENTNUM) 

INTEGER  LTRSI7E, WIOTH,Hf IGHT , NSPACE, NUMROWS , NUMCOL S, ALPHNUM , L T RNU  1 
X ,SENTNUM 

LOGICAL  ADJUST  $ REAL  WRAT IO,HRATI 0 
LOGICAL  VEXPANO, XFORM , IN AGE , INVERSE 
INTEGER  LETTER (3 2, 32) , NUMLTRS ,RECNUM 
COMMON/*.  ETT  ERS /LET  TER,  NUMLTRS,  RECNU" 

INTEGER  RFCFLTR 

PEAL  FLTPLTR  (180  ) , PROTO (180) 

COMMON /PROFLTR /Ft TRlTk,FROTO, RFCFLTR 

COMMON/X FORMS/ FOR TLTR 




COMPLEX  WORKER  (f-'.) 

INTEGER  POWERS  (7) 

INTEGER  NN(2) , NU MCOLS , NUMROWS, NC3L , NROW , FLAG, TOP, SOTTON  , ITEW=> , 
X COL, ROW,  LEFT, RIGHT ,IFORM,ISIGN 

DIMENSION  TRPLTR  (32,  32) 

DATA  POWERS/2,4, 8,16,32, 64,128/ 

FL AG*0 

ICO  CONTINUE 

CALL  SI7EIT(LTRS17E,WIOTH,HEIGHT) 

CALL  TPPMLT (LETTFP,TRPLTR,SENTNUW) 

NUMROW  S*  NUMC  OL  S=  3 2 
NRINGS*MSPACE**.S 

IF (NUMROWS.LT. NR TNGS)  NUMR0WS=NRI NGS 
IF ( NUMCOLS.LT. NR  I NGS)  NUMCOL S*NRINGS 
IF  (M00( NUMROWS, 2)  .EQ.i)  NUMR0WS=NJMR0 WS ♦ ! 

IF  (M00(NUMC0LS,2> .EQ.I)  NUMC01S=NUMC0LSM 
TOP  *1 ABS (NUMROWS -LTRSI7E)/2 
IF  (NUMROWS. LT.LTRSI7E)  GO  TO  ISO 
BOTTOM *T OP +1 ♦LTRSI ZE 
ITEMP*  TOP 

00  140  NC0L*1,LTRSI7E 
TOP»ITEMP 
00  110  NROW* 1 , TOP 
110  FORTLTR(NROW,NCOL)  =0 
00  120  NR0W=1,LTRSI7E 
T0P*T0P*1 

120  F0RTLTR(T0P,NC0L ) =TRPLTF (NROW,NCOL) 

00  130  NROW* BOTTOM, NUMROWS 
130  FORTLTR(NROW,NCOL>=0 

140  CONTINUE 
GO  TO  17  0 
150  CONTINUE 
T0P«T0P*1 

BOTTOM=T  OP ♦NUMROWS -1 
DO  160  NC0L=1,LTRSIZE 
ROW«0 

00  160  NROW* TOP, BOTTOM 
ROWsRO  W>1 

160  FORTLTR(ROW,NCOL)=TRPLTF (NROW,NCOL) 

170  CONTINUE 

LEFT*I ABS ( NUMCOL S-LTRSI7E)/2 
IF ( NUMCOLS .LT • LTRSI7E)  GO  TO  220 
FIGHT* LEFT ♦LTRSI 7E+1 
OC  210  NROW* 1, NUMROWS 
00  180  NCOL* RIGHT , NUMC 01 S 
180  FORTLTR(NROW,NCOl)=0 
ITEMP*  RIGHT 
00  190  NCOL* 1 , LTRSI ZE 
ITEMP* IT  EMP-1 
COL*LTRSIZE-NCOL+1 

190  FORTLTR( NROW, I TEMP) =FORTLTR( NROW, COL) 

00  200  NC0L*1,LEFT 
200  FORTLTR( NROW, NCOL) *0 
210  CONTINUE 
GO  TO  232 
2 20  CONTINUE 

LEFT*LFFT+1 
RIGHT*LEFT4NUMC0LS-1 
00  230  NR0W*1, NUMROWS 
COL*0 

00  230  NCOL*LEFT, RIGHT 
COL*COL* 1 

2 30  FORTLTP( NROW ,COL  > *FORTLTR (NROW, NCOL) 

2 32  CONTINUE 

CALL  ARNGCPX(F0RTLTR,64,64,F0RTLTR,Nim0WS,NU*C0LS) 

MMmm 
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240 

245 


250 

C 

C 

C 

' C 
255 
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IF  (FLAG.GT.O)  GO  TO  2^1 

NN(1)  = NUMR0WS  S NN  ( 2)  =NUMC0LS  $N0IM=2  % TSIGN=-1  3 IFDRM=P 
COL =R0  W=  0 
00  240  ITEMP=1,7 

IF  (NUMROWS.  EQ  .POWERS  (ITEMP)  ) RDW=1 
IF ( NUMCOLS .EQ. POWERS (ITEHP) ) C0L=1 

CONTINUE 

IF  ( ROW. EG. 1. AND.  COL.EO.l)  GO  TO  2F0 

CALL  FOURT  ( FORTLTR , NN, NDIM, ISIGN, IFORM, WORKER) 

GO  TO  25  5 

CALL  FOJP.T  <F0RTLTR,NN,ND1M, ISIGN, IFORM, 0.0) 

THIS  IS  THE  PLACE  TO  PUT  IN  AN  INVERSE  PRINT  ROUTINE  THAT 
FILTERS  THE  FOURIER  TRANSFORM  FROM  THE  MIDOLE  AND  THEN  CALL  SlPPT 

CONTINUE 

IF  (FLAS.GT • 0)  RETURN 

flag*flag*i 

RETURN 

END 

SUBROUTINE  FOURT (DATA ,NN , NOI M, ISIGN, IFORM, WORK ) 

THE  COOlEY-TUKEY  FAST  FOURIER  TRANSFORM  IN  USASI  SASIC  FORT 5 A N 

TRANSF  0RM(Ki,K2»  . ..)  = SUM  (DATA  ( J1 , J 2,  . ..)4EX®  (ISIGN*  2*  PI4SQU  (-1) 
* ((J1-1)4(K1-1)/NN<1)+<J2-1)4(K2-1)/NN(2)+...))) , SUMMED  FOP  ALL 
'Ji,  Ki  BETWEEN  1 ANC  NN(1),  J2,  <2  BETWEEN  1 AND  NN(2),  ETC. 

THERE  IS  NO  LIMIT  TO  THE  NUMBER  OF  SUBSCRIPTS.  DATA  IS  A 
MULTIDIMENSIONAL  COMPLEX  ARRAY  WHOSE  REAL  AND  IMAGINARY 
PARTS  ARE  ADJACENT  IN  STORAGE,  S'JCH  AS  FOE  TRAN  IV  PLACES  THEM. 

IF  ALL  IMAGINARY  PARTS  ARE  ZERO  (DATA  ApE  DISGUISED  REAL),  SE: 
IFORM  TO  ZERO  TO  CUT  THE  FUNNING  TIME  BY  UP  TO  FORTY  PERCENT. 
OTHERWISE,  IFORM  = *1.  THE  LENGTHS  QC  ALL  DIMENSIONS  ARE 
STORED  IN  ARRAY  NN,  OF  LENGTH  NOIM.  THEY  MAY  BE  ANY  POSITIVE 
INTEGERS,  THO  THE  PROGRAM  RUNS  FASTFR  ON  COMPOSITE  INTEGERS,  AND 
ESPECIALLY  FAST  ON  NUMBERS  RICH  IN  FACTORS  OF  TWO.  ISIGN  IS  +1 
OR  -1.  IF  A -1  TRANSFORM  IS  FOLLOWED  BY  t +1  ONE  (OR  A *1 
BY  A -1)  THE  ORIGINAL  DATA  REAPPEAR,  MULTIPLIED  BY  NTOT  (=NN  ( 1)  * 

NN  ( 2)  * ...)  • TRANSFORM  VALUES  ARE  ALWAYS  COMPLEX,  AND  ARE  REIUPNEO 
IN  ARRAY  DATA,  REPLACING  THE  INPUT.  IN  ADDITION,  IF  ALL 
DIMENSIONS  ARE  NOT  POWERS  OF  TWO,  A 5 RAY  WORK  MUST  BE  SUPPLIED, 
COMPLEX  OF  LENGTH  EQUAL  TO  THE  LARGEST  NON  2**  K DIMENSION. 
OTHERWISE,  REPLACE  WORK  BY  ZERO  IN  THE  CALLING  SEQUENCE. 

NORMAL  FORTRAN  OATA  ORDERING  IS  EXPECTED,  FIRST  SUBSCRIPT  VARYING 
FASTEST.  ALL  SUBSCRIPTS  BEGIN  AT  ONE. 

RUNNING  TIME  IS  MUCH  SHORTER  THAN  T«E  NAIVE  NTOT4*  2,  BEING 

GIVEN  BY  THE  FOLLOWING  FORMULA.  DECOMPOSE  NTOT  IN^O 

2**K2  4 3*»K3  * 544K5  * ....  LET  SJM2  * ?*K  2,  SUMF  = 34K3  ♦ 5*  K? 

...  AND  NF  = K?  ♦ K5  ♦ ....  THE  TIME  TAKEN  BY  A MULTI- 
DIMENSIONAL TRANSFORM  ON  THESE  NTOT  DATA  IS  T * TC  ♦ NTCTMT1  + 
T2*SUM2«-T3*SUMF+TA*NF)  . ON  THE  CDC  3300  (FLOATING  POINT  ADD  TIME 
OF  SIX  MICROSECONOS) , T = 3000  ♦ NTOT4 (5  Oi +4  34SUM2+F64SUMF* 

3204NF ) MICROSECONDS  ON  COMPLEX  DATA.  IN  AOOITION,  THE 
ACCURACY  IS  GREATLY  IMPROVED,  AS  THF  RMS  RELATIVE  FPPOR  IS 
SOUNOED  BY  34 244  (-B) 4SUM (FACTORt J) 4 * 1.5 ) , WHERE  B IS  THE  NUM3E& 

OF  BITS  IN  THE  FLOATING  POINT  FRACTION  ANP  FACTOR ( J)  ARE  THF 
PRIME  FACTORS  OF  NTOT. 

PROGRAM  BY  NORMAN  BRENNER  FROM  THE  BASIC  PROGRAM  BY  CHARLES 
RAOER.  RALPH  ALTER  SUGGESTED  THE  IDEA  FOR  THE  OIGIT  REVERSE. 

MIT  LINCOLN  LABORATORY,  AUGUST  19*7.  THIS  IS  THE  FASTEST  AND  MOST 
VERSATILE  VERSION  OF  THE  FFT  KNOWN  TO  THE  AUTHOR.  SHORTER  PRO- 
GRAMS FDUR1  AN0  FOUR2  RESTRICT  DIMENSION  LENGTHS  TO  POWERS  0s  TWO. 
— fcgc  ******  (JUNE  1967),  SPECIAL  ISSUE  ON  -FT. 


rME  DISCRETE  FOURIER  TRANSFORM  places  THKr  e rf  st  kic  t tons  upon  thf 

DATA, 

1.  THF  NUMBER  OP  INPUT  DATA  AND  THf  NUMBTR  OF  TRANSFORM  VALJES 
MUST  BE  THE  SAME. 

2.  BOTH  THE  INPUT  DATA  ANO  THE  TRANSFORM  VALUES  MUST  REPRPSFNT 
fQUISPACED  POINTS  IN  THEIR  RESPECTIVE  OOMAINS  OF  TIME  ANO 
FREQUENCY.  CALLING  THESE  SPACINGS  OElTAT  ANO  OELTAF,  IT  MUST  BE 
TRUE  THAT  OELTAF=2*PI/(NN(I)*QELTAT> . OF  COURSE,  OElTAT  NEED  NOT 
ME  THE  SAME  FOR  EVERY  DIMENSION. 

3.  CONCEPTUALLY  AT  LEAST,  THE  INPUT  nATA  ANO  THE  TRANSFORM  OUTPUT 
REPRESENT  SINGLE  CYCLES  OF  PERIOOIC  FUNCTIONS. 

EXAMPLE  1.  THREE-DIMENSIONAL  FORWARD  FOURIER  TRANSFORM  OF  A 
COMPLEX  ARRAY  OIMFNSIONfO  32  BY  25  BY  13  IN  FORTRAN  IV. 

DIMENSION  DATA  (3?,  2r  ,13)  , WORK  (5  0)  ,NN(3) 

COMPLEX  OATA 
OATA  NN>32,25,13/ 

00  1 1*1,32 
00  1 J *1 , 25 
00  1 K*l, 1 3 

OATA(I ,J,K) *C0MPIEX  VALUE 
CALL  FOURT(OATA,NN, 3,-1, 1, WORK) 

EXAMPLE  2.  ONE-DIMENSIONAL  FORWARO  TRANSFORM  OF  A RFAL  ARRAY  OF 
LFNGTH  54  IN  FORTRAN  IT. 

DIMENSION  OATA  (2,64) 

00  2 1*1,64 

DAT  A (1 ,1 ) =RE  AL  PART 

PATA(2,I)*0. 

CALL  F0URT(DATA,64, 1,-1, 0,0) 

DIMENSION  OATA (1) , NN ( 1) , I FACT ( 32) ,WDRK(1> 

TWOPI*fe.  2*3135307 
NR*0. 

MI*  0 • 

MSTPI*  0. 

WSTPR*  0. 

IF  (NOI  M- 1)920,1,  1 
NT0T*2 

DO  2 1 01  M*  1 » NOIM 
IF (NN( IOIM) ) 920,920,2 
NTOT»NTOT*  NN(IOIM) 

MAIN  LOOP  FOR  EACH  DIMENSION 

NP1*2 

DO  910  IDIM*1,NP7M 
NsNN(IOIM) 

NP2«NP1»N 

IF (N- 1)920,900, 5 

FACTOR  N 


‘ 31 
32 

1, 


to 

50 

51 


t 
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l 
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60 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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70 


71 

72 

73 
7 t 


i 

i 


c 

c 

c 


95 

C 

C 

c 

c 

100 

110 


• * 1 20 

C 

1 

t 


IREM*M-IOIV*IQUOT 

IF(IQUOT-IOIV)  60 , 31, 31 

IF(IREM)  tO,  32,  tO 

1FACT(IF)*I0IV 

IF«IF*i 

M«IQUOT 

GO  TO  30 

I0IV*I0IV*2 

GO  TO  30 

IF  { IRE  M)  60,51,60 
NTWO*NTWO*NTWO 
GO  TO  70 
IFACT(IF)*M 

SEPARATE  FOUR  CASES— 

1.  COMPLEX  TRANSFORM  OR  REAL  TRANSFORM  FOR  THE  • IH,  E TH, f f C. 
DIMENSIONS. 

2.  REAL  TRANSFORM  FOR  THE  2N0  OR  3RD  DIMENSION.  METHOD- 
TRANSFORM  HALF  THE  DATA,  SUPPLYING  THE  OTHER  HALF  9Y  CON- 
JUGATE SYMMETRY. 

3.  REAL  TRANSFORM  FOR  THE  1ST  DIMENSION,  N OOD.  METHOD- 
TRANSFORM  HALF  THE  DATA  AT  EACH  STAGE,  SUPPLYING  THE  OTHER 
HALF  BY  CONJUGATE  SYMMETRY. 

t.  REAL  TRANSFORM  FOR  THE  1ST  DIMENSION,  N EVEN.  METHOD- 
TRANSFORM  A COMPLEX  ARRAY  OF  LENGTH  N/2  WHOSE  REAL  PARTS 
ARE  THE  EVEN  NUMBERED  REAL  VALUES  AND  WHOSE  IMAGINARY  3 ARTS 
ARE  THE  OOn  NUHPEREO  REAL  VALUES.  SEPARATE  AND  SUPPLY 
THE  SECOND  HALE  BY  CONJUGATE  SYMMETRY . 

NON2«NPl*  (NP2/NTWO) 

ICASE* 1 

IF  CIOI  M-4 >71 ,90,90 
IFEIFORMW  2,72,90 
ICASE*  2 

IF(IDI  M-i)  73,73,90 
ICASE*  3 

IF(NTWO-NP1)90  ,9f  ,7t 

TCASE*t 

NTWO*NTWO/2 

N=N/2 

NP2*NP2/2 

NT0T*NTDT/2 

1*3 

00  90  J*  2,  NTOT 
OAT  A ( J )*OAT  A(I) 

1*1*2 
I1RNG*  NP1 

IF  (ICASE -2)  100  ,9 5*100 
I1RNG*NP0*  (l*NPREV/2) 

SHUFFLE  ON  THE  FACTORS  OF  TWO  IN  N.  AS  THE  SHUFFLING 

CAN  BE  DONE  BY  SIMPLE  INTERCHANGE,  NO  WORKING  ARRAY  IS  NEEDED 

IF (NTWO-NF1)  600,600,110 

NP2HF*NP2/2 

J*1 

no  150  1 2* 1 » NP2, N0N2 
IF(  J-I  2)  120,130,  130 
IlMAX*I2*NON2-2 
00  125  11*12, UMAX, 2 
DO  125  13*11, NTOT, NP2 
U3*  J*I  3-12 
TFMPR*OATA  (13) 

TEMPI* DAT A (13*1) 

DATA  ( 1 3)  *OAT  A ( J3  ) 


OAT  A ( J 3)  =TEMPR 
125  DATA  ( J3M)  *T  EMPI 

130  M*NP2HF 
140  rF(J-N)15C, 150,145 

145  J*J-M 
M*M/2 

rF(M-NOM2)  15  0,14  0,14  0 
150  U=J*M 
C 

C MAIN  LOOP  FOR  FACTORS  OF  TWO.  PERFORM  FOURIER  TRANSFORMS  OF 

C LENGTH  FOUR,  WITH  ONE  OF  LENGTH  TWO  IF  NEF OEO.  THC  THIDOLE  -SCTOR 

C M* EXP( IS IGN# 2* PI  SQRT(-1)*M/(4*MMAX> ).  CHECK  FOR  W=ISI  GN*S Q*T ( -1) 

C AND  REPEAT  FOP  W= ISIGN*SQRT ( -1) *CON JUG AT F < M) . 

C 

N0N2T=N0N2*N0N2 

IPAR*NTHO/NPl 

310  IF(IPAR-2)  350,330,320 
320  IPARsIPAR/4 
GO  TO  310 

330  OO  340  11=1,  I1RNG,  2 

00  340  J3=I1,N0N2,NP1 
00  340  <1=J3,NT0T,N0N2T 
K*2=K1*N0N2 
TEMPRsOATA  <K2) 

TEMPI=OATA(K2M> 

DATA(K2)  =0  AT  A ( K1 ) -T  EMPR 
DATA (K  2*1) =D AT  A ( K1 M> -TEMPI 
OATA(Kl)  =PATA(K1)  ♦TEMPR 
340  0ATA(K1*1)=DATA(K1>1) ♦TEMPI 

350  MMAX*NOM2 

360  IF  ( MMA  X-NP2HF)  370, 600,600 
370  LMAX*MAX0 (N0N2T, MMAX/2) 

IF  ( MMA  X-NON  2)  4 05  ,4  05,38C 
380  THETA=-TMOPI*FLOAT (N0N2) /FLO AT (4*MM8 X) 

IF  (ISI  GN)  4 0 0 , 390 , 39  0 
390  THETA=-THETA 
400  NR=C0S (THETA) 

WI=SIN (THET  A) 

HSTPR* -2.»  HI*MI 
MSTPI=2,*NR*NI 

405  00  570  L=NON2,LMAX, NON2T 

M=L 

IF (MMA X-NON2) 4 20, 420,410 
410  W2R*WR*HR-WI*MI 
*121  = 2. *HP*WI 
N3R=V42R*  WR-W2I*WI 
*»3I=M2P*MI*H2I*WP 
420  00  530  11=1, IlRNG, 2 

00  530  J3=I1,N0N2,NP1 
KMIN=J3*IPAR*M 
IF ( MMA  X-NON2)430, 430,440 
4 30  KMIN*  J 3 

440  KOIF=IP#R*MMAX 
450  KSTEP=4»K0IF 

00  520  <1=KMIN,NT0T,KSTEP 

K2«K1*K0IF 

K3=K2*K0IF 

K4*K3*KDIF 

IF (MMA X-N0N2) 460 ,460,460 
460  U1R*0ATA (Kl) +OATA (K2) 

U1I*0ATA ( Kl+1) ♦OATA(K2*l) 

U2R*OATA(K3) ♦0ATA(K4) 

U2I*OATA (K3  + 1) +0 AT  A (K4*l> 

U3R»OATA  (Kl)  -DATA  (K2) 

U3I*DATA ( K 1 ♦ 1 > -DATA(K2*1) 


470  O4R«0ATA  <K3*1>  -0ATA(K4*1> 

U4  I *DA  Tl  (K4>  -DATA  (K3> 

GO  TO  510 

475  U4R*0ATA<K44l)  -OATA(K3*l> 

U4I»0ATACK3)-0ATA<K4) 

GO  TO  510 

400  T2R*W2R*0ATA<K2>  -W2I* DAT  A <K2*1> 

T2I*W2R*  DAT  A (K  2* 1) ♦W21*OATA ( K2) 

T 3R*WR*D AT  A < K3)*WI*  DAT  A < K3+1 ) 

T3I*WR*DATA  0(3  *1 ) ♦MI*  DATA (K3 ) 

T4R*W3R*DATA(K4>  -W3I»0AT A (Kh *1) 

T4I«W3R»0ATA(K4*1> ♦M3I*DATA (K4) 

U1R«0ATA  (Kl)  4T2R 
U1I»0ATA  (K 1 ♦ 1)  4T2I 
U2R»T3R*T<*  R 
U2I*T3I*T4I 
U3R*0ATA<K1> -T 2R 
U3I*OATA  <Xl4l>  -T2I 
IF  ( I SI  GNH-90  ,500  ,500 
490  U4R«T3I-T4I 
U4I-T4R-T3R 
GO  TO  510 
500  U4R»T4I-T3I 
U4I*T3R-T4R 

510  OATA(Kl)  «U1R*U2R 

DATA <K 1*1 ) *U1I*U2I 
DAT  A CK  2)  =U3R*U4R 
OA  T A ( K 2*  1 ) »U  31  HM>  I 
DATA  (K  3)  «U1R-U2R 
DATA<K3*1)*U1I-U2I 
DAT  A(K4) «U3R-U4R 
520  0ATA1K4M)  »U3I-U‘ I 
fMIN»4*<KMlN-J3>  ♦ J 3 
KOIF*KSTEP 

IF ( KOI  F-NP2 ) 45  0, r 30,530 
530  CONTINUE 
M*MMAX-M 

IF  (I SI  GN)R40  ,550  ,550 
54  0 TEMPR*  W^ 

WR*-MI 
NI*-TEMPR 
GO  TO  550 
550  TEMPR*  MR 
*R«WI 
*I»TEMPR 

560  IF  (M-L  MAX)  565, 565, 410 
565  TEMPR*  WR 

WR*WR*WSTPR-WI*WSTPI+WR 
570  NI*WI* W$TPR*TEMPP*WSTPI*W1 
IPAR*3-IPAR 
HMAX«MMAX*MMAX 
GO  TO  350 
C 

C MAIN  LOOP  FOR  FACTORS  NOT  EQUAL  TO  TWO.  APPLY  THE  TWIDDLE  FACTOR 

C W»EXP<  ISIGN*2*PI*SQRT  <-l>*  ( J2-1)  M Jl -J2)  / (NP2*  IFPl)  > , THEN 

C PERFORM  A FOURIEP  TRANSFORM  OF  LENGTH  IFACT(IF),  MAKING  USE  DF 

C CONJUGATE  SYMMETRIES. 

C 

600  IF(NTWO-NP2)605,700,7u0 

605  IFPl»NON2 
IF*1 

NPlHF*NPl/2 

610  IFP2*IFM/IFACTCIF> 

J1RNG*NP? 

If  (ICASE-3)  612,611,612 


V A*  \ * <1  w » * « « *#  # i. 

J2STP*NP2/IFftCT(IP) 

J1RG2=  (J2STP*IFP2)  /2 
612  J2MIN*1»IFP2 

IF (IFP1-NP2) 615,64 0,6^0 
615  00  635  J2*J2HIN,IFP1,IFF? 

THETA* -T H0PI*FL0AT ( J2-1) /FLOAT (NP2) 

IF  ( ISI GN)  b 25  ,6  20  ,6  2C 
620  THETA*-THFTA 
625  SIMTH*  SI  N ( THET  A/  2 . ) 

WSTPR*  -2«*SINTH*FINTH 
«STPI=SIN(THETA) 

WR=WST  PR^l. 

WI*WSTPI 

J1MIN=J2*IFP1 

00  635  J1=J1MIN, JlRNGjIFPl 

11HAX*  J1  ♦IlRNG-2 

OO  630  I1  = J1,I1MAX,  2 

DO  630  13=11, NTOT ,NP2 

J3MAX=I3*IFP2-NP1 

00  630  J3=I3, J3MAX,NP1 

TENPR= DATA ( J3) 

DAT  A ( J 3)  =0ATA(J3)*MF-DATA(J3*1)*M1 
6 30  OATA  ( J3*l)  =T  EMPR*  VI  ♦DATA  ( J3*1)*WR 
TEMPR= MR 

HRs HR* WSTPR-MI*MSTPI*WR 
635  WI  = TEHPR*HSTPI*WI*HSTPR+W1 
640  THETA=-TWCPI/FLOAT  (IFACT  (IF)  ) 
IF(ISIGN)650,645 ,645 
645  THETA*-THFTA 
650  SINTH=SIN(THETA/2.) 

MSTPR=  -2.*  SI  NTH*  SI  NTH 
MSTPI=  SI  N ( THET  A) 

*STEP=2*N/IFACT(IF) 

KRANG=KSTEP* (IFACT (IF)/?) >1 
00  698  I1=1,I1RNG,2 
00  698  13=11, NTOT, NP2 
00  690  KMIN=1,KRANG,KSTEP 
J1MAX=I3*J1RNG-IFP1 
00  680  J1=I3,J1MAX,IFP1 
J3NAX= J1 ♦IFP2-NP1 
OO  680  J3=  J1 , J3MAX, NPl 
J2MAX= J3*IFP1-IFP2 

K=KMIN ♦( J3-J1+ (J 1-13) /IFACT (IF) )/N°) HF 
IF(KMI N-l) b?5,65r ,665 
655  S(JMR=Q  • 

SUNI=0 . 

00  660  J2=J3,J2MAX,IFP2 
SUMR=S  UMR+DATA (J2) 

660  SUNI=SUMI+OATA (J2+1) 

M0RK(K)=SUMR 
M0RK(K«-1)  =SUMI 
GO  TO  630 

665  KC0NJ=K*2*  (N-KMIN  + 1) 

J2= J2M  AX 
SUMR=D AT  A ( J2) 

SUMI=OATA(  J2  + 1) 

0L0SR=  0, 

OLOSI=0. 

U2* J2-IFP? 

670  TEHPR*  SUMR 
TEMPI* SUMI 

SUMR*T  WOMR* SUMR-OL DSR  +DAT  A ( J2) 

SUMI*TMOMR*SUHI-OLOSI+OATA(J2*1) 

OL0SR*TEMPR 

OUOSI-IIMPI  _ _ 
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IF(  J2-  J3)  675,675  ,67P 
675  TEMPR=M**SUMR-OLOSR*OATA (J2> 

TEMPI* MI* SUMI 
WORKOO*TFMPR-TFMPI 
NORK(KCONJ)  = TEMPR*TEMPI 
TEMPR*WR*SUMI-OLOSI*OAT  (J2+1) 

TEMPI* MI *SUMR 
N0RK(K«-1)=TEMPR*TEMPI 
MORK(KCDNJM)  *TEMPR-TFMPI 
680  CONTINUE 

IF ( KMI N“l) 6 8 5 » 66  F , 6 66 
685  MR*MSTPRM. 

MI*HSTPI 
GO  TO  69  0 
6 86  TEMPR*MR 

WR= MR*  MST°R-MI*WFTPI*WR 
NI*TEMPR*MSTPI+MT*WSTPR»MI 

690  TMOMR=MR*MR 
TF(ICASE-3>692,691,692 

691  IF(IFPl-NP2)695,f  92,692 

692  K*1 
I2MAX=I3+NP2-NP1 

OO  693  12=13, I2MAX,NP1 
OAT  A (I  2)  *MORK ( K) 

OAT  A (I 2*i) =WORK(K*l) 

693  K=K*2 

GO  TO  698 
C 

C COMPLETE  A REAL  TRANSFORM  IN  THE  1ST  DIMENSION, 

C JUGATE  SYMMETRIES  AT  EACH  STAGE. 

C 

695  J3HAX=I3*IFP2-NP1 

DO  697  J3=I3, J3MAX,NP1 
J2MAX*  J3  6NP2-J2STP 
TO  697  J2=J3, J2MAX, J2STP 
J1MAX=  J2  + J1RG2-IFP2 
U1CNJ*J3*J2MAX*J2STP-J2 
OO  697  J1=J2,J1MAX,IFP2 
K=1+J1-I3 
nATA(Jl)=WORK(K) 

DATA(J1*1)  =MORK(K4-l> 

IFCJ1-J?)  6 9 7,697  , 6 9 6 

696  OATAtJIDNJ)  =MORK(K) 

DATA  (J1DNJ  + 1)  = -WORK  (K  + l) 

697  U1CNJ=  J1CNJ-IFP2 

698  CONTINUE 
IF=IF+1 
IFPi=I F°2 

IF (IFP1-NP1)  70  0,7  0 0,610 
C 

C COMPLETE  A REAL  TRANSFOFM  IN  THE  1S1  DIMENSION, 

C JUGATE  SYMMETRIES. 

C 

700  GO  TO  (900,800  ,900,701)  ,ICASE 

781  NHALF*  N 

N=N+N 

THETA* -TMOPI /FLO AT (N> 

IF(ISIGN)703,7  02,?02 

702  THETA* -THET  A 

703  SINTH*SI  N ( THET  A/  2 , ) 

HSTPR*  -2.-  SINTH*  SINTH 
MSTPI*SI  N (THET  A) 

MR*MST PR+1. 

MI*HST  PI 
IMIN*3 

^ 


N OJD,  BY  2""i- 


N EVEN, 


BY 


D N 


GO  TO  725 
710  J= JMIN 

no  720  I=IMIN, NTOT,NP? 

SUMR=<DATA<I)+0ATA(J>)/2. 

SUMI= ( OATA(I+l)+DATA(J+l))/2. 

OIFR=<DATA (I)-OATA  (J)  )/2. 

OIFI=  (DATA  <I+l)-OATA<J+l>>/2. 

TEMPR=  WR+SUMI  + WI*  DIFR 
TEMPI=WI  + SUMI-WR*DIFR 
DATA (I ) = SUMR+T EMPR 
OATA(I+l)=DIFI+TEMPI 
OAT  A ( J ) = SUMR-T  EMPR 
OATA(J+l)=-DIFI+TEMPI 
720  J=J+NP2 

IMIN=I  MIN+2 
JMIN=JMIN-2 
TEMPR=  M* 

MR= HR*  WSTPR-WI+WSTPI+WR 
MI*TEMPR*  WSTPI +MI+  WSTPR  + WI 
725  IF(IMIN-JMIN)710  ,730,740 

730  IF(ISIGN>731,740,740 

731  OO  735  I=IMIN,NTOT,NP2 
735  OATA(I+l)=-DATA<I  + l) 

740  NP2=NP2+NP2 

NTOT=NTDT+NTOT 
U=NTOT  +1 
IMAX=NTDT/2+l 
745  IMIN=I  MAX-2+NHAL  F 
Is I MIN 
GO  TO  755 

750  OATA(J)  = OATA  (I) 

OATA(J+1)=-OATA(I+1) 

755  1=1+2 

J= J -2 

rF(I-IMAX)750, 76  0,760 
760  OATA(J)  = OATA (I  MIN) -DATA (IMIN+1) 

OATA(J+1)=0. 

IF(I-J >770,780,780 
765  OATA(J)sOATA (I) 

OATA(J+l)=OATA  (1+1) 

770  1=1-2 

J=  J-2 

If  ( I-I  MIN)775,77f,765 
7 75  OATA(J)  = DATA(IMIN) +DAT A ( I MIN +1) 

DAT  A ( J +1 ) = 0 • 

IMAX  = I MI  N 
GO  TO  7*5 

780  OAT  A(l)  = OATA (1)+DATA(2) 

OAT  A (2 ) a 0 • 

GO  TO  900 
C 

C COMPLETE  A REAL  TRANSFORM  FOR  THE  2ND  OR  3RD  DIMENSION  BY 

C CONJUGATE  SYMMETRIES. 

C 

800  1F(I1RNG-NP1)605,900,900 

805  OO  860  I 3= 1 , NTOT , NP2 
I2MAX= I3+NP2-NP1 
OO  860  12=13, 12MAX,NP1 
IMIN=I 2+I1RNG 
I*AX=I2+NPl-2 
'JMAX  = 2*I  3 + NP1-IMIN 
IF  (1 2- 13)820, 8 20, 810 
810  UMA  X= JMA  X +NP2 

820  If  ( I OI  M-2  >850,850,830 

830  U=  JMAX  +NP0 

nn  auo  t tma* .z 


OATA(I  >«OATA  C J > 

OATA  (I  ♦1)"-0ATA(J*1) 


* 840  J»J-2 

850  U>JMAX 

DO  860  I«IMIN,IMAX,NP0 
DATAd )»OATA ( J> 

OATA(I  d)  *-OATA(  J*l) 

860  U* J-NPP 

C 

C FND  OF  LOOP  ON  FACH  DIMFNSION 

C 

900  NP0«NP1 
NP1»NP  2 
910  NPREV-N 
9 20  RETURN 
FNO 

SUBROUTINE  TPPMl  T < LETTEF  , TRPLTR,  SENT  NUN) 

C THIS  ROUTINE  OOES  TR APF701 O BUILD  AND  MULTIPLY 
INTEGER  LFTTFR(3?»3?) fSFNTNUM 
OIHENS ION  TRPLTR (3?, 3?) , TRPWOW (32) 

1 FORMAT (1H0) 

2 FORMAT  (1H  , 24  («•  " ) , /1H  TRAPEZOID  PARAMETERS  * THIS  TRAPEZOID 
X MAS  APPLIED  TO  WINDOWS  OF  SENTENCE  NUHBEP  ",12, /1H  ,2«("  ** ) I 

3 FORMAT  <1H  LENGTH  OF  TRAPEZOID  PLATEAU  * **,12) 

4 FORMAT  (1H  ,/,**  VARIANCE  ON  NORNALI’f D MINrOW  * M,F9.2) 

5 FORMAT  (1H  ,/,**  TPAPEZUIT  WINDOW  ARRAY  IS*  **) 

* 6 FORMAT  C"  ’M0F6.3) 

7 FORMAT (IH  THIS  IS  A NORMAL  WINDOW  APPLIED  TO  SENTENCE  NJMBER 

X-.IZ) 

IF  CPSKIP.EQ.-OR")  GO  TC  45 

COAT  A NORMALIZED  WINDOW  SET  IPLEN=0  AND  TRPHT*VARIANCE 
IPLEN  * 0 $ IRFHT  = S. 

C 

C THIS  PART  BUILDS  THE  TRAPE  701 C 
C 

INCLN* (32-IPLEN) /2 
IF  (IPLEN. NE.O)  GO  TO  9 


10 

20 


30 


35 

40 

43 


* 

L 


44 


DO  8 I >1 f 16 

TRPWOW  (17-1)  *TRPWOW(16d)*EXP(-(  (I -1  )’*?)/(  2."  ( (TRPHT) ' * 2))  ) 
WRITE  1 

WRITE  7 ( SENT NUM 
GO  TO  44 

IF  (INCLN. EQ.0)  GO  TO  3! 

DO  10  1*1  » INCLN 

TRPWOW (I ) * TRPWOW (33-1) * ( I / ( I NCLNfl . ) ) * TRPMT 

CONTINUE 

PO  20  1*1,  IPLEN 

TRPWOW (T*INCLN )® 1 • 

CONTINUE 

IF((32-(2'  INCLN)  ) .E0.1P(  EN)  GO  TO  43 

iooo«in:ln*i 
00  30  I«ltIOOO 

TRPWOW (33-1) * ( 1/ (INCLN*?. )>* TRPHT 
CONTINUE 
GO  TO  43 
DO  40  I* 1 « 16 

TRPWOW (I) « TRPWOW ( 33-1) *1 . 

CONTINUE 
WRITE  1 
WRITE  1 

WRITE  ?t  SFNTNUM 
WRITE  3,  IPLEN 
WRITE  4.TRPHT 
WRITE  5 

WRITE  6, TRPWOW 


o o 


THIS  PART  ODES  THF  MULTIPLICATION 

kb  DO  60  1*1,32 
00  50  J*  1 » 3 2 

50  TRPLTR (J,I> *TRPWOW(I>* (FL0A1  (LETT FR ( J, I ) ) ) 
60  CONTINUE 
PSKIP«"OKM 
RETURN 
END 


ROYOOSV  ////  END  OF  LIST  //// 
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